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Abstract 



This thesis presents the results of a high spectral resolution (/l/A/l = 49000) study of the 
circumstellar (CS) gas around the intermediate mass, pre-main sequence UXOR stars 
BF Ori, SV Cep, UX Ori, WW Vul and XY Per. The resuhs are based on a set of 38 
echelle spectra covering the spectral range 3800-5900 A, monitoring the stars on time 
scales of months, days and hours. 

All spectra show a large number of Balmer and metallic lines with variable blueshif- 
ted and redshifted absorption features superimposed to the photospheric stellar spectra. 
Synthetic Kurucz models are used to estimate rotational velocities, effective tempera- 
tures and gravities of the stars. The best photospheric models are subtracted from each 
observed spectrum to determine the variable absorption features due to the circumstellar 
gas; those features are characterized, via multigaussian fitting, in terms of their velocity, 
V, dispersion velocity, Av, and residual absorption, i?max- 

The absorption components detected in each spectrum can be grouped by their sim- 
ilar radial velocities and are interpreted as the signature of the dynamical evolution of 
gaseous clumps. Most of the events undergo accelerations/decelerations at a rate of 
tenths of m s~^. The typical timescale for the duration of the events is a few days. The 
dispersion velocity and the relative absorption strength of the features do not show dras- 
tic changes during the lifetime of the events, which suggests that they are gaseous blobs 
preserving their geometrical and physical identity. 

A comparison of the intensity ratios among the transient absorptions suggests a solar- 
like composition for most of the CS gas. This confirms previous results and excludes a 
very metal-rich environment as the general cause of the transient features in UXOR stars. 
These data are a very useful tool for constraining and validating theoretical models of 
the chemical and physical conditions of the CS gas around young stars; in particular, it 
is suggested that the simultaneous presence of infalling and outflowing gas should be 
investigated in the context of detailed magnetospheric accretion models, similar to those 
proposed for the lower mass T Tauri stars. 

WW Vul is unusual because, in addition to infalling and outflowing gas with proper- 
ties similar to those observed in the other stars, it shows also transient absorption features 
in metallic lines with no obvious counterparts in the hydrogen lines. This could, in prin- 
ciple, suggest the presence of CS gas clouds with enhanced metallicity around WW Vul. 
The existence of such a metal-rich gas component, however, needs to be confirmed by 
further observations and a more quantitative analysis. 

All these results have been published by Mora et al. (2002, Chapter 5) and Mora 
et al. (2004, Chapter 6). 

Rotational velocities for a large sample of pre-main sequence and Vega-type stars 
have been determined from high resolution (/l/AA = 49000) echelle spectra. The first 
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Abstract 



minimum of the Fourier transform of many photospheric line profiles has been used to 
estimate the velocities. The resulting velocities have been published, along with spectral 
types determined by other co-authors, by Mora et al. (2001, see Chapter 4). 

Key words. Stars: formation - Stars: pre-main sequence - Stars: circumstellar matter 
- Accretion: accretion disks - Lines: profiles - Stars: individual: BF Ori, SV Cep, 
UX Ori, WW Vul, XY Per - Stars: rotation 



Resumen 



En esta tesis se presentan los resultados de un estudio, basado en espectros de alta reso- 
lucion (/l/A/l = 49000), del gas circunestelar (CircumStellar, CS) en las estrellas UXOR 
de masa intermedia BF Ori, SV Cep, UX Ori, WW Vul y XY Per. Las observaciones 
empleadas son un conjunto de 38 espectros echelle, obtenidos en el rango espectral 3800- 
5900 A, con los que se ha efectuado un seguimiento de las estrellas en escalas temporales 
de meses, dias y boras. 

Todos los espectros muestran un gran numero de componentes de absorcion des- 
plazadas al rojo y al azul superpuestas sobre el espectro fotosferico estelar, en Imeas 
de Balmer y metalicas. Se han utilizado espectros sinteticos, generados a partir de los 
programas y modelos de atmosfera de Kurucz, para estimar velocidades de rotacion, 
temperaturas efectivas y gravedades estelares. Se ban sustraido los mejores modelos 
fotosfericos de cada espectro observado para determinar las componentes de absorcion 
variables debidas al gas circunestelar. Dichas componentes ban sido caracterizadas, me- 
diante ajustes gaussianos multicomponente, en terminos de su velocidad, v, dispersion 
de velocidades, Av, y absorcion residual, ^max- 

Las componentes de absorcion detectadas en cada espectro se pueden agrupar en 
eventos de acuerdo a la similitud de velocidades radiales, lo cual es interpretado como 
la buella de la evolucion dinamica de condensaciones de gas. La mayorfa de los eventos 
experimentan aceleraciones/desaceleraciones del orden de decimas de ms"^. El tiempo 
de vida tfpico de estos eventos es de unos pocos dias. Ni la dispersion de velocidades ni 
la intensidad relativa de las absorciones de las componentes muestran cambios drasticos 
durante el tiempo de vida de los eventos. Esto sugiere que son originados por pequenas 
nubes de gas que mantienen su identidad geometrica y ffsica. 

El estudio de las relaciones de intensidad entre distintas Imeas para los distintos 
eventos sugiere una composicion similar a la solar para la mayor parte del gas CS. Esto 
confirma algunos resultados previos y excluye un medio ambiente muy rico en metales 
como el origen comun de las componentes transitorias en estrellas UXOR. Los datos 
obtenidos suponen restricciones o validaciones observacionales que pueden ser apli- 
cadas en los modelos teoricos de condiciones ffsico-qufmicas del gas CS en estrellas 
jdvenes. En particular, se sugiere que la presencia simultanea de gas en caida y eyectado 
por la estrella deberia ser investigado en el contexto de modelos detallados de acrecion 
magnetosf erica, similares a los propuestos para las estrellas T Tauri de baja masa. 

Se ha descubierto que WW Vul es una estrella peculiar porque, ademas de mostrar 
gas en cafda y eyeccion con propiedades similares a las observadas en las otras estrellas, 
presenta tambien componentes de absorcion transitorias en Imeas metalicas sin contra- 
partida obvia en las Imeas de hidrogeno. Este hecho podria, en principio, sugerir la 
presencia de nubes de gas CS de metalicidad elevada alrededor de WW Vul. La exis- 
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Resumen 



tencia de una componente de gas rico en metales debe, sin embargo, ser confirmada 
mediante observaciones adicionales y un analisis cuantitativo detallado. 

Todos estos resultados ban sido publicados por Mora et al. (2002, capitulo 5) y Mora 
et al. (2004, capitulo 6). 

Tambien se ban determinado velocidades de rotacion para una gran muestra de estre- 
llas PMS y de tipo Vega mediante espectros echelle de alta resolucion (/l/A/l = 49000). 
Para determinar las velocidades se ha utilizado el primer mmimo de la transformada de 
Fourier de diversos perfiles de Imea fotosfericos. Las velocidades obtenidas ban sido 
publicadas por Mora etal. (2001, ver capitulo 4), junto con determinaciones de tipos 
espectrales realizadas por otros coautores. 

Palabras clave. Estrellas: formacion - Estrellas: pre-secuencia principal - Estrellas: 
materia circunestelar - Acrecion: discos de acrecion - Lmeas: perfiles - Estrellas: indi- 
viduales: BF Ori, SV Cep, UX Ori, WW Vul, XY Per - Estrellas: rotacion 
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Chapter 1 
Introduction 



This thesis, entitled "Kinematics of the circumstellar gas around UXOR stars", has been 
developed in the scientific framework of the evolution of circumstellar (CS) disks around 
Pre-main sequence (PMS) stars. Its specific topic is the detection and characterization of 
Transient Absorption Components (TACs) in high resolution spectra of five intermediate 
mass PMS Herbig Ae/Be (HAeBe) stars. These objects belong to the HAe subgroup of 
stars, which comprises HAeBe stars with masses lower than 5 M©. The stars in this thesis 
have photopolarimetric behaviour similar to the HAe star UX Orionis (included in the 
sample), therefore they are called UXORs. It is believed that the UXORs are surrounded 
by large CS protoplanetary disks seen edge-on. 

The current status of the research projects most relevant for this thesis is reviewed in 
Sections 1.1 to 1.6. The evolution of the research carried out in this thesis is exposed in 
Sections 1.7 and 1.8. Finally, the structure of this dissertation is shown in Section 1.9. 

1.1 Discovery of the first extrasolar planets 

The search for extrasolar planets is a very interesting research field, both from the point 
of view of fundamental astrophysics (formation and evolution of stars and planetary 
systems) and the philosophical viewpoint (possible existence of extraterrestrial life). 

The first detection of extrasolar planets was made by Wolszczan & Frail (1992), 
who discovered three bodies orbiting around the pulsar PSR1257+12 via very precise 
measurements of the emission time of the pulses. Two of these planets are similar to the 
Earth in mass and orbital period. The third body is somewhat smaller. However, Mayor 
& Queloz (1995) were the first to discover an extrasolar planet around a Main Sequence 
(MS) star: 5 1 Peg. Their method relies in precise measurements of the radial velocity 
of the stars with accuracies down to ~10 ms~^ This work has become a fundamental 
milestone in the star and planetary formation field. Many exoplanets around MS stars 
have been discovered so far with the radial velocities technique, at least 120^ at the date 
of writing up (April 6* 2004). 

As soon as the number of discovered planets around MS stars was significant, it be- 
came apparent that the "classic" paradigm of planetary system formation was inadequate 

'An updated list of all discovered exoplanets can be found in the "Extrasolar Planets Encyclopaedia", 
http : //www . obspm . f r/encycl/encycl . html 
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to describe the observations. This should not be surprising, as the theory was developed 
from only one particular case: the Solar System. The main lack of the theory is its inabil- 
ity to explain the common appearance of giant planets very close to the star, even nearer 
than the orbital distance of Mercury. These bodies, which are called "Hot-Jupiters", are 
the easiest to detect via the radial velocities method. 

1.2 Evolution of circumstellar disks 

The theories of the formation of planetary systems have been very much improved since 
the discovery of the first extrasolar planets. However, there is not yet a theory which 
explains, in a self-consistent way, the physical mechanisms involved in the generation 
of the observed large diversity of planetary systems. It is also believed that the theories 
will be reformed again when methods capable of detecting telluric planets at distances 
~1 AU are developed. 

However, there is agreement among the theoreticians that the planets are formed in 
the CS accretion disks (Ruden, 1999). These disks are originated during the gravitational 
collapse undergone by a molecular cloud core prior to the generation of a star. The 
angular momentum accumulated by a cloud core is very large because of its very big 
initial dimensions (i?core ~ 0.1 pc). In this way, it is impossible for the cloud material to 
fall directly onto the centre of gravity. The matter is accumulated in a disk, which rotates 
around the central condensation (protostar), where the stellar mass is being accreted 
(Hartmann, 1998). The evolution of a disk is regulated by the matter that arrives at its 
outer edge from the collapsing envelope and several internal processes (mainly viscous 
evolution and gravitational instabilities). All these processes imply a parameter, namely 
mass accretion rate, M, from the disk to the central protostar. 

The evolution of protoplanetary disks in low mass T Tauri stars is rather well un- 
derstood (Hartmann, 1998). The observations suggest that, for a Classical T Tauri Star 
(CTTS) with a typical mass of IMq, the accretion rate from the envelope to the disk is 
about 10~^Mo/yr at the beginning of the core collapse. This rate steadily decreases until, 
eventually, it vanishes after 0. 1-0.2 Myr. However, the observed accretion rates, M, from 
the disk to the protostar are two orders of magnitude lower. Hartmann (1998) suggests 
that most of the matter is transferred to the protostar during some FU Orionis outbursts 
(he admits, however, that the number of detected FU Ori objects is not large enough, 
compared to the number of identified PMS stars). These phenomena are characterized 
by a huge increase in the stellar luminosity during brief lapses of time (? ~ 100 yr) and 
are generally associated with dramatic changes in the accretion rate M. According to 
Hartmann (1998), if the accretion rate from the envelope to the disk is higher than that 
from the disk to the protostar, the disk mass grows steadily. Disks may become gravita- 
tionally unstable if their mass reaches a certain threshold value. These instabilities could 
allow large amounts of matter to fall quickly onto the star. The potential energy released 
would give rise to the outburst, characterized by a large increase in brightness. 

When the envelope eventually collapses, the result is a CTTS surrounded by a CS 
disk, whose evolution is essentially dominated by viscous friction processes (Hartmann, 
1998). In the CTTS stage, the accretion rate M experiments a much smoother evolution: 
it decreases steadily as the disk is being depleted of matter (there is room for possible 
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Age (yr) 

Figure 1.1: Temporal evolution of the accretion rate in a typical CTTS (taken from 
Calvet et al., 2000). Three different regimes can be identified in the figure. The first 
one extends from to 0.1-0.2 Myr and corresponds to the gravitational collapse of the 
molecular cloud core. It can be appreciated that the accretion rate from the envelope to 
the disk (dashed line) is higher than the accretion rate, M, from the disk to the protostar. 
This difference of material generates sporadic gravitational instabilities in the disk that, 
eventually, give rise to FU Ori outbursts. During the second stage, which spans from 
0.1-0.2 Myr to a few Myr, the CS disk gives matter to the star and expands via viscous 
friction processes. Finally, after a few Myr, the accretion disk dissipates and no more 
matter is transferred to the star. 



EXor outbursts, which are less violent than FU Ori ones). Even though the disk becomes 
less massive, its size increases due to angular momentum transport in the outward direc- 
tion. This scenario is only valid in isolated stars. If the star is in a binary system, the 
companion would stop the disk expansion, and the disk would be consumed in a much 
shorter time. It is thought (Calvet et al., 2000) that this is the origin of Weak T Tauri 
Stars (WTTSs), which have similar ages than CTTSs but lack traces of accretion or CS 
disks. 

A schematic representation of the temporal evolution of the mass accretion rate for 
a typical CTTS is shown in Figure 1.1, taken from Calvet et al. (2000). Three different 
regimes can be identified in the figure: molecular cloud core collapse, viscous accretion 
of the disk and disk dissipation. 

The protoplanetary disks around HAeBe stars are much less known. Moreover, there 
are differences between the stars of spectral types B9 and later (HAe) and those of earlier 
types (HBe). In general, most of the HAe stars have disks similar to those of CTTSs 
during the whole PMS phase and up to about 10 Myr (Natta et al., 2000a), when second 
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generation disks (see below) begin to be detected. On the other hand, the HBe stars 
usually lack protoplanetary disks. This has been interpreted by Natta et al. (2000a) as a 
consequence of the intense radiation fields present in the neighbourhood of these objects. 
The radiation could destroy the disk before the primordial molecular cloud remnants are 
dissipated and the star becomes detectable. 

1.3 Grain growth and planetesimals 

As we said in the previous section, it is believed that the formation of planets takes place 
in protoplanetary disks (Ruden, 1999). The whole process of planetary formation can 
be divided into three general stages (Beckwith et al., 2000): 1. The dust grains at the 
disk, which have initial sizes typical of the interstellar (IS) medium (J ~ 1 yum), grow 
up to sizes d ~ I km in about ? ~ 10"* yr. 2. The solid bodies of J ~ 1 km, called 
planetesimals, grow in a very fast runaway process until the telluric planets and giant 
planet rocky cores are formed. 3. The rocky cores accrete the gas located in nearby 
orbits by gravitational attraction. The final mass of the planet depends on the remaining 
reservoir of gas available when the rocky core is created. In this way, telluric, ice giant 
or gas giant planets can be formed. It is expected that the whole planet formation process 
finishes after t ~ 10^ yr. Figure 1.2, (taken from Beckwith et al., 2000), shows an scheme 
of the three described phases of planet formation. 

Many numerical simulations (e.g. Weidenschilling et al., 1997; Kokubo & Ida, 2002) 
show that, once the planetesimals have been formed, the formation of a planetary system 
is almost unavoidable. The reason is that the gravitational field generated by the plan- 
etesimals (and not the dust grains) makes the captures of material much easier. Safronov 
(1969) first demonstrated that the gravity makes the coUisional cross section of a plan- 
etesimal to grow with the fourth power of the radius cr ~ r^, which is a much higher 
increment than that of the geometrical cross section, which raises with the second power 
of radius cTg ~ r^. This effect favours the captures of matter by the more massive plan- 
etesimals and originates a runaway accretion phase of the biggest objects. This stage 
will eventually finish when the planetesimals fuse to form the telluric planets and rocky 
cores of the planetary system. 

Planetesimal formation is the most critical and worst known phase during the whole 
planetary formation process. It was at first thought that these bodies are formed via 
gravitational instabilities in the disk, which could generate the gravitational collapse of 
a significant amount of mass (e.g. Goldreich & Ward, 1973). However, this hypothesis 
has been discarded because the instabilities cannot be set up until the grains have grown 
to sizes d ~ 10-100 m (Cuzzi et al., 1993; Weidenschilling, 1995). The only alternative 
mechanism for generating planetesimals is the gradual growth of dust grains by inelastic 
collisions (Beckwith et al., 2000). Initially, the dust grains have a size typical of those in 
the IS medium {d ~ \ jjm). It is thought that the grains can grow by inelastic collisions 
until they reach macroscopic dimensions of J ~ 1 m in about t ~ 10'^ yr. The settling 
of grains in the midplane of the disk is one of the most important factors in the growth, 
because it rises the frequency of collisions. The bodies continue growing until they 
become planetesimals {d ~ \ km) after t ~ \QI^ yr. 

The main problem about the planetesimals is the time needed for their formation. 
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Figure 1.2: Different phases in the process of planet formation (taken from Beckwith 
et al., 2000). It is currently believed that the planets are formed via dust grain growth 
processes inside the CS disks. At first (left), the grains, which have sizes of J ~ 1 //m 
typical of the IS medium, grow by inelastic collisions. After about t ~ 10"* yr, the grains 
reach macroscopic sizes of <i ~ 1 m. These processes continue until the grains grow up 
to J ~ 1 km after t ~ IQf^ yr. These bodies, called planetesimals, have a non-negligible 
gravitational field that accelerates their growth (centre). When a massive rocky core, 
M ~ IMgj, is formed, it will accrete all the CS gas in its surroundings. It will eventually 
become a gas giant planet if there is enough gas available. The whole process is supposed 
to finish after t ~ \ Myr. 
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because if this interval is longer than the lifetime of the gas disk, it would be impos- 
sible to generate giant planets. This has motivated several scientists to study both the 
physics involved and the observational evidences of grain growth. Many experiments 
have been performed to understand the physics of the process. Poppe et al. (2000) mea- 
sured the adhesion efficiency between micrometric particles of different materials and 
sizes. Supulver et al. (1997) studied the influence of the outer frost layers in the ad- 
hesion, which is important for particles of sizes around d ~ 1-100 cm. Kouchi et al. 
(2002) analyzed the influence of the organic surface layers in particles of size J ~ 1 mm. 
Paraskov et al. (2003) have developed an experimental setup, which will allow them to 
study the role of erosion in high velocity impacts with grains of size d ~ \ mm. In ad- 
dition, many theoretical studies have been carried out. For example, Dominik & Tielens 
(1997) calculated critical junction velocities between grains of several materials, while 
Komet et al. (2001) simulated numerically the growth of grains in the disks. 

Currently, there are not undisputable proofs of grain growth in CS disks. The most 
promising indications arise from the detailed study of the spectra and Spectral Energy 
Distributions (SEDs) in the millimetre and submillimetre ranges (Beckwith et al., 2000). 
The simplest theoretical models suggest that the accumulation of mass in big size opaque 
bodies lowers the extinction. This phenomenon must be carefully modelled in order to 
make accurate predictions, because the vast diversity of molecules present in CS disks, 
the unknown spatial distribution of the millimetre emission and the complexities of light 
reprocessing in the disk do not allow a direct interpretation of the observations. For 
example, D'Alessio et al. (2001) have predicted millimetre fluxes according to self- 
consistent models of protoplanetary disks. It is hoped that the definitive proofs will 
be obtained as soon as the millimetre interferometer ALMA becomes operational. This 
instrument will join a high sensitivity together with high spatial resolution (Beckwith 
et al., 2000). 

1.4 Search for planetesimals: p Pictoris 

It has been seen in the previous section that the computational modelling of planetesi- 
mals is feasible with the facilities already available. However, the direct observation of 
planetesimals is extremely complicated. The surface of one of such bodies is several or- 
ders of magnitude lower than the sum of the surface of all the dust grains in the disk. In 
this way, both the thermal emission and the scattered light of a planetesimal are hidden 
by the dust. Even though direct detections are impossible, Lagrange-Henri et al. (1988) 
proposed the existence of planetesimals orbiting around the star ^ Pictoris as the best 
explanation to its spectroscopic activity. 

Pic is a MS star of spectral type A5V, approximately 20+10 Myr old (Barrado y 
Navascues et al., 1999). 15 Pic is one of the prototypical Vega stars. Those stars are MS 
objects with prominent infrared excesses, named after Vega {a Lyr), the first identified 
star of the class. The Vega phenomenon was discovered after some routine calibrations 
were performed with the satellite IRAS (Aumann et al., 1984). It has been interpreted as 
the signature of possible circumstellar disks. This hypothesis was confirmed by Smith 
& Terrile (1984), who obtained a coronagraphic scattered light image of the /3 Pic disk. 
This was the first image of a CS disk, which was oriented edge-on. 
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It was soon apparent (Backman & Paresce, 1993) that the yS Pic disk is formed primar- 
ily by dust grains of size d <: ljum. Due to dissipation by collisions, radiation pressure 
and Poynting-Robertson drag, the grains have a lifetime much shorter than the stellar 
age. In this way, it should exist a source which replaces the grains as they evaporate, 
otherwise the disk could not survive for a long time. Protoplanetary disks are very dif- 
ferent from the Pic disk, because they are mainly composed of gas, they are accreted 
or dissipated in much earlier stages and they are unable to provide a continuous supply 
of dust grains. In other words, the jS Pic disk is an object of different origin, formed after 
the gravitational collapse of the molecular cloud core, so it is thought to be of "second 
generation". 

Many studies about Pic were conducted after the discovery of its CS disk. A central 
gap clear of dust, up to a radius of about 25 AU was discovered (Lagage & Pantin, 1994), 
several asymmetries in the disk, including a warp, were detected (Burrows et al., 1995). 
Amorphous and crystalline silicates have been discovered in the disk (Telesco & Knacke, 
1991; Knacke et al., 1993). Finally, small amounts of gas have been detected by means 
of Transient Absorption Components (TACs) and stable absorptions superimposed over 
several photospheric spectral lines (Hobbs et al., 1985). 

The stable components have the same radial velocity as the star and are relatively 
narrow, Av ~ 2 kms"\ (Lagrange et al., 1998), so the gas responsible for the absorption 
is placed at a constant distance from the star. This seems a paradox, because many of the 
ions showing stable components have a p rate (quotient between the radiation pressure 
and the gravitational attraction) much greater than 1. For example, yScan = 35, ySpen = 
4.87 (Lagrange et al., 1998). The same authors propose that the CS gas responsible for 
the absorption is hold by a neutral hydrogen torus (this element is little affected by the 
radiation pressure), so the radial velocity of the gas is zero. The amount of H i needed 
is small and compatible with the actual upper limits (see below). If this model is true, 
the gas causing the absorptions needs to be continually replenished, because the torus 
delays the outward migration, but does not stop it. This explanation has the additional 
property of being compatible with the theory the same authors have developed to explain 
the transient components (see below). 

Most of the absorption components (TACs) observed in p Pic are redshifted, so they 
are called RACs (Redshifted Absorption Components). There are a few detections of 
blueshifted components (Bruhweiler et al., 1991; Crawford et al., 1998), called BACs 
(Blueshifted Absorption Components). All the detections of CS gas in Pic have been 
made in metallic lines and never in hydrogen. This surprising result motivated Freudling 
et al. (1995) to search for neutral hydrogen in the p Pic disk via the 21 cm Hi line. The 
detection was negative, so they could establish an upper limit to the hydrogen column 
density of A^(H) < W^rn^. 

The RACs in yS Pic display a wide range of radial velocities, from v - 10 km s"^ to 
V <: 300 kms '. It has been observed a correlation between the radial velocity and the 
width of the RACs (Lagrange et al., 1996). In order to achieve a high radial velocity, 
the gas needs to get very close to the star (e.g. at distances of r ~ 0.1 AU for gas with 
velocity v ~ 20 kms"^) in high ellipticity orbits. However, the high radiation pressure 
exerted by [5 Pic prevents much of the observed ions to approach so near the star. In this 
way, a source that provides the gas at distances r ~ 0. 1 AU with infall radial velocities is 
needed. 
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The model developed by Lagrange-Henri et al. (1988) provides the best explanation 
of the RACs in p Pic. They assume that the RACs are generated by the evaporation of 
solid bodies (FEBs, Falling, Evaporating Bodies). These bodies would be around 1 km 
in size (i.e. they are planetesimals) and would approach the star in highly excentrical 
orbits, as the Solar System comets do. In this way, the stellar radiation would evaporate 
the outer layers of the planetesimals, thus generating a coma made by dust and gas in 
different ionization stages. The dust would continuously replenish the CS disk and the 
evaporated gas would generate the transient absorptions and the stable ones when the 
outward gas migration would be stopped by the gas torus. The tail would not have any 
influence in the RAC generation, because it would be highly coUimated and would not 
be able to obscure a significant fraction of the stellar disk. 

Beust et al. (1990) performed a series of numerical simulations, later improved by 
Beust et al. (1996) and Beust et al. (1998), which confirmed the ability of the FEB mech- 
anism to produce RACs in p Pic. An illustration of such simulations can be seen in Fig- 
ure 1.3, taken from Beust et al. (1998). The FEB model has two fundamental problems: 
first, an efficient physical mechanism to increase the orbital excentricity of the planetes- 
imals is needed, in order to explain the high number of RACs detected per year; second, 
there is a much large number of RACs than BACs but, if the system was symmetrical, 
the same number of comets approaching and moving away from the star should be ob- 
served. Beust & Morbidelli (1996) possibly provided the best explanation assuming the 
existence of a Jupiter-like giant planet in a slightly elliptical orbit (M ~ IM/, r ~ 20 AU, 
e ^ 0.05) and the presence of a certain amount of asteroids in the disk. This planet would 
perturb the orbits of the planetesimals, by means of the (4:1) mean motion resonance, 
and would put them in high excentricity orbits, so they could become the comets re- 
quired by the FEB model. The excentricity of the planetary orbit would reproduce, for 
selected lines of sight, the observed statistical predominance of RACs over BACs. 

Summarizing, the observation of RACs in jS Pic is nowadays the best indication of 
the possible existence of planetesimals in MS stars (see the review paper by Lagrange 
et al., 2000). Furthermore, yS Pic was the strongest candidate star to harbour an extrasolar 
planet until the discovery of 51 Peg B by Mayor & Queloz (1995). 

1.5 Search for planetesimals: UXORs 

The discovery of the yS Pic disk motivated many astronomers to look for other CS disks 
in PMS and MS stars. The disk of yS Pic has some characteristics that make it quite easy 
to detect: it is nearby, large and seen edge-on. The following detections of CS disks 
were conducted about a decade later by O'dell et al. (1993), who identified CS disks 
using WFPC images taken with the Hubble Space Telescope (HST). They found about 
ten PROtoPLanetarY DiskS (proplyds) in the Orion nebula. The disks appeared as dark 
regions over the bright background of the H ii region created by the strong radiation field 
of the Trapezium cluster. This technique has been extended, using adaptive optics, to 
terrestrial telescopes. Another method used to detect CS disks is optical and infrared 
coronagraphy. It has been successfully performed with the HST and adaptive optics 
equipped terrestrial telescopes. The coronagraph is no longer needed in the mid-infrared 
and longer wavelength regimes. Some disks have been detected in direct images obtained 
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Figure 1.3: Simulation of a RAC via a FEB (taken from Beust et al., 1998). Top: a plan- 
etesimal approaches yS Pic (r ~ 0. 1 AU) in an elliptical orbit. Its outer layers evaporate 
and form a coma made of gas and dust which covers a significant fraction of the stellar 
disk. The observer is located at the right and the line of sight is represented by two 
horizontal dashed lines. Bottom: Two graphics show the computed absorption profiles 
for the K (left) and H (right) lines of Ca ii. [5 Pic is a high speed rotator, so its photo- 
spheric lines are much broader than the RACs and correspond to the pseudocontinuum 
upon which the absorptions are formed. Each line shows two CS components. The first 
one is a very narrow, intense and saturated absorption which has the same heliocentric 
radial velocity of the star (v - 20 km s"^). This is the stable absorption and cannot be re- 
produced by the FEB simulations, so it is introduced by hand in the spectra. The second 
absorption is originated by the comet coma. It is broader and shallower, but it is also 
saturated, as can be seen from the relative intensities of the two lines in the Ca n doublet. 
This component appears at a radial velocity of v ^ 45 km s~\ which is the radial velocity 
of the star plus the velocity of the planetesimal projected on the line of sight. 
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in the submillimetre and mid-infrared ranges. The highest spatial resolution has been 
reached with interferometric techniques, in the infrared and millimetre domains. 

The improvements made in the available instrumentation has allowed the observa- 
tion of many CS disks. Examples of primordial disks are those of TW Hya (Weinberger 
et al., 1999b; Krist et al., 2000), p Oph (Brandner et al., 2000), LkHor 101 (Tuthill et al., 
2001), MBM 12 (Jayawardhana et al, 2002) and Carina (Smith et al., 2003). Several 
second generation disks have been found in main sequence stars, e.g. HR 4796 (Koemer 
et al., 1998; Jayawardhana et al., 1998), Vega and Fomalhaut (Holland et al., 1998) and 
e Eri (Greaves et al., 1998). Some HAeBe stars near the Zero Age Main Sequence 
(ZAMS) also have second generation disks, e.g. AB Aur (Marsh et al., 1995; Naka- 
jima & Golimowski, 1995), HD 163296 (Mannings & Sargent, 1997) and HD 141569 
(Weinberger et al., 1999a; Augereau et al., 1999; Merm et al., 2004). 

It has been shown in the previous section that Lagrange-Henri et al. (1988) pro- 
posed that both the disks and the RACs in yS Pic are originated by the evaporation of 
planetesimals. This motivated some astronomers to look for CS gas in other stars as a 
complementary alternative to the direct detection of the disks (these searches were spe- 
cially intense when no other detection of CS disks was achieved). If the FEB hypothesis 
is applicable to other stars, the detection of TACs could suggest the existence of a CS 
disk edge-on or nearly edge-on. Hobbs (1986) and Lagrange-Henri etal. (1990) first 
conducted systematic searches of CS gas (TACs) in Vega, A-shell and F-shell MS stars. 
They used high resolution optical spectra and identified a bunch of candidate stars (e.g. 
HR 10 y 51 Oph). Some authors (e.g. Grady et al., 1991; Lecavelier Des Etangs et al., 
1997) successfully extended the search for TACs to the ultraviolet (UV) range. 

Later on, Grinin et al. (1994) detected TACs in visible spectra of the HAe star UX Ori, 
whereas Grady et al. (1995) found spectral variability in some UV lines. These results 
stimulated many authors (e.g. Grinin et al., 1996; Grady et al., 1996; de Winter et al., 

1999) to extend the search to other HAeBe stars (see the review paper by Grady et al., 

2000) . It was apparent that the HAe stars classified as UXORs were the most prone to 
present TACs in their spectra. As it will be explained below, the UXORs can be consid- 
ered the progenitors of 13 Pic stars, due to the similarity in mass, spectral type and disk 
orientation with respect to the line of sight. The study of the possible relation between 
TACs, CS disks and planetesimals in MS and PMS stars became a fundamental problem 
in the stellar formation and evolution field. 

The UXORs are PMS stars, mainly HAe, which show photometric variations similar 
to the star Algol. Hoffmeister (1949) first created a specific class to group all these 
objects together. The initial name, RW Aur stars, eventually became UXOR, because it 
was soon realized that the star UX Ori presented the Algol behaviour in a much stronger 
way. The UXORs alternate periods of low activity and almost constant brightness with 
severe dimming episodes (AV ~ 2.5). During the initial stages of darkening the stars 
redden, until there is a (turnaround) point after which the star, instead of getting redder, 
begins to turn bluer. Grinin et al. (1988) found out that the decreases in brightness were 
followed by increases in polarization up to values of Py ~ 6% in the deepest photometric 
minima. Figure 1.4, taken from Grinin et al. (1991), shows sample light curves of the 
UXOR star WW Vul, including a photometric minimum. 

Grinin et al. (1988) proposed the first model that explained successfully the photopo- 
larimetric behaviour of UXORs. According to Grinin et al. (1988), the light received by 
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Figure 1.4: Light curves for the UXOR star WW Vul, taken from Grinin et al. (1991). 
Top: Photopolarimetric data of WW Vul obtained by Grinin et al. (1991) during 300 
nights are shown. During this period, the star experienced an UXOR deep minimum 
Ay ^ 2. It can be appreciated that the drops in brightness are always matched by in- 
creases in polarization, up to values of Py - 6% in the photometric minimum. Bottom: 
Several colour-magnitude diagrams are shown for the star. In a darkening event, the star 
departs from the top of the diagrams (maximum brightness), goes down to lower bright- 
ness and, eventually reaches the photometric minimum. It then goes up until it recovers 
its original luminosity. At first the brightness follows a linear extinction law, until the 
extinction reaches AV ^ 1. For greater dimming, the reddening becomes increasingly 
lower than that given by the linear law, i.e. the slope of the curve grows or, even, changes 
its sign. The biggest variations are found in the colour indices U-B and B-V, for which 
the star becomes even bluer when the brightness decreases. 
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an observed is the sum of two components, namely the stellar radiation and the light scat- 
tered by an (almost) edge-on CS disk. The photopolarimetric variations are attributed 
to the passage of an opaque dust cloud in front of the line of sight. The cloud dims the 
stellar light according to an interstellar extinction law, i.e. the star becomes redder when 
fainter. On the other hand, the disk luminosity is assumed to be constant. As the disk 
is seen in scattered light, it is supposed to be bluer than the star and to have a non-zero 
polarization. 

According to Grinin et al. (1988), the maximum brightness corresponds to instants 
where no cloud is crossing the line of sight, so the contribution of the disk to the total 
luminosity is negligible and only the stellar light is observable. When a big cloud crosses 
the line of sight, the stellar brightness fades following the extinction law, so the star 
draws a straight line in the upper part of the colour-magnitude diagrams. If the cloud 
is large enough, the amount of light scattered by the disk can be comparable to, or 
greater than, that emitted by the star. This happens during deep minima, AV > 1, and 
implies both a blueing of the observed light and an increase in the polarization. On the 
other hand, the luminosity of the disk places a limit on the depth of the photometric 
minima (AV^ax - 2.5), because the disk, due to its large dimensions, remains visible 
even when the star has been fully obscured (i.e., the role of the cloud is similar to that of 
a coronagraph). 

Grinin et al. (1991) proposed that dust clouds are common in all HAe dense proto- 
planetary disks. Therefore, the UXOR phenomenon would be observed in all PMS stars 

with edge-on disks, i.e. the UXOR phenomenon is purely geometrical. This hypothesis 
receives additional support from the fact that, apart from the deep photometric minima, 
the UXORS are very similar to other non-UXOR HAe stars (Natta et al., 2000a). 

Natta & Whitney (2000) made some numerical simulations of the UXOR behaviour. 
They used disk models more realistic than those utilized by Grinin et al. (1988) (e.g. they 
considered flared disks). Natta & Whitney (2000) found out that the UXOR phenomenon 
should be observable for those stars whose disks have an inclination respect to the line 
of sight of i ~ 45-68°. Lower inclinations do not allow a full occultation of the stellar 
light and larger angles cause too much extinction for the star to be optically detected. 
A simple statistical calculation shows that about a half of the stars should present the 
UXOR phenomenon. Currently, the best non-biased photometric sample of HAeBe stars 
is that of the Hipparcos catalogue, van den Ancker et al. (1998) carefully studied the 
variability presented by each of the HAeBe stars observed with the satellite. It can be 
drawn from their study that about 30% of the analysed HAe stars show a photometric 
variability greater than 0.5 mag. This figure is roughly compatible with the prediction of 
50% given by Natta & Whitney (2000). 

The TACs in UXORs were initially interpreted as the undisputable trace of the pres- 
ence of planetesimals (Grady et al., 2000, see references therein). However, this inter- 
pretation has been recently questioned because the UXORs, contrary to yS Pic, present 
spectroscopic activity in the Balmer and Nai lines. Nai atoms cannot survive more than 
a few seconds in the strong radiation field of an A type star. This is hardly compatible 
with a cloud of a size comparable to the stellar disk, because the atoms should have been 
travelling about a day since their ejection from the comet. Sorelli et al. (1996) found that 
this would only be possible if the planetesimals were significantly larger than in yS Pic or 
if they fully evaporate in a single approach to the star. 
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On the other hand, the detection of TACs in the Balmer series is not consistent with 
a cometary composition for the CS gas. Natta et al. (2000b) used a Non-Local Ther- 
modynamical Equilibrium (NLTE) code to show that the abundances of the gas cloud 
which originated some TACs in UX Ori are solar or nearly solar. This is not compati- 
ble with the high metallicities (m > 500) expected in highly hydrogen-depleted bodies 
such as planetesimals. Besides, Beust et al. (2001) showed that the FEB mechanism is 
not capable of generating TACs in HAeBe stars, because the strong stellar winds highly 
coUimate the comets comae, so they could not cover a significant fraction of the stellar 
disk. The evaporation of comets would only be efficient if they approach to the star in 
wind free cavities. Finally, Hartmann et al. (1994) suggest that the RACs in CTTSs are 
formed in a completely different scenario: magnetospheric accretion of gas. The latter 
possibility is explored in the next section. 

Summarizing, it can be said that prior to this thesis the relation between TACs and 
planetesimals was not clear except for p Pic. This interesting question has been one of 
the main guidelines followed in the research conducted in the present work. 

1.6 Magnetospheric accretion 

Classical T Tauri stars present an intense photometric and spectroscopic activity. Several 
of these phenomena can be explained in terms of the magnetospheric accretion of gas 
from the CS disk: brightness increments modulated by the rotation period (hot spots), 
excess of UV radiation (veiling), significant emission in hydrogen and metallic lines and 
presence of RACs (the simultaneous observation of extended line emission and a RAC 
is called an inverse P Cygni profile). In this Section, the fundamentals of the magneto- 
spheric accretion theory and its relation with TACs in HAe stars will be reviewed. 

Hartmann et al. (1994) proposed the first self-consistent theoretical model capable 
of computing emission line profiles generated by magnetospheric accretion in CTTSs. 
These calculations have been later improved (MuzeroUe et al., 1998a, 2001) and the 
resulting profiles have been compared to real CTTSs spectra (MuzeroUe et al., 1998b, 
2001). The main hypothesis in the model is that the material accumulated in the circum- 
stellar disk falls onto the star along the stellar magnetic field lines. The kinetic energy 
acquired by the infalling matter is released, after the collision with the stellar surface, 
as UV radiation. This model qualitatively explains the presence of hot spots (impact 
regions in the stellar surface), the veiling of the spectra and the strong emission present 
in some lines. The material gains a significant velocity during the infall, so it could be 
expected the presence of RACs, due to the absorption of the continuum by the gas, for 
lines of sight corresponding to large inclinations of the disk. Other phenomena, like the 
possible appearance of high velocity BACs, are out of the scope of the accretion models 
and can only be explained if high velocity winds are assumed. 

Hartmann et al. (1994) postulated the existence of a perfectly dipolar magnetic field 
whose force lines are bound to the disk and rotate in a rigid body-like picture. It is 
assumed that the magnetic field destroys the structure of the disk and channels the CS 
matter to the stellar surface. In this way, the gas undergoes a free fall following a tra- 
jectory given by the field lines. The magnetosphere has axial symmetry and extends 
from an initial radius to an outer radius. This symmetry forces the profiles to be non 
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Figure 1.5: Geometry assumed in the magnetospheric accretion models (taken from 
Hartmann et al., 1994). The figure represents a section in the radial direction, which 
includes the star (half circle in the left), the disk (bar in the right) and the magnetosphere 
(field lines). The stellar magnetic field is assumed to be purely dipolar and its intersection 
with the disk is bound to the range R„ = 2.2-3.3 R^,. The CS gas located inside the 
magnetosphere experiences a free fall to the star along the field lines. The kinetic energy 
acquired during the fall is released as UV and visible radiation when the material hits 
the stellar surface. This is the origin of the line veiling. A detailed study of the radiative 
transfer reproduces the emission profiles and RACs observed in many lines in CTTSs. 



time- variable, because the accretion is evenly distributed along the whole disk. The disk 
is assumed to be thin and completely opaque. Figure 1.5, taken from Hartmann et al. 
(1994), shows the geometry assumed by the theory. 

The gas is considerably heated during its free-fall (T ~ 8000 K). Currently, this pro- 
cess is not fully understood. This motivated Hartmann et al. (1994) to use semiempirical 
heating laws that relate the gas temperature with the radial distance and the departure 
point from the disk. These temperature distributions are relatively constant and are char- 
acterized by the highest temperature achieved by the gas (T ~ 6000-10000 K). The 
density along each trajectory can be calculated if the accretion rate is known. Once the 
density, temperature and velocity of the gas are calculated for every point in the magne- 
tosphere, a detailed radiative transfer calculation can be performed in order to determine 
the contribution of each parcel of accreted material to the circumstellar profile of several 
lines. This computation is complex and uses the extended Sobolev method (MuzeroUe 
et al., 2001). In order to make the calculation easier, the star and the shock region in 
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the stellar surface are replaced by black bodies of appropriate temperature and area. The 
free parameters in the model are the stellar mass, the stellar radius, the rotation velocity, 
the inner and outer magnetosphere radii, the accretion rate, the maximum temperature 
achieved by the gas and the inclination of the rotation axis with respect to the line of 
sight. 

MuzeroUe et al. (2001) have performed a thorough study of the parameter space, 
in order to find what regions are compatible with the observations. They have used a 
typical CTTS with = 0.5Mq, i?^ = 2Rq. The remaining parameters are in the fol- 
lowing ranges: Accretion rate M : 10~^-10~^Mo/yr, maximum temperature of the gas 
Tmax : 6000-10000 K, inclination of the line of sight / : 10-75°, rotation velocity : 0- 
25 kms"^ Four different prototypical magnetospheres were considered: small/wide 
(Rm = 2.2-3 R^), small/narrow (R„ = 2.8-3 R^), large/wide (R„ = 5.2-6 R^) and 
large/narrow = 5.8-6 7?^). 

In order to confirm the validity of the calculations, MuzeroUe et al. (2001) tried to 
reproduce the observed profiles and the veiling factors for Hor and the Na i D doublet for 
several CTTSs. They found a reasonable agreement for the stars with intermediate-low 
accretion rates. However, the Ha models fail to reproduce the profiles for the higher 
accretion rate objects (e.g. DR Tau, M = 5 x lO^^Mg/yr), which have so strong winds 
that the lines display classical P Cygni profiles that hide the effects of accretion. 

The best results are obtained for the star BP Tau and are shown in Figure 1.6. It 
can be appreciated that both the shape of the profiles, including the large wings due to 
the Stark broadening, and the approximate emission fluxes have been reproduced. The 
Nai line models reproduce the observed inverse P Cygni profiles. The parameters used 
in the calculation are M = 10 *^Mo/yr, R^ = 2.8-3 R^, T^^ = 8000 K, i = 70° and 

= 10 kms~^ The Hfi and Hy spectra were obtained 4 years after those of Ha and 
Na I D. However, the agreement is good, which suggests that the accretion rate of the 
object has not substantially varied during that interval. 

It has been shown that the magnetospheric accretion models reproduce properly the 
main features of the CTTSs spectra. In particular, the RACs of the inverse P Cygni pro- 
files observed in many CTTSs are obtained both in hydrogen and metallic lines. This 
achievement motivated Sorelli et al. (1996) and Natta et al. (2000b) to claim that the 
presence of TACs in UXORs could be explained in terms of non axisymmetric magne- 
tospheric accretion. Therefore, the infall of material onto the star would be restricted to 
a kind of magnetic field tubes (funnels). This hypothesis explains the spectral variability 
in terms of the creation, destruction and rotation of accretion funnels. 

Despite its success in CTTSs, the profiles calculated by MuzeroUe et al. (2001) are 
not directly applicable to intermediate mass UXOR stars, i.e. the kind of objects studied 
in this thesis. First, it is not clear that these stars can sustain a stable magnetosphere, 
with fields strong enough and rigidly anchored to the CS disk. Second, the high rotation 
velocity of these objects, namely ~ 200 kms"\ would imply a very quick revolu- 
tion of the magnetosphere, which would substantially modify the synthetic line profiles. 
Third, the line profiles studied in this thesis do not present much emission, which is re- 
stricted to a few lines in the Balmer series (Hfi, Hy and H5). That is, the profiles are 
composed of different overlapping CS absorptions without emission, so they are very 
different from the inverse P Cygni profiles observed in typical CTTSs. Finally, the TACs 
exhibit a fast evolution, in times of days and even hours without any apparent periodicity. 
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Figure 1.6: Comparison between synthetic (dashed lines) and observed (solid lines) ac- 
cretion profiles in the CTTS BP Tau (taken from MuzeroUe et al., 2001). It can be appre- 
ciated a good agreement in the Balmer lines, both in the profile shape and the total emit- 
ted flux. Inverse P Cygni profiles, composed of extended emission and a RAC, are repro- 
duced in the Nai D lines. The parameters used in the calculation are M = IO'^Mq/jt, 
Rm = 2.8-3 R^, Tniax = 8000 K, / = 70° and = 10 km s'K The spectra of Bfi and Hy 
were obtained 4 years after those of Ha and Na i D. However, the agreement is good, 
which suggests that the accretion rate of the object has not substantially varied during 
that interval. 
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Figure 1.7: Emission and absorption H/? CS profiles generated by magnetospheric ac- 
cretion, (taken from MuzeroUe et al., 2001). These profiles have been calculated for a 
typical CTTS with = O.SM©, i?^ = 2Rq, and a magnetosphere of type small/wide 
Rm = 2.2-3 Ri,. The rotation axis is inclined 60°, a typical value for UXOR stars. It 
can be seen that it is necessary to combine a high (low) accretion rate with a low (high) 
maximum temperature in order to achieve the strong emissions observed in T Tauri stars. 
The profiles with greater resemblance to the ones observed in UXOR stars are those ob- 
tained with the lowest temperatures and accretion rates, which are characterized by an 
almost absence of emission combined with narrow absorption components. 



This temporal variability cannot be explained in terms of axisymmetric steady models. 

In Figure 1.7, taken from MuzeroUe et al. (2001), the models with the closest re- 
lationship to the UXOR stars are shown. The displayed profiles correspond to the 
line and have been computed for a typical CTTS with = O.SMq, R^, = 2Rq and 
Vi, = kms"^ The magnetosphere is of type small/wide (Rm = 2.2-37?*). The line of 
sight is inclined 60° respect to the rotation axis. This is a typical value for HAe stars 
with UXOR activity (Natta & Whitney, 2000). The free parameters in the simulation are 
the accretion rate, M, and the maximum temperature, T^ax- The profiles shown only rep- 
resent the CS contribution to the lines, excluding the stellar photosphere and the veiling 
continuum. 

The parameter space can be divided into three loose regions with different types of 

line profiles: high accretion rate and temperature, low accretion rate and temperature 
and, finally, high (low) accretion rate and low (high) temperature. The inverse P Cygni 
profiles are obtained when one parameter is high and the other one is kept low. This 
implies some restrictions to the models that have been extensively studied by Muze- 
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roUe et al. (2001). If a high accretion rate and a high temperature are combined, the 
material flux produces a greater level of continuum veiling than that generated in the 
shock region. This produces strong absorptions, including an extremely broad compo- 
nent (Av ~ 500 kms~^). Finally, for low accretion rates and temperatures, low emission 
profiles with redshifted narrow absorption components (Av ~ 100 kms~^) are obtained. 
Only the later case is similar to the observed profiles in UXOR stars, which will be 
thoroughly studied in Chapters 5 and 6. If the models explained in this Section were 
applicable to HAe stars, it would be deduced from the observations that the physics of 
accretion is less violent in UXORs than in CTTSs. 

1.7 Evolution of the subject of the present PhD thesis 

The structure and contents of the present PhD thesis considerably differ from the initial 
approach followed by the author and his advisor. This a natural consequence of the 
research process, because only after the analysis of the data began, it was realized that 
the initial approach had to be modified and the real power of the observational data could 
be established. In this Section the evolution that the thesis has followed until it reached 
its actual structure and contents will be exposed, as well as the reasons that motivated 
those changes. 

The initial objective was the study of the evolution of the properties of planetesimals 

in CS disks, from the beginning of the PMS stage until the MS. It was expected to find 
new observational physical restrictions that could be included in the theoretical models 
of star formation. It was initially assumed, as a working hypothesis, that the TACs are 
the footprints of the presence of cometesimals, both in PMS and MS stars. In this way, 
the thesis was defined as the search and study of TACs in a significative sample of stars 
(broad range of ages and masses) with CS disks. 

The adopted observational strategy was to obtain a large amount of high resolution 
echelle spectra for a big sample of stars. 12 observing nights were used for this program 
in the framework of the EXPORT collaboration (see details in Chapter 2). The WHT 
telescope (4.2m, La Palma) equipped with the UES spectrograph was used. After the 
spectra were reduced, a database composed by 198 echelle spectra of 49 stars with pos- 
sible CS disks was available, both primordial (PMS stars) and second generation (MS 
stars). The PMS sample of stars includes UXORs, HAeBe with observed disks, HAeBe 
with spectroscopic behaviour similar to {3 Pic, CTTSs, WTTSs and Early T Tauri Stars 
(ETTSs). The observed MS stars are of types Vega, A-shell with possible Pic activity 
and Post T Tauri Stars (PTTSs) taken from the catalog of Lindroos (1986). 

The initial EXPORT sample of objetcs included in the observing proposal was biased 
towards intermediate mass objects: only a 30% of the stars were T Tauri or PTTSs. On 
the other hand, the sensitivity of the instrument together with the low brightness of the 
late spectral type stars made the observation of low mass stars very difficult: only about 
20% of the total number of spectra obtained correspond to T Tauri stars or PTTSs, and 
many of them have a low Signal to Noise Ratio (SNR). These restrictions, imposed by 
the sample studied and the available instrumentation, have limited the research in this 
thesis to the intermediate mass stars, mainly of spectral type A. 

A preliminary analysis, made after the reduction of the data, revealed that most of the 
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MS stars do not present any TAC, with the noteworthy exception of HR 10, an A-shell 
star similar to yS Pic (see Section 1.4). However, there were only two spectra of HR 10 
available, taken within an interval of three months. These two spectra were not enough 
to make a thorough study of the properties and interpretation of TACs in MS stars, so 
this kind of objects was excluded from this thesis because, on the one hand, the sample 
would not be significant and, on the other hand, the time interval between the two spectra 
is excessively large. It will be seen in Chapters 5 and 6 that the TACs vary in timescales 
of days and even hours. 

Summarizing, after the preliminary analysis was conducted, three unavoidable re- 
quirements were identified in order to study the presence of planetesimals in CS disks: 
presence of TACs in the spectra, moderate to high SNR and good temporal coverage 
(many spectra per star with an approximate time delay of one day, during many nights). 
This restricted the list of possible objects of study basically to a few HAe stars, most 
of them of UXOR type. Finally, the five objects with the best temporal monitoring and 
larger spectral variability of the sample were selected: BF Orionis, SV Cephei, UX On- 
onis, WW Vulpeculae and XY Persei. The star HD 163296 was also studied for a while. 
However, its strong stellar winds produce considerable alterations of the line profiles, 
so it was not possible to make an adequate characterization of the TACs present in its 
spectra. 

BF Ori, UX Ori, WW Vul (e.g. see Grinin et al., 1991) and SV Cep (Rostopchina 
et al., 2000) have been classified as UXORs. Oudmaijer et al. (2001) found a strong 
photopolarimetric variability in XY Per, but they could not undoubtedly establish its 
membership to the UXOR class. Since the spectral variability of XY Per is, according 
to the TACs, identical to that of the remaining objects studied, from now on it will be 
considered an UXOR for all practical purposes. 

Once the initial expectations of this study were delimited to the detection and char- 
acterization of planetesimals in the EXPORT subsample of UXOR stars, the question 
of the relation between TACs and planetesimals already remained. As it has been said 
before, Natta et al. (2000b) demonstrated that, for a series of strong TACs detected in 
UX Ori, the chemical abundances were approximately solar and not compatible with the 
evaporation of planetesimals. The simulations made by Beust et al. (2001) also led to 
the same conclusion, so the link between TACs and planetesimals became doubtful. Due 
to the high relevance these results have for the present thesis, it was initiated a collabora- 
tion with Dr. Natta (Astrophysical Observatory of Arcetri, Italy) to thoroughly study the 
connection between TACs and planetesimals. We put in common the extensive experi- 
ence with theoretical models of Dr. Natta and our collection of excellent quality echelle 
spectra and our capability of synthesize photospheric spectra. 

A substantial part of the collaborative efforts, devoted to an exhaustive study of the 
TACs in UX Ori, were made during two short stays in the Astrophysical Observatory of 
Arcetri in 2000 (2 months) and 2001 (3 months) It was verified that the TACs are not 
generally related to planetesimals but to CS gas clumps of solar metallicity, probably 
related to magnetospheric accretion phenomena. It was also discovered that the extraor- 
dinary temporal coverage of the EXPORT spectra allowed to study, with an accuracy 
without precedents, the kinematics of the gas originating the TACs in UX Ori. These 
results have been published by Mora et al. (2002, from now on Paper I). 

Once it was realized the success of the analysis of UX Ori, a similar study for the 
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remaining selected stars was started: BF Ori, SV Cep, WW Vul and XY Per. This new 
work has been successfully finished and has been published by Mora et al. (2004, from 
now on Paper II). The conclusions obtained in Paper I and Paper II do not substantially 
support the link between TACs and planetesimals, so it was needed to change again 
the subject and title of the thesis to its final form: "Kinematics of the circumstellar gas 
around UXOR stars". 

1.8 Objectives and carried out work 

The research carried out in this PhD thesis has not been developed in a linear way (see 
previous Section). That is, there existed some initial objectives, according to the relevant 
part in the EXPORT observational proposal, but they were modified many times as a 
consequence of the partial results obtained during the research. Therefore, in this Section 
we do not describe the (obsolete) primitive objectives, but the significant milestones 
achieved during the whole PhD research process. 

• The first duty accomplished in this thesis was the obtention of the required obser- 
vational data: a set of echelle spectra. The observations were carried out in the 
framework of the EXPORT collaboration. The student joined the collaboration 
when the observational proposal was written, the observing time allocated and 
many of the observations (2 out of 12 nights) performed (see Chapter 2). The raw 
data were available from the very beginning of the research. 

• Thanks to the good weather during the spectroscopic observations (100% of clear 
nights), a total amount of 198 echelle spectra of many objects was obtained. The 
reduction of the spectra required, along with the observing runs, the first year and 
a half of this thesis, and is described in Chapter 3. 

• Once the data were reduced, a preliminary inspection of the spectra, in order to 
classify the observed objects and estimate their interest was started. The first re- 
sults of this study were presented by Mora et al. (2000). The selection process 
of interesting objects has been long, because of the diff'erent changes in the thesis 
subject, due to the partial results that were obtained. 

• In order to characterize properly the TACs in the reduced spectra, it was necessary 
to determine and subtract the photospheric spectrum of the stars studied. It was 
decided to use synthetic spectra to estimate the photospheric spectra. The codes 
SYNTHE and SYNSPEC were used for this purpose, so they had to be properly 
installed and, sometimes, modified. The spectral synthesis process is described in 
detail in Sections 5.3.1 and 6.3.1. 

• A fundamental parameter in the computation of stellar spectra is, along with the 
effective temperature and the gravity, the rotation velocity projected on the line 
of sight, vsinz. A working group by Dr. Montesinos, Dr. Solano (who was the 
coordinator) and the student was constituted with the only purpose of systemati- 
cally determine v sin / for all the stars in the EXPORT sample. In this way, besides 
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the obtention of fundamental parameters for the uherior analysis of the data, a re- 
search of a great scientific value on its own was carried out. The results have been 
published by Mora et al. (2001), along with the spectral type determinations for 
all the stars in the sample made by Dr. Merm. In Chapter 4 a detailed exposition 
of the method employed and the results obtained for the stars in this thesis can be 
found. 

• After the best photospheric spectrum is obtained for each star, it is subtracted 
from the observed spectrum. The study of the subtraction residual R has been 
done in several metallic and hydrogen lines and revealed the presence on many 
TACs in the spectra. In general, there appeared several overlapped TACs in each 
spectrum. In order to characterize simultaneously and self-consistently all the 
TACs multicomponent gaussian fits were used. Estimations of the radial velocity 
V, velocity dispersion Av and intensity R^^ were obtained for each TAC. This 
process is described in Section 5.3.2 

• The study of the the multigaussian deconvolution data revealed that absorption 
components detected simultaneously in different lines can be grouped according 
to their radial velocity. The origin and variability of these absorptions have been 
attributed to the dynamical evolution (acceleration or deceleration) of gas clouds 
present in the CS disks. The kinematics of CS gas clumps in UXOR stars has been 
traced with a detail without precedents. The method used in the identification of 
these clouds is exposed in Sections 5.4 and 6.3.2. 

• Apart from the infall and outflow velocities of the CS gas, the dispersion velocity 
and the intensity of the absorptions have also been studied in detail. The possible 
cross-relationships between different parameters have also been analysed. Highly 
significant correlations between v vs. Av and Av vs. i^max have been found (see 
Sections 5.4.1 and 6.4.1). 

• The study of the Rjaax values allowed us to examine what line multiplets are sat- 
urated for each gas blob and each star. These results impose restrictions to the 
physical parameters of the cloud, remarkably the density of the CS medium, which 
should be included in any realistic model explaining the existence of TACs in 
UXOR stars (e.g. FEBs or magnetospheric accretion). In some stars (UX Ori and 
XY Per) the intensity ratios in the Balmer lines have been related to the possible 
emission of radiation by the CS cloud. The most relevant results can be found in 
Sections 5.4.2, 5.5.1 and 6.4.2. 

• Finally, a metallicity analysis similar to that of Natta et al. (2000b) has been per- 
formed for the whole sample of stars in this thesis. The results show that, in 
general, the CS gas clouds in UXOR stars have a solar or nearly solar metallic- 
ity. Otherwise, WW Vul seems to be a different object because, in addition to 
solar metallicity TACs, it shows metallic absorption components without an evi- 
dent counterpart in hydrogen lines. The significance of these results is discussed 
in Sections 5.5.1, 5.5.2, 6.4.3 and 6.4.4. 
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1.9 Brief description of the structure of this thesis 

Chapter 2 explains all the information related to the observations used in this thesis: 
The initial observing proposal, the observations themselves and the final set of spectra. 
Chapter 3 shows the details of the reduction procedures used in this thesis. Chapter 4 
explains the method used by Mora et al. (2001) in the measurement of the projected 
rotation velocities for all the EXPORT stars. Chapter 5 includes the complete study 
of UX Ori, performed by Mora et al. (2002) and originally published in Astronomy & 
Astrophysics. Chapter 6 consists of the the work carried out by Mora et al. (2004) about 
BF Ori, SV Cep, WW Vul and XY Per, which was originally published in Astronomy & 
Astrophysics. Finally, Chapter 7 contains the conclusions obtained in this PhD thesis. 



Chapter 2 
Observations 



In this Chapter, the observational campaigns carried out to obtain the spectra used in 
this thesis will be described. Section 2.1 explains the scientific case of the observing 
proposal, composed of three different subprograms. The data used in the present research 
are only a fraction of those obtained in the subproposal devoted to the study of CS disks, 
whose main features are discussed in Section 2.2. Section 2.3 contains some general 
remarks about the observations and the database generated with the reduced spectra. The 
detailed information relative to the observations of the stars studied in this thesis (number 
and distribution of the spectra, SNR, etc.) are not discussed here but in Chapters 5 and 6 
(Sections 5.2 and 6.2). 



2.1 The EXPORT observational proposal 

The spectra used in this thesis were obtained during the observing runs awarded to the 
EXPORT collaboration^ during the International Time (see below) in 1998 and 1999. 
One of the main objectives of the EXPORT proposal was to perform a precursor scien- 
tific study for the DARWIN mission of the European Space Agency (ESA). DARWIN 
is one of the most ambitious projects currently managed by ESA. The mission involves 
the development of an infrared interferometer in space capable of detecting and charac- 
terizing extrasolar telluric planets. The DARWIN project, apart from being extremely 
expensive, requires the development of several technological precursor missions, which 
would address the viability of the new technologies to be implemented in the real mis- 
sion. Therefore, there is not currently a realistic schedule to carry out the project. The 
EXPORT proposal (Eiroa et al., 2000b), titled "Planetary systems: their formation and 
properties", tried to contribute to the scientific case for DARWIN, conducting astronom- 
ical research of high scientific value on their own. 

The observing proposal was submitted in 1997. At that time, the detections of extra- 
solar planets that came after the discovery of 51 Peg B by Mayor & Queloz (1995) were 
questioning the validity of the accepted star formation paradigm. The EXPORT collabo- 
ration, composed of about twenty European and American astronomers, was constituted 

'The web page of the EXPORT collaboration, http://laeff.esa.es/EXPORT/, contains detailed 
information about its objectives, proposal, observations and results. 
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Table 2. 1 : EXPORT observations schedule, excluded the microlensing alert time avail- 
able in the IAC-80 telescope. WHT (4.2 m): high resolution spectroscopy. INT (2.5 m): 
intermediate resolution spectroscopy. NOT (2.5 m): UBVRI photopolarimetry. CST 
(1.5 m): JHK photometry. JKT (1.0 m): CCD images. 



Telescope 


May 1998 


July 1998 


October 1998 


January 1999 


WHT 


u-n 


28-31 


23-26 


28-31 


INT 


14-17 


28-31 


24-28 


29-31 


NOT 


14-17 


28-31 


23-27 


29-31 


CST 


15-17 


28-31 


23-26 


28-31 


JKT 




25-31 


24-1 





with the main objective of applying for the "International Time"^ available at the Canary 
Islands telescopes in 1998. The proposal, which was awarded with the total International 
Time, consisted of the following research lines: 

• Formation and evolution of planetary systems via observational studies of pro- 
toplanetary disks around PMS and MS Vega-type stars. The allocated telescope 
time was: 12 nights with the WHT (4.2 m, UES high resolution echelle spectro- 
graph), 16 nights with the INT (2.5 m, IDS intermediate resolution spectrograph), 
16 nights with the NOT (2.5 m, TURPOL photopolarimeter) and 15 nights with 
the CST (1 .5 m, infrared photometer and CAIN infrared camera). 

• Search for planetary atmosphere signatures in the spectra of r Boo and 51 Peg. 
This subprogramme was assigned 4 nights with the WHT telescope equipped with 
the UES spectrograph. 

• Search for new exoplanets by means of photometric transits in stellar clusters and 
microlensing techniques: 16 nights with the JKT telescope (1 m, CCD camera) 
were allocated for the quest for planets via transits. A certain amount of alert time 
for the microlensing searches for planets was allocated in the IAC-80 telescope 
(0.8 m, CCD camera). 

Table 2. 1 shows a summary of all the observing time granted to the EXPORT collab- 
oration (excluding the microlensing alert time with the IAC-80 telescope). The weather 
was quite good: 100% of the nights were adequate for the spectroscopy programmes and 
more than half of the nights were photometric. 

2.2 Formation and evolution of circumstellar disks sub- 
programme 

The objectives of this subprogramme were the following: 

^The "International Time" amounts the 5% of the total observing time available 
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1. Study of the evolution of CS disks, from the PMS (dense protoplanetary disks) to 
the MS stages (second generation dust disks). 

2. Characterization of the planetesimals during the whole disk lifetime: appearance- 
disappearance, absorption profiles and chemical abundances. As a working hy- 
pothesis it was assumed that the TACs are the signature of the evaporation of 
planetesimals (as it probably happens in Pic) both in MS and PMS stars. 

The selected sample consists of stars with possible CS disks: UXOR with detected 
TACs, PMS with possible UXOR photopolarimetric behaviour, HAe near the ZAMS, 
Vega-type, A-shell with detected TACs and PTTSs. In order to be as more representative 
as possible, the sample studied includes objects with a broad range of ages and stellar 
masses. However, the sample is biased towards intermediate mass stars. 

It has been shown in the previous Section that the observations allocated to this 
subprogramme consist of optical high and intermediate resolution spectroscopy, optical 
photopolarimetry and infrared photometry. Many of the objects studied have also been 
observed with the infrared space telescope ISO, whose data archive was released in 1999. 
Because of their great interest, the available mid-infrared ISO spectra, obtained with the 
SWS spectrograph, were studied. 

The study of the disks required the photopolarimetric data and the ISO spectra. The 
initial objective was to combine the construction of excellent quality SEDs (due to the 
simultaneity of the photometric data) with the study of the solid state features present 
in the ISO spectra. The SED construction for all the stars in the sample is one of the 
most interesting issues of the PhD thesis of Dr Merm, which has been recently defended 
(Merin, 2004). The study of the ISO spectra has been restricted to the characterization 
of the A ~ lOyum silicate emission band (Palacios et al., 2000). 

The characterization of the TACs for the EXPORT sample is exactly the subject of 
this thesis. In principle, all the available observations should have been used. The in- 
termediate resolution spectra were expected to be the most important tool, because they 
should have allowed to perform the majority of the detections and statistics of TACs. 
The high resolution spectra would have been dedicated to the detailed characterization 
(kinematics, profiles, abundances) of the most interesting TACs. The visible photopo- 
larimetry and the infrared photometry should have been used to study the relation be- 
tween the photometric variability (e.g. UXOR minima) and the spectroscopic activity. 
Finally, the study of the ISO spectra should have revealed traces of the CS solid material. 

It has been shown in Section 1.7 that, after an initial inspection of the data, it was 
decided to study the TACs only in UXOR stars. The variability showed by the stars 
studied was small, with the remarkable absence of deep photometric minima. No corre- 
lation between the photopolarimetric and spectroscopic variabilities was detected. It was 
also realized that, in order to characterize properly TACs (number, velocity, width and 
intensity), high resolution was absolutely needed. Therefore, the intermediate resolution 
spectra were excluded from the study. Finally, when it was confirmed that the TACs 
in UXOR stars are not generally related to the evaporation of planetesimals but to gas 
accretion, it was decided to defer the study of the ISO spectra for a future research. In 
summary, the only observations used in the detailed analysis of the TACs are the high 
resolution echelle spectra. 



26 



Chapter 2. Observations 



2.3 EXPORT observations of echelle spectra 

It has been said in Section 2.1 that 12 nights of telescope were allocated for the echelle 
observations devoted to the evolution of CS disks. The nights were scheduled in 4 blocks, 
separated by 3 months intervals: 2 nights in May 1998, 4 in July 1998, 2 in October 1998 
and 4 in January 1999 (Mora et al., 2000). This strategy permitted to perform long-term, 
short-term and very- short-term monitoring (months, days and hours, respectively) of the 
most interesting objects. 

The instrument used was the "Utrecht Echelle Spectrograph" (UES), which was per- 
manently located in a Nasmyth focus of the "William Herschel Telescope" (WHT, 4.2 m) 
during all its operational life. The selected configuration of UES provided spectra in the 
range 3800-5900 A with a spectral resolution R = A/AA = 49 000. A CCD detector 
of 2048 X 2048 pixels was used to record the spectra. The overlap between contiguous 
orders rendered a continuous coverage in the mentioned spectral range. The projection 
on the sky of the selected entrance slit was 1.15 arcseconds. 

The reduced spectra consist of a set of 198 spectra of 49 astronomical objects. Ta- 
ble 2.2 shows a detailed list of them. The classification in the table corresponds to that 
made by Mora et al. (2001), according to the spectral type, SED and photopolarimetric 
behaviour of each star. The integration times range from 10 minutes for the brightest 
objects (e.g. HR 10, V = 6.23, SNR ~ 2400) to 45 minutes for the faintest stars (e.g. 
VV Ser, V ~ 11.9, SNR ~ 14). The observed limiting magnitude is about V - 12 and 
corresponds to 45 minutes exposures. Larger integrations have not been performed, be- 
cause of the cosmic rays (~2% of affected pixels for 45 minutes exposures) and the need 
to obtain spectra for a large sample of stars. 

Figure 2. 1 shows a histogram of the number of UES spectra obtained per star. It can 
be appreciated that a significant number of objects was observed only once. This was 
due to two reasons. The first one is that the intermediate resolution INT spectra obtained 
simultaneously were real-time analysed to evaluate the relative interest of the objects. 
In this way, those stars not showing significant traces of CS material in the INT spec- 
tra were rejected for subsequent observations. This approach has been generally valid. 
However there are exceptions, e.g. RY Ori (Eiroa et al., 2000a) presents a normal inter- 
mediate resolution spectra, but shows a large amount of TACs in its only high resolution 
spectrum. The second reason is that some stars were too faint to obtain enough high 
SNR UES spectra. Using the real-time information provided by the INT and observing 
with the highest possible frequency those objects of greater apparent interest allowed us 
to optimize the observational strategy. 

Currently, all the data obtained by the EXPORT collaboration are available to the 
scientific community. All the photopolarimetric data have been published and analysed 
(Oudmaijer et al., 2001; Eiroa et al., 2001, 2002). Many studies have been done to ex- 
ploit the spectroscopic data obtained by the collaboration (Mora et al., 2001, 2002, 2004; 
Merm et al., 2004; Baines et al., 2004). Because of the large amount of data involved 
and the peculiarities of the reduction, the spectra are only available under request. 
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Table 2.2: List of high resolution echelle spectra obtained by the EXPORT collaboration 
with the WHT telescope (4.2 m, La Palma). The classification shown corresponds to that 
made by Mora et al. (2001), according to the spectral type, SED and photopolarimetric 
behaviour of each star. Column 3 includes the total number of exposures obtained per 
star. 
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Number of exposures 



Figure 2.1: Histogram of the number of spectra obtained per star. A significant number 
of stars was observed only once, because of their apparent lack of CS material (evalu- 
ated in real-time from the inspection of the intermediate resolution INT spectra) or their 
low brightness. Thanks to the information gathered from the intermediate resolution 
spectra, the observations could be optimized and the most interesting objects were more 
frequently observed. 



Chapter 3 

Reduction of the observational data 



In this chapter, the method used in the reduction of the high resolution EXPORT spectra 
is described. The EXPORT database is composed of 198 echelle spectra (see Chapter 2). 
38 of them have been used in this thesis (4 of BF Ori, 7 of SV Cep, 10 of UX Ori, 8 of 
WW Vul and 9 of XY Per). The reduction of echelle spectra is mostly a standardized 
process, which is described from Section 3.2 to Section 3.9. However, the normalization 
of the extracted spectra is not trivial. In Section 3.10, the normalization procedure ap- 
plied to the present set of data is reviewed. In Section 3.11, the quality of the reduced 
spectra and the limitations shown by the data, in order to be used in other researches 
different from this thesis, are discussed. 

3.1 Characteristics of an echelle spectrum 

In general, echelle diffraction gratings are used to obtain high resolution spectra R = 
A/AA ^ 50000 (Schroeder, 2000). These gratings are designed to operate in high inter- 
ference orders (m ~ 100). Thus high dispersion is obtained, but the free spectral range 
severely diminishes for each order (Amax - ^ mm ~ 100 A) and a spatial overlap occurs 
between contiguous orders. Special filters that only transmit the light of a spectral order 
can be used if only a small wavelength range, completely contained in an spectral order, 
is required. However, if a large spectral coverage is wanted, a low dispersive power op- 
tical element, in addition to the echelle grating, can be used. Thus the different orders 
can be spatially separated and a large portion of the spectrum (imax - ^min ^ 2000 A) can 
be recorded in a single CCD detector. 

The spectra used in this thesis were obtained with the UES echelle spectrograph 
placed on the WHT telescope (see Chapter 2). An optical diagram of UES is shown in 
Figure 3.1, taken from Unger et al. (1993). The echelle grating and the cross disperser 
element, composed of three different prisms, can be clearly identified. The instrument 
was set up to obtain spectra in the wavelength range 3800-5900 A with a resolution 
A/AA = 49 000. 49 diffraction orders (m = 95-153) of the grating were used. Each order 

o 

covers an approximate range of 100 A. The spectra were recorded on a square CCD 
detector of 2048 x 2048 pixels. 

The UES raw spectrum of a hot star is shown in Figure 3.2. The different spectral 
orders are shown as narrow approximately horizontal strips. The wavelength increases 
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Figure 3.1: Optical diagram of the UES echelle spectrograph, taken from Unger et al. 
(1993). The high dispersive power element (echelle diffraction grating: E) and the low 
dispersive power component (set of three different prisms: PI, P2 and P3) can be easily 
identified. 

from left to the right. The spectral direction is parallel to the x axis whereas the spatial 
direction runs through the y axis. The interference order increases from top to bottom. 
The orders with the highest interference number have the lowest wavelength, the highest 
width and the greatest inter-order spatial separation. 

The whole process of echelle reduction processes can be classified into two cate- 
gories. On the one hand, those corrections derived from the use of CCD detectors: bias, 
fiat-field and cosmic rays. On the other hand, those specific tasks of the echelle reduc- 
tion: Extraction of the orders, scattered light subtraction, wavelength calibration, and 
continuum normalization. 



3.2 Software 

The astronomical environment MIDAS has been used both for the reduction and the 
analysis of the data. This software package contains basic image manipulation com- 
mands and high level echelle reduction procedures (ESO, 1999a,b). The IRAF routine 
continuuin was also used during the normalization of the extracted spectra. In addi- 
tion, many programs have been developed by the author in the MIDAS programming 
language. They provided a certain degree of automation in the reduction (i.e. they con- 
stitute a pipeline) and rendered the tools to normalize the continuum. 
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Figure 3.2: Raw echelle spectrum of the star UX Ori. The spectrum is divided into 
several approximately horizontal strips. Each band corresponds to an interference order 
of the diffraction grating. The interference order increases from top (m = 95) to bottom 
(m = 153). For each individual order, the wavelength increases from left to the right. 

3.3 Bias correction 

The first step in the whole reduction procedure is the bias subtraction of the raw images. 
The bias is an electrical potential, which is added to that created by the electrons stored 

in each pixel of a CCD prior to their passage through the analog to digital converter. 
Thus it is assured that all the pixels will supply a positive number of counts, because the 
read-out noise may randomly generate an incorrect negative number of counts in pixels 
with very little or no stored photoelectrons. Two beneficial effects arise if negative counts 
are avoided. On the one hand, one encoding bit per pixel is saved. On the other hand, 
the possibility that a pixel with a large number of stored electrons displays wrongly a 
negative number of counts is prevented (Howell, 2000). 

If the temperature of the CCD is kept fixed during the whole night, the bias value 
should be constant. It is enough to take zero seconds exposures with the aperture closed 
to estimate the level of counts introduced by the bias. The bias frames are almost con- 
stant and their median is the bias value. This number must be subtracted from all the 
spectra before any other correction is performed, because the counts do not correspond 
to photons captured by the CCD. 

It was verified during the reduction that the bias does vary during the night. Perhaps 
it was due to problems in the thermal isolation of the CCD, or to the use of the maximum 
read-out speed available by the camera electronics. This effect is unimportant, since the 
dispersed light subtraction corrects these variations afterwards. 

The bias values were around 800 counts. It was observed that this value strongly 
depends on the CCD read-out speed. 
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Figure 3.3: Detection of echelle orders. A raw spectrum, similar to that of Figure 3.2 
is shown. The detected orders polynomial function, obtained from the fit of a standard 
star, is overplotted as a set of white lines. Each order is identified by a number, which 
is not equal to the interference order of the echelle grating. It can be seen that the 
orders detection is excellent, because it is almost impossible to distinguish the observed 
spectrum from the fit. 

3.4 Detection of echelle orders 

Each of the echelle orders is an independent spectrum which has to be reduced. Each 
order has a diff'erent curvature, width and position in the CCD, as it can be seen in 
Figure 3.2. The first step in the echelle reduction procedure is the characterization of all 
these parameters for each order. As they are nearly horizontal, the shape of the orders 
can be expressed via the following equation. 

y = f(x,m) (3.1) 

where y is the position of the order centroid in the ordinate axis for the abscissa x value. 
It has been verified that / can be properly approximated by low order polynomials (ESO, 
1999b). The polynomials used for the UES spectra are of 3'^'' degree in x and 4"^ degree 
in m. It was verified that the position of the orders did not show significant variations 
during each night. Therefore, it was enough to perform a single order detection for each 
observing night. The detection was made using very bright standard stars, which gave a 
good SNR in the whole set of orders. In Figure 3.3 a raw stellar spectrum is shown. The 
polynomial function obtained from the fit to a standard star is also displayed. It can be 
appreciated that the agreement is excellent. 

The width of the orders in the spatial direction is needed to extract correctly the 
spectrum. Vertical cuts were made in the observed spectra. It was verified that the orders 
have spatial profiles approximately gaussian in shape. In this way, the total width of an 
order is defined as 6 cr, being cr the dispersion value obtained from a gaussian fit to the 
spatial profile of the order. This value maximizes the SNR of the spectrum, because all 
the stellar light is extracted and the less possible amount of noise is added. 
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Figure 3.4: Dispersed light subtraction. Left: the figure shows a cross section of an 
echelle spectrum in the spatial direction. The orders can be identified as narrow profiles, 
approximately gaussian in shape, superimposed over a light continuum of smooth spatial 
variability. In order to remove the dispersed light, the continuum is approximated by a 
bidimensional cubic spline polynomial. This function is created from a suitable sampling 
of the inter-order space. Right: a cross section of the same echelle spectrum in the spatial 
direction after the dispersed light has been subtracted. 



As MIDAS uses the same extraction slit for all the spectral orders, the value adopted 
for the slit is the width of the bluest orders, which are the broadest ones. The width of 
the spatial profile strongly depends on the seeing, so it was necessary to compute (using 
a program) a different slit size for each spectrum. 



3.5 Scattered light subtraction 



The spectrographs have many imperfections that originate additional sources of illumi- 
nation on the CCD. Therefore, the light collected by the CCD is the sum of the astro- 
nomical spectrum and a scattered light continuum. Figure 3.4 shows a cross section of an 
echelle spectrum in the spatial direction. In the figure, the orders appear as narrow pro- 
files superimposed over a light continuum of smooth spatial variability. The dispersed 
light is not proportional to the intensity of the stellar spectrum, so it must be removed 
prior to the extraction of the spectrum. Otherwise, the equivalent width of all the spectral 
lines would be underestimated. 

The dispersed light presents a smooth variation across the CCD, so its intensity can 
be estimated for all the pixels using the values obtained from the inter-order space. This 
space is sampled at regular intervals, in order to create a dense network of data points. A 
point for every 100 pixels was used in the spectral direction of the inter-order space for 
the UES spectra. These points were used to generate a bidimensional cubic spline poly- 
nomial that reproduces accurately the dispersed light continuum. Once the dispersed 
light is estimated for every pixel, its subtraction is immediate. An example of this pro- 
cedure is displayed in Figure 3.4. 
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Figure 3.5: Removal of cosmic rays. Left: the figure shows the spectral region of the 
Ca II doublet for WW Vul. The spectrum has been reduced without applying the sub- 
traction of cosmic rays process. Right: the same spectrum is shown, but the algorithm 
for identification and removal of cosmic rays has been applied during the reduction pro- 
cedure. It can be verified that the process has subtracted all the cosmic rays but the 
absorption lines remain unaltered. 

3.6 Removal of cosmic rays 

The cosmic rays interact with the CCDs during the observations. When a cosmic ray 
hits a pixel, many photoelectrons are generated, because the particle is charged and it 
has a large amount of kinetic energy. Thus, the information stored in the pixels that have 
interacted with the cosmic ray (usually n ~ 1-5) is lost. If a spectrum polluted by cosmic 
rays is extracted, there randomly appear very narrow peaks (1-3 pixel wide), apparently 
in "emission" (see Figure 3.5). 

The high resolution of the echelle spectra makes the spatial profile of each order 
to slowly vary along the spectral direction. This effect can be used to identify the af- 
fected pixels and to reconstruct the stellar spectrum, to a certain extent. Thus, the spatial 
profiles along all the pixels in the spectral direction of each order are extracted, and an 
"observed" profile is obtained for each pixel along each order. For each pixel and each 
order, the median of the 30 nearest observed profiles along the spectral direction is com- 
puted, so two sets of "median" and "RMS error" (Root Mean Square, RMS) profiles 
can be obtained. The method for detecting and rejecting cosmic rays is a simple sigma- 
clipping peak rejection algorithm. If for any pixel in the observed profile, the diff'erence 
between the observed and the median values is greater than 4 times the RMS error (4 cr), 
then the observed value is replaced by the median value. 

This method fails when the fraction of pixels affected is so high that it is not possible 
to remove them from the median profiles. Thus a practical limit to the duration of the 
CCD integrations is imposed, because the number of cosmic rays grows linearly with 
time. The longest exposures taken lasted for 45 minutes and yielded fractions of bad 
pixels of about ~2%. 

It has been verified a posteriori that this method generally produces very good results. 
An example is shown in Figure 3.5. However, the algorithm has sometimes ruined the 
continuum of the spectrum, because periodic spurious absorptions have been introduced 
in the reddest orders. The width of these artificial absorptions is about ~5 A and the 
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Figure 3.6: Spatial profile of a flat-field (top) and an astronomical object (bottom). It 
can be appreciated that the orders of the astronomical object have an spatial profile ap- 
proximately gaussian in shape and they are contained inside those of the flat-field, whose 
shape is similar to a square wave. The origin of the difference is that the observed object 
is a point source (a star), while the flat-field is similar to the spectrum of an extended 
object. The flat-field is originated when the spectrograph flat-field entrance slit (which 
is larger than that of the astronomical objects) is uniformly illuminated with a tungsten 
lamp. 

depth is around ~5% of the continuum. 

3.7 Flat-field correction 

The flat-field correction has to be taken into account because of the variations of sensi- 
tivity between different pixels in a CCD. These differences are generally small. How- 
ever, dangerous defects can be present, for example cool or hot pixels, which present 
a substantially lower or higher quantum efficiency, respectively. Those bad pixels can 
be isolated or may belong to defective columns. These columns are specially damaging 
because they can destroy the information contained in a whole order. 

Flat-field exposures are used to identify and correct this undesirable effect. The name 
comes from the techniques developed for the observation and reduction of astronomical 
images. For this kind of observations it is possible to obtain images of fields uniformly 
illuminated (flat-fields), so any deviation from flatness can be attributed to imperfections 
in the CCD or the instrument (Howell, 2000). However, it is not possible to illuminate 
uniformly all the pixels in the detector while doing spectroscopy. The spectroscopic flat- 
field is the spectrum obtained when the spectrograph entrance slit is illuminated with an 
incandescent filament lamp (a tungsten lamp in our case). The entrance slit used for the 
flat-fields is larger than that used for the astronomical objects, the reason being to assure 
that the spatial profile of any order in any astronomical object is contained inside the 
spatial profile of the flat-fields (see Figure 3.6). 

The flat-field correction involves the division of the spectra of astronomical objects 
between the median of all the flat-fields obtained during a night. Thus, the pixel-to- 
pixel differences in sensitivity are eliminated However, the shape of the extracted orders 
is distorted, because the observed spectrum has been divided by the curvature function 
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of the flat-field. This effect can be removed if the stellar spectra are multiplied by the 
blaze function of the flat-field after the extraction of the spectrum. The curvature func- 
tion of the flat-field can be obtained by extracting the median flat-field as if it were an 
astronomical object. 

The MIDAS echelle routines perform the flat-field correction after the dispersed light 
has been subtracted. It has been confirmed a posteriori that this approach is wrong, 
because the photoelectrons generated by the scattered light in a pixel are proportional 
to its sensitivity. Because of this error in the software, the reduction was more difficult 
than expected. On the one hand, the detection of the orders was done twice (objects and 
flat-fields). On the other hand, the dispersed light had to be subtracted from the median 
flat-field, this being a CPU time consuming procedure. 

The flat-field correction could not be used for two observational campaigns, namely 
October 1998 and January 1999, because the slits used (objects and flat- fields) were not 
well aligned, so the spatial profiles of the stellar objects were not contained inside the 
flat-field spatial profiles. Despite this fact, it has been afterwards verified that the absence 
of flat-field correction has not significantly degraded the SNR of those spectra. 

3.8 Extraction of the spectra 

Once the corrections for bias, scattered light and cosmic rays have been applied, the spec- 
tra are ready for being extracted. The extraction is the obtention of a one-dimensional 
spectrum for each order of a bidimensional echelle spectrum recorded on the CCD. For 
any order and any position (pixel) in its spectral direction, the stellar light is distributed 
across the whole corresponding spatial profile. The extraction is the assignment to the 
pixel in that order of a value equal to the sum of all the counts stored across the spatial 
profile. The size of the profile is defined by a suitable extraction slit (see Section 3.4). 
It needs to be long enough to gather all the light from the object, but it does not have to 
be exceedingly large in order to minimize the noise. An example of spectrum extraction 
can be seen in Figure 3.7. 

3.9 Wavelength caUbration 

The wavelength calibration is a procedure that allows the extracted spectra to be ex- 
pressed in terms of the wavelength and not of the pixel number along the spectral direc- 
tion. The spectra produced by arc lamps filled with vapour of diverse atomic elements 
is used. These lamps generate a discrete emission spectrum composed of the lines cor- 
responding to the electron transitions of the atoms in the gas. A thorium-argon (Th-Ar) 
lamp was used for the reduction of the UES spectra. An example of arc lamp spectrum is 
displayed in Figure 3.8. The wavelength and the relative intensity of the strongest lines 
is tabulated for these lamps. Thus, if a sufficient number of lines is identified in each 
order (~20), the relation between the pixel numbers and the wavelengths can be found. 
In general, it is enough to use low order polynomial fits (n < 4 in our case). 

The wavelength calibration is fully integrated in MIDAS. Thus it is only needed to 
perform interactively a few identifications and the program is then capable to identify 
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Figure 3.7: Extraction of a spectrum. If the extraction procedure is applied to the bidi- 
mensional spectrum on the left, the spectrum on the right is obtained. The extracted 
spectrum is composed of 59 rows, corresponding to each of the orders of the spectrum 
on the left. 




Figure 3.8: Th-Ar arc lamp spectrum. It can be seen that it is a discrete spectrum, 
composed of a huge amount of emission lines whose wavelengths and relative intensities 
are tabulated for the strongest transitions. 
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iteratively a sufficient number of lines for each order without any further input. In gen- 
eral, the calibrations are excellent. The typical RMS errors are ~20 mA (0.39 pixels), 
i.e. 5 times lower than the spectral resolution (~100 mA). 

The Th-Ar lamp got damaged during one observing night (January 29* 1999) and 
the reserve iron-argon lamp was used instead. The available identification tables for the 
Fe-Ar lamps did not have a sufficiently large number of identifications. A new table was 
created from the existing tables of different lamps, included in the MIDAS and IRAF 
packages. The calibration results were acceptable, but worst than those obtained with 
the Th-Ar lamp. 

3.10 Continuum normalization 

The wavelength calibrated spectra have a continuum defined by the curvature function 
of the orders, which is the product of the spectrograph response function by the spec- 
tral distribution of the star. If the curvature function is obtained, the observed spectrum 
can be normalized dividing it by the aforementioned function. The continuum normal- 
ization is the last step in the reduction procedure carried out in this thesis, because a 
spectrophotometric calibration was not needed. 

The continuum of a high SNR spectrum obtained from a low rotation velocity (v sin / 
< 50 kms~^) F spectral type (or earlier) star is quite easy to determine. The higher the 
temperature, the fewer and shallower are the photospheric lines. Moreover, low rotation 
velocities generate less overlaps between neighbouring lines. If the previous conditions 
are given, the characterization of the envelope function of each extracted order can be 
easily performed via a moderate order polynomial (~4-7). However, it is impossible 
to fit properly the orders containing the B aimer series lines, because these lines are as 
broad as the orders. The Balmer series is essential for this thesis, so it was made the 
approximation that the curvature function in an order containing a Balmer line is equal to 
the logarithmic average of the curvatures in the contiguous orders. This approximation, 
being inexact, produces reasonable results, because it has been verified that the curvature 
changes slowly from order to order. 

Only a few of the stars in the whole sample strictly comply with all the stated re- 
quirements, specially that about the spectral type, because the EXPORT observations 
cover the full range B9-K7. However, it is enough to determine the curvature function 
for a single (normalization standard) star each night. The curvature function obtained is 
used to normalize the remaining stars. If this procedure is performed to a (target) star of 
the same spectral type as the normalization standard star and the seeing of both obser- 
vations is the same, then a target star spectrum normalized to the quotient between the 
fluxes of the target and the normalization standard stars is obtained. If the target star has 
a diff'erent spectral type or it has been observed under diff'erent seeing conditions, orders 
without curvature will be obtained, but they will not be normalized to any fixed quantity, 
because the quotient of the fluxes depends now on the wavelength. 

Once the spectrum of a single non- variable star has been normalized to unity, it can 
be used to determine the curvature function of the orders of any other spectrum of the 
same star. It is enough to take into account the changes in the radial velocity of the 
normalized spectrum before dividing it by the observed spectrum. Thus the curvature 
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function calculation can be "propagated" towards other observing nights and the number 
of polynomial fits needed is minimized. Moreover, once the spectral orders of the target 
stars are normalized by means of the curvature function of the normalization standard 
star, they can be merged. Then the complete normalized spectrum of any object can be 
generated and new normalization standard star spectra can be generated. 

3.11 Validity of the reduction method 

It has been verified that the extraction and wavelength calibration of the spectra are 
excellent. It has also been verified that the subtraction of cosmic rays has degraded the 
signal to noise ratio of a few stars to levels about <20. Fortunately, the stars in this thesis 
have not been significantly affected by this effect. 

The reduction method used in this thesis has been compared with those used by 
two astronomers experienced in the use of echelle spectra: Dr. Cameron (University of 
St. Andrews, United Kingdom) and Dr. Grinin (Crimean Astrophysical Observatory, 
Ukraine and Pulkovo Astronomical Observatory, Russia). The spectra reduced for this 
thesis are similar in quality to those obtained by Dr. Cameron for his research activities. 
However, it became apparent that the normalization of our spectra is worst than that 
achieved by the group of Dr. Grinin for a series of echelle spectra of UXOR stars (Grinin 
et al., 2001). Two are the reasons of this difference. On the one hand, this group has a 
collection of, properly normalized, intermediate resolution spectra obtained in the same 
wavelength range than the echelle spectra for all the stars studied. On the other hand, 
they heavily rely on synthetic spectra to adequately normalize the Balmer lines by means 
of an iterative fit process. Even though the use of artificial spectra from the very moment 
of the reduction may introduce an additional error source, the goodness of the results 
obtained justify the validity of the procedure. 

The reduced spectra can be used without reservation for any purpose which requires 
small spectral ranges (AA :$ ± 10 A, Av ^ + 600 kms~ ). However, the continuum 
is not accurately defined if wider wavelength ranges are used. This was made clear in 
an attempt made by Dr. Montesinos to use the UES spectra to determine stellar effec- 
tive temperatures and surface gravities using the Balmer lines (AA ^ + 100 A). It was 
verified that such large intervals led to erroneous results. In this thesis no TAC with 
radial velocity greater than 500 kms~^ has been analysed, so the validity of the spectral 
analysis, performed in Chapters 5 and 6, is guaranteed. 



Chapter 4 



Calculation of projected rotational 
velocities via the Fourier transform of 
line profiles 

The fit of stellar parameters is an iterative process which will be explained thoroughly 
in Sections 5.3.1 and 6.3.1. This procedure needs some initial estimates for the follow- 
ing stellar parameters: effective temperature (Teff), surface gravity (logg) and rotational 
velocity projected on the line of sight (vsin/). These values have been obtained from 
the spectral types and v sin / rotational velocities derived by Mora et al. (2001) for the 
whole EXPORT sample. The initial values of v sin i can be directly used in the fit of 
stellar parameters. However, in order to obtain Tefr and logg from the spectral type, 
some empirical calibrations are needed. In this thesis, the relations obtained by de Jager 
& Nieuwenhuijzen (1987) and Habets & Heintze (1981) have been used. In this Chap- 
ter, the method employed by Mora et al. (2001) for the measurements of v sin / will be 
described in detail. 



4.1 Theoretical basis 

Many methods for the determination of stellar rotational velocities have been developed 
so far. The most important procedures are briefly described by Gray (1992). Some 
are based in spectroscopic measurements. Two of them are particularly suited for their 
application to echelle spectra: simultaneous fitting of v sin / and other stellar parameters 
via synthetic spectra and the use of the Fourier Transform of isolated line profiles. The 
first method is extensively used in Section 6.3.1, whereas the second one has been used 
by Mora et al. (2001) (strictly speaking, the procedure used is somewhere in between 
both generic methods, because a small amount of synthetic spectra is also utilized). 

As it can be seen in Gray (1992), the first spectroscopic determinations of rotational 
velocities were performed by visual comparison of the spectra of target stars with some 
template spectra of rotational standard stars, which were grouped in catalogues (e.g. 
Slettebak et al., 1975). A more advanced technique is based on the use of empirical 
relations between the Full Width at Half Maximum (FWHM) of some lines and v sin /. 
An example of such relations is given by Gray (1992) from the data of Slettebak et al. 
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(1975). These calibrations give precisions in the measurements around 10-20%. 

More advanced procedures came from the use of the whole line absorption profile 
(D(A) = Fy/Fc), and not only the FWHM. In this way, the Fourier transform of the 
profile {d{(T)) can be computed and the spectrum changes from the wavelength space 
(/I, A) to the spatial frequency domain (cr, A~^). In what follows, the fundamentals of 
the simplest, although very efficient, method of measurement of v sin / from the Fourier 
Transform will be described. This method is very precise. Errors as low as 5% can be 
achieved in favourable cases. 

The main hypothesis of the method is that the star rotates as a rigid body and the 
emitted spectrum H{X) is the same for all points in the stellar disk. Therefore, the ob- 
served integrated spectrum is found to be the spectrum emitted by any part of the disk, 
H{A), convolved with a certain rotational profile, G{A): 



D(A)= f H(A-AA)G(AA)dAA (4.1) 

U —CO 

A linear law of limb darkening is also postulated. This law relates the intensity 
emitted from the disk surface element exactly located under the line of sight (7°) with 
the light emitted from any other element (Ic) whose radius vector R subtends an angle 9 
with the observer's direction. 

/J7° = 1 -6(1 -cos^) (4.2) 
Thus, an analytical formula for the rotational profile: G{Au,^ + A/l) can be found: 

2(1 - 6)[1 - {AAlAAif]^ + \n6{\ - {AAI AAif^ 
G{AA) = ' (4.3) 

nAAiil - 6/3) 

where AAi = /iune x v sin i/c. Carroll (1933) found the following approximate expression 
for the line profile Fourier transform g(cr), assuming a limb darkening coefficient e = 0.6. 

, , Ji(x) 3cosx 3sinx 

where x = IncrAnneV sini/c. This function is the sum of three oscillatory components. 
The first component contains a first-order Bessel function. Being an oscillating fuction, 
the profile function has several zeroes. The first one (ctj) is located at: 

0.660 c 

o-i = — (4.5) 

/tlineV Sm I 

It is known that the transform of the observed spectrum d(<T) is the product of the 
transform of the profile g(cr) by the transform of the spectrum emitted by every infinites- 
imal region of the stellar disk /i(cr): 

d(cr) = h(cr)Xg((r) (4.6) 

So every frequency cr which makes the rotational profile g{o-) zero will also make 
any line profile transform d{cr) to vanish. If the zeroes of h{cr) correspond to higher 
frequencies than those of g(cr), it would be enough to measure the first zero of J(cr) to 
obtain a value of v sin /. If equation 4.5 is inverted: 
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vsmi = 



0.660 c 

'tlineC 1 



(4.7) 



4.2 Operational procedure 



According to the previous Section, the rotational velocity can be immediately computed 
if the first zero of the Fourier transform of a suitable line profile is known. Actually, 
the line profiles are not only affected by rotation but by many other factors, such as 
the presence of nearby faint lines, low SNR, defective reduction, etc. Therefore, the 
transform does not completely vanish, so instead of zeroes we will have minima in the 
power spectrum (the modulus of the transform). The higher the quality of the selected 
lines and the spectra, the deeper the minima will be. The method consists of three 
steps: search for strong isolated lines whose main broadening mechanism is rotation, 
normalization of the line profiles and computation of the Fourier transforms. Synthetic 
spectra have been used to perform properly the first two stages in the method, as it will 
be seen below. 

It is impossible to find true isolated lines in a stellar spectrum, due to the huge amount 
of actual optical transitions. It will be assumed that an isolated line is that with "higher 
intensity" than those in the wavelength range where the line exists, i.e. the neighbour- 
hood of the line absorption profile, which will be rotationally broadened. It can be easily 
determined from photospheric zero rotational velocity spectra if any line is of "higher 
intensity". First, a synthetic spectrum with a temperature and gravity similar to those in 
the star should be generated. The spectra have been calculated using the stellar model 
atmospheres given by Kurucz (1993a) and the code SYNTHE (Kurucz, 1993b). Since 
the depth of the lines does not vary quickly with changes in the stellar parameters, very 
accurate initial values for Tgff and logg are not needed. It is enough to estimate them 
from the spectral types. Initial guesses of the spectral types were obtained from the 
astrophysical database Simbad (http : //simbad . u-strasbg . fr). 

The Fourier transform of both the observed absorption profile (in order to measure 
V sin /) and the non-broadened synthetic spectrum were evaluated in the same wavelength 
range. If the first minimum of both transforms match, then the value obtained for v sin / 
is wrong, because the frequency used cti is not related to rotation but to the distance 
between the line studied and another one of significant intensity. On the other hand, if the 
first zero of the observed profile happens at a lower frequency than that of the synthetic 
spectrum, it can be undoubtedly affirmed that the target line is of "higher intensity" than 
those in its neighbourhood. Two examples of line selection can be seen in Figure 4.1. 
One of the examples is correct, while the another one is wrong. 

Once the right lines are selected, the profiles must be properly normalized. In prin- 
ciple it should be enough to fix the left and right edges to unity. However, if the rotation 
velocity is high (v sin / ^ 100 km s~^), the lines will generally be broad and shallow. This 
favours neighbouring lines to blend, so the continuum will not be given by the apparent 
edges of the lines. If the spectrum is normalized according to the procedure described, 
the equivalent width of the line will be underestimated and the determination of the rota- 
tion velocity will be wrong. Synthetic spectra can also be used to minimize this problem. 
If the borders of the lines are not set to unity but fitted to a continuum given by an arti- 
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Figure 4.1: Line selection for vsinz estimations. Left: two portions of the spectrum 
of BF Ori are shown (solid black line). The top spectrum was obtained in JD 1209.55 
and the bottom one in JD 1112.63. Synthetic spectra generated using Teff = 8750 K, 
log g = 3.5 and v sin i = are overplotted (blue dashed line). The upper and lower bounds 
of the interval used for the Fourier transform computations are shown by two dotted 
lines. The ions responsible for the most important electron transitions in the interval are 
also shown in each graph. Right: next to each line selection graph, the power spectrum 
of both the observed profile (black solid line) and the synthetic profile (blue dashed line) 
are displayed. The placement of the first minimum in the observed profile is shown in 
each graph as a dotted vertical line. It can be appreciated that the selection made in the 
upper left panel is correct, because the first minimum of the Fourier transform of the 
synthetic spectrum is located at a higher frequency. However, the selection made in the 
lower left graphic is wrong because both minima match in frequency. 
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Figure 4.2: Continuum placement for v sin i calculations. A portion of the spectrum of 
XY Per, obtained in JD 1 1 13.55, is shown as a black solid line. The blue dashed line is 
a synthetic spectrum used to place the continuum of the observed spectrum. The same 
spectrum but without broadening is shown in green. It can be seen that most of the 
observed broad absorption component, which will be used for the v sin / determination, 
is mainly composed by a strong broadened line of Mg i. The dotted lines indicate the 
interval used to calculate the power spectrum. If the apparent edges of the observed line 
had been used to normalize the line profile, the equivalent width of the line would have 
been underestimated by about 50%. 

ficial spectrum of the star, then a more reasonable value for the equivalent width of the 
line will be used. The artificial spectra required are those used in the line selection pro- 
cess. The value of vsinz is needed to broaden the synthetic spectra, but this magnitude 
is exactly what is being looked for. Once again, very precise values are not needed, so 
V sin i data previously published by other authors or, when not available, the typical val- 
ues relative to the spectral types, have been used. The procedure described is explained 
in Figure 4.2. 

Fourier transforms were computed by means of the MIDAS routine FFT/POWER, 
which provides the real part, the imaginary part and the power spectrum of the transform. 
Once the power spectrum is obtained, the determination of its first minimum location 
and, therefore, v sin i are trivial. 

The rotational velocity has been measured for the largest possible number of avail- 
able lines in the spectra of each star. Thus the average values and statistical errors finally 
presented by Mora et al. (2001) have been calculated. In general, the statistical errors 
obtained are about 10-20 km s"^ 

4.3 Uncertainties in the method 

There are more precise and complete methods for the determination of the rotational 
velocity based on the Fourier transform (Gray, 1992). These alternative methods do 
not only use the position of the first minimum, but the whole shape of the function. 
Therefore, additional information can be obtained, as the value of e, the identification 
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of rotation or macroturbulence as the dominant line broadening mechanism, and even a 
measurement of the differential rotation in the star. The great advantage of using only 
the position of the first minimum for the determination of v sin i relies in its absolute 
independence from theoretical models (e.g. the variation of e with the temperature and 
gravity) and experimental measurements (e.g. the use of rotational standard stars or 
instrumental and thermal broadening profiles). 

The assistance of synthetic spectra for the selection of lines and continuum fitting has 
allowed the method to be applied to very high rotational velocity (v sin i ^ 200 km s"^) 
stars. This is undoubtedly an advantage (4 out of 5 stars in this thesis have v sin i > 
200 kms~^). However, the results obtained are less model-independent, because errors 
on the computation of the artificial spectra have been introduced in the final values ob- 
tained. 

The resolution of the observations places a lower limit to the range of measurable ve- 
locities. This can be expressed in terms of the Fourier transform. The 6 km s^ resolution 
of the UES spectra can be represented by a narrow instrumental profile which is con- 
volved with the stellar spectrum. This instrumental profile is a function, approximately 
gaussian in shape, with a FWHM of 6 km s ' and a minimum in its power spectrum lo- 
cated at the frequency corresponding to this velocity. The observed line profile transform 
is the product of the three individual transforms (rotational profile, instrumental profile 
and spectrum emitted by every region of the stellar disk). Thus, another minimum in the 
power spectrum will be found exclusively due to the instrumental profile, so, if v sin / is 
lower than the resolution, the first minimum of the transform will not be generated by 
rotation but by the instrumental profile. 

The limb darkening law used is another source of systematic errors. The method 
relies in the use of a linear function with e = 0.6. This approximation is valid for the 
Sun and stars of similar spectral types. Many authors have questioned the validity of this 
empirical relation. Diaz-Cordoves & Gimenez (1992) demonstrated that the difference 
between a linear law and other models of limb darkening is very small for the wavelength 
ranges explored in this thesis, at most ~5% in the less favourable conditions. However, 
the specific value of e used in the v sin i measurements does matter. 

The influence of 6 will be maximum when the rotational profile fully dominates the 
observed line shape. This happens when v sin / is much greater than both the resolution 
of the observations (6 kms ') and the macroturbulence (~2 kms"^). Solano & Fernley 
(1997) estimated the error associated with the selection of a wrong limb darkening co- 
efficient. We have improved this method as follows. A grid of rotational profiles with 
vsin/ = 100 kms~^ and e values ranging from 0.0 to 1.0 with a step of 0.05 has been 
generated. A set of artificial spectral lines has been generated from these profiles with 
an arbitrary central wavelength of 5000 A. The line shape is exclusively regulated by 
the rotational profile, v sin i has been calculated for these artificial lines with the method 
described in this section (i.e. assuming e = 0.6). Thus the errors in the measurement of 
V sin / arising from a wrong value of e have been estimated and are shown in Figure 4.3. 
Diaz-Cordoves et al. (1995) provide theoretical predictions for the limb darkening coef- 
ficients. It can be verified that the maximum range of variation for the whole EXPORT 
sample is e = 0.35-0.95. The maximum error will be of 7% for the cooler stars (K spec- 
tral type). If only the five stars studied in this thesis are considered, the limb darkening 
coefficients will be restricted to the range e = 0.50-0.65, which implies a maximum error 
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Figure 4.3: Error made in the the calculation of v sin i according to the linear limb dark- 
ening coefficient of the star. In the figure, the maximum error made when vsin/ is 
measured for stars with limb darkening coefficient e ^ 0.6 is shown. This error is negli- 
gible (< 1.6%, 6 = 0.50-0.65) for the stars studied in this thesis, but could be significant 
(< 7 %) for the coolest objects in the EXPORT sample, which present the highest e 
values. 

of 1.6%. It can be concluded that the limb darkening law used is not a significant source 
of error, except for, possibly, the cooler stars in the EXPORT sample. 

4.4 Comparison with values obtained by spectral syn- 
thesis 

Table 4.1 shows the values of vsin / calculated via the Fourier transform and those ob- 
tained in Section 6.3.1 by means of synthetic spectral fits to photospheric lines. The 
values relative to the star HD 163296 have also been included, because its stellar pa- 
rameters were also calculated by the spectral synthesis method. That star was initially 
included in this thesis, but its strong stellar winds rendered the TACs deconvolution val- 
ues unusable. It can be seen that the agreement between both completely independent 
methods is remarkable for the whole set of stars. 
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Table 4. 1 : Comparison between two v sin / measurement methods. The values of v sin / 
for the stars in this thesis are shown in the table. Column 2 includes the results obtained 
searching for the first minimum in the power spectrum. Column 3 gives the values 
retrieved from the simultaneous fit of Tgff , log g and v sin / by means of spectral synthesis. 
The errors of these measurements are discussed in Section 6.3.1. The rotational velocity 
was not determined by both methods for UX Ori, because the value obtained from the 
Fourier transform provided an excellent agreement between the observed and synthetic 
spectra. 

fThe values obtained for HD 163296, not included in this thesis, are given for compari- 
son purposes only. 



Star 



Fourier transform Spectral synthesis 



BFOri 
SV Cep 
UXOri 
WWVul 
XYPer 



37 + 2 37 
206 + 13 225 
215 + 15 



220 + 22 210 
217 + 13 200 
133 ± 6 125 



HD 163296^^ 
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Abstract 

We present the results of a high spectral resolution (i/A/l = 49000) study of the cir- 
cumstellar (CS) gas around the intermediate mass, pre-main sequence star UX Ori. The 
results are based on a set of 10 echelle spectra covering the spectral range 3800 - 5900 
A, monitoring the star on time scales of months, days and hours. A large number of 
transient blueshifted and redshifted absorption features are detected in the B aimer and 
in many metallic lines. A multigaussian fit is applied to determine for each transient 
absorption the velocity, v, dispersion velocity, Av, and the parameter R, which provides 
a measure of the absorption strength of the CS gas. The time evolution of those param- 
eters is presented and discussed. A comparison of intensity ratios among the transient 



49 



50 



Chapter 5. A dynamical study of the circumstellar gas in UX Orionis 



absorptions suggests a solar-like composition of the CS gas. This confirms previous re- 
sults and excludes a very metal-rich environment as the cause of the transient features in 
UX Ori. The features can be grouped by their similar velocities into 24 groups, of which 
17 are redshifted and 7 blueshifted. An analysis of the velocity of the groups allows us 
to identify them as signatures of the dynamical evolution of 7 clumps of gas, of which 4 
represent accretion events and 3 outflow events. Most of the events decelerate at a rate 
of tenths of ms"^, while 2 events accelerate at approximately the same rate; one event 
is seen experiencing both an acceleration and a deceleration phase and lasts for a period 
of few days. This time scale seems to be the typical duration of outflowing and infalling 
events in UX Ori. The dispersion velocity and the relative absorption strength of the 
features do not show drastic changes during the lifetime of the events, which suggests 
they are gaseous blobs preserving their geometrical and physical identity. These data 
are a very useful tool for constraining and validating theoretical models of the chemi- 
cal and physical conditions of CS gas around young stars; in particular, we suggest that 
the simultaneous presence of infalling and outflowing gas should be investigated in the 
context of detailed magnetospheric accretion models, similar to those proposed for the 
lower mass T Tauri stars. 

Key words. Stars: formation - Stars: pre-main sequence - Circumstellar matter - Ac- 
cretion, accretion disks - Lines: profiles - Stars: individual: UX Ori 

5.1 Introduction 

The detection of planetesimals is highly relevant for the study of the formation and evolu- 
tion of planetary systems, since it is nowadays accepted that planets form from CS disks 
via the formation of planetesimals (Beckwith et al., 2000). Several lines of evidence sug- 
gest that the young main sequence A5V star Pic (20 Myr, Barrado y Navascues et al., 
1999) hosts planetesimals inside its large CS disk. The main argument is the presence of 
transient Redshifted Absorption Components (RACs) in high resolution spectra of strong 
metallic lines, like Can K 3934 A. These spectroscopic events have been interpreted as 
being caused by the evaporation of comet-like highly hydrogen-depleted bodies. The 
interpretation is known as the Falling Evaporating Bodies (FEB) scenario (Lagrange 
et al., 2000, and references therein). However, the presence or absence of planetesimals 
during the pre-main sequence (PMS) phase is a controversial observational topic. The 
time scale for the formation of planetesimals (~10'^ yr, Beckwith et al., 2000) is shorter 
than the duration of the PMS phase (~1-10 Myr), which suggests that they should exist 
during PMS stellar evolution. 

UX Ori-like PMS objects (UXORs) are characterized by a peculiar photo-polarimetric 
variability, which has been interpreted as the signature of massive, almost edge-on, 
CS disks (Grinin et al., 1991). Most UXORs have A spectral types and therefore are 
the PMS evolutionary precursors of [3 Pic. UXORs have been reported to show RACs 
(Grinin et al., 1994; de Winter et al., 1999); they also show Blueshifted Absorption Com- 
ponents (BACs) in their spectra. In this paper, the acronym TAC (Transient Absorption 
Component) will be used to denote both RACs and BACs. In analogy to [3 Pic, RACs 
observed in UXORs have been interpreted in terms of the FEB scenario (Grady et al.. 
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2000). However, this interpretation has recently been questioned by Natta et al. (2000b), 
who used the spectra obtained by Grinin et al. (2001) to analyze the dynamics and chem- 
ical composition of a very strong, redshifted event in UX Ori itself, an A4 IV star (Mora 
et al., 2001) - 2x 10^ year old (Natta et al., 1999). Gas accretion from a CS disk was 
suggested in Natta et al. (2000b) as an alternative to the FEB scenario to explain the 
observed RACs in UX Ori. In addition, Beust et al. (2001) have found that the FEB hy- 
pothesis cannot produce detectable transient absorptions in typical HAe CS conditions, 
unless the stars are relatively old (> 10^ yr). 

A detailed observational study of the kinematics and chemistry of TACs, e.g. by 
means of echelle spectra which simultaneously record many metallic and hydrogen lines, 
can discriminate between the two scenarios. For instance, in a FEB event large metal- 
licities are expected, while gas accreted from a CS disk would have approximately solar 
abundances. A strong observational requirement is set by the time scale of monitoring 
of the TACs. In Natta et al. (2000b) spectra were taken 3 days apart. Since UX Ori is a 
highly variable star, there is some ambiguity in identifying transient spectral components 
of different velocities detected over this time interval as having the same physical origin. 
A better time resolution is needed in order to ensure that the TACs observed at different 
velocities are due to the dynamical evolution of the same gas. The EXPORT collab- 
oration (Eiroa et al., 2000b) obtained high resolution echelle spectra of a large sample 
of PMS stars (Mora et al., 2001). About 10 PMS stars showed TACs which were inten- 
sively monitored (At < 1 day). The study of the kinematics and chemistry of these events 
provides important tools for identifying their origin, or at least to put severe constraints 
on it. 

This paper presents an analysis of the TACs observed in UX Ori by EXPORT and 
shows that the results are not compatible with a FEB scenario. The layout of the paper 
is as follows: Section 5.2 presents a brief description of the EXPORT observations. 
Section 5.3 presents the procedures followed in the analysis of the spectra. Section 5.4 
presents the kinematic and chemical results of the detected TACs. Section 5.5 gives 
a discussion of the dynamics and nature of the gas. In Section 5.6, we present our 
conclusions. 

5.2 Observations 

High resolution spectra of UX Ori were taken in October 1998 and January 1999 using 
the Utrecht Echelle Spectrograph (UES) at the 4.2m WHT (La Palma Observatory). We 

o 

collected 10 echelle spectra in the wavelength range 3800 - 5900 A, with resolution 
/1/A/l = 49000 (6 km/s); exposure times range between 20 and 30 minutes. The ob- 
serving log is given in Table 5. 1, in which the long-term (months) and short-term (days, 
hours) monitoring is evident. Standard MIDAS and IRAF procedures have been used for 
the spectroscopic reduction (Mora et al., 2001). Final typical signal to noise ratio (SNR) 
values are ~150. 

EXPORT also carried out simultaneous intermediate resolution spectroscopy, optical 
photo-polarimetry and near-IR photometry of UX Ori during the nights when the UES 
spectra were taken (Mora et al., 2001; Oudmaijer et al., 2001; Eiroa et al., 2001). Ta- 
ble 5.1 also gives photo-polarimetric data taken simultaneously with the UES spectra; 
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Table 5.1: EXPORT UES/WHT UX Ori observing log. The Julian date (-2450000) of 
each exposure is given in column 1. Columns 2 to 5 give simultaneous V, H, K and %Py 
photopolarimetric data, when available, taken from Oudmaijer et al. (2001) and Eiroa 
et al. (2001). 



Julian date 


V 




H 


K 


1112.5800 


9.94 


1.33 


8.19 


7.37 


1113.6034 


9.82 


1.32 


8.20 


7.37 


1113.7194 










1207.5268 










1208.5072 






8.41 


7.56 


1209.4204 


9.80 


1.15 


8.31 


7.62 


1210.3317 


9.89 


1.26 






1210.3568 


9.86 


1.34 






1210.4237 


9.88 


1.30 


8.34 




1210.5156 


9.90 


1.31 


8.34 





they show that UX Ori was always in its bright state with little variation, but significant 
(Oudmaijer etal., 2001). 

5.3 Analysis of the spectra 

A detailed analysis of TACs requires high SNR (> 25) and high spectral resolution 
(< 10 km/s) in order to resolve the kinematic components observed simultaneously in 
diff'erent lines, as is shown in Fig. 1. This figure presents the Nai D lines in one of the 
UES spectra together with the simultaneous intermediate resolution spectrum (R ~ 6000) 
in the same lines. The UES spectrum clearly shows a sharp, IS, narrow zero velocity 
component and a RAC in both lines of the doublet while both components, though visi- 
ble, cannot be cleanly separated in the lower resolution spectrum. This is why we restrict 
our further analysis of the absorptions to the UES spectra. 

All the UX Ori UES spectra show CS spectral features in a variety of lines. The fea- 
tures are seen in absorption and can be either blueshifted or redshifted; some underlying 
emission is also present (see below). Their profiles are complex and blended compo- 
nents are directly seen in many cases. Thus, the analysis of the CS contribution to the 
observed spectra requires a careful subtraction of the stellar photospheric spectrum and a 
method to characterize the blended kinematical components. In this section we describe 
the procedure we have followed, which is the one used in Natta et al. (2000b), although 
we have improved it by considering all detected TACs simultaneously. 

5.3.1 Subtraction of the UX Ori photospheric spectrum 

Firstly we estimate the radial velocity of the star. The heliocentric correction of each 
observed spectrum is computed using MIDAS. Then, the radial velocity of the star is 
estimated by fitting the Na I D sharp IS absorption. This absorption is stable in velocity 
(velocity differences in our spectra are less than 1.0 km/s) and is related to the radial 
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Figure 5.1: Comparison between intermediate (top, J.D. 2451210.3098, 10 min. ex- 
posure time) and high (bottom, J.D. 2451210.3317, 30 min. exposure time) resolution 
spectra. The Na i D spectral region is shown. A sharp IS zero velocity component and 
a much broader RAC (indicated in the figure) can be seen in both lines of the doublet in 
the high resolution spectrum. They are also detected in the intermediate resolution spec- 
trum, but the big difference in the width of the components is not clear. For the sake of 
clarity an artificial vertical displacement of 0.2 units has been applied to the intermediate 
resolution spectrum. 

velocity of the star (Finkenzeller & Jankovics, 1984). A value of 18.3 + 1.0 km/s is 
obtained, in good agreement with the 18 km/s estimate by Grinin et al. (1994). 

Kurucz (1993) model atmospheres assuming solar metallicity and turbulence veloc- 
ity of 2 km/s have been used to synthesize the photospheric spectra, which are later used 
for a comparison with the observed ones. Atomic line data have been obtained from the 
VALD online database (Kupka et al., 1999), and the codes SYNTHE (Kurucz, 1979; Jef- 
fery et al., 1996) and SYNSPEC (Hubeny et al., 1995) are used to synthesize the metallic 
lines and Balmer hydrogen lines, respectively. In this way we compute a large grid of 
synthetic spectra with the effective temperature, gravity and rotation velocity as free 
variable parameters. 

The best set of free parameters is estimated by comparing a number of lines among 
the observed spectra and the synthetic ones. This is not an easy task since most UX Ori 
lines are variable to some extent, and the choice of pure photospheric lines is not trivial. 
Our choice is to consider a number of weak lines that show a minimum degree of vari- 
ability (less that 1%) in the UES spectra. The blends located at 4172 - 4179 A (mainly 
Fe II and Ti ii) and 4203 A (Fe i) have been used to estimate Tgff (effective temperature) 
and log g (logarithm of the surface gravity), since they are sensitive to changes in Tgff 
and logg for early type A stars (Gray & Garrison, 1987, 1989). The synthetic spectra 
have also been broadened allowing for stellar rotation. The best set of free parameters 
we find is Tgff = 9250 K, logg = 4.0, v sin i = 215 km/s. Fig. 5.2 shows the comparison 
between the synthetic and the average observed spectrum in the above-mentioned lines; 
the agreement is excellent. The stellar parameters of the synthetic spectrum are in very 
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Figure 5.2: Comparison of the observed average spectrum of UX Ori (solid line) 
to the synthetic spectrum with stellar parameters Teff = 9250 K, logg = 4.0 and 
V sin / = 215 km/s (dashed line). The spectral region shown includes the blends at 4172 - 
4179 A (mainly Fe n and Ti n) and 4203 A (Fe i). 

good agreement with the spectral type and rotational velocity derived elsewhere (Mora 
et al, 2001). 

The next step consists of the subtraction of the best synthetic photospheric spectrum 
from each observed spectrum. This has been carried out by fitting the continuum with 
a linear law between two points selected at both sides of those lines showing TACs. As 
an example Fig. 5.3 shows the Can K line results, where the synthetic spectrum and 
the J.D. 2451209.4204 UX Ori observed one are plotted together with the residual after 
the subtraction. The residual corresponds to the CS contribution to the observed spec- 
tra. As a more convenient way to represent this contribution we define the R parameter 
(normalized residual absorption): 

= 1 - FobJFsyn 

R = means no CS absorption, and R = \ denotes complete stellar light occultation, 
i.e. R quantifies the CS absorption strength and could be directly compared to the models 
in Natta et al. (2000b). R is also plotted in Fig. 5.3. 

5.3.2 R multicomponent gaussian fitting 

The CS contribution presents a very complex line profile caused by the blending of 
diff"erent transient features. As an example. Fig. 5.3 shows the components of the Can K 
line which are clearly present in that particular observation. For the sake of analysing 
the profile we assume that it is caused by gas with difi'erent kinematics and that each 
component can be represented by means of a gaussian function. This choice provides a 
direct interpretation of the fit parameters in terms of a physical characterization of the 
gaseous transient absorption components: the gaussian center gives the central velocity 
(v) of each component, the gaussian FWHM (Full Width at Half Maximum) provides 
the velocity dispersion (Av) and finally the gaussian peak gives the value of the residual 
absorption R. Since the features are blended a multigaussian fit is needed, which has 
been carried out by means of the IRAF ngaussfit routine. Fig. 5.3 shows the results of 
the multigaussian fit for the Ca ii K line. 4 gaussians, i.e. 4 kinematical components, each 
with their corresponding parameters, provide an excellent agreement with the residual 
CS contribution profile. 
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Figure 5.3: Subtraction of the photospheric spectrum. The spectral range corresponds 
to the Can K line. Top: observed (solid line) and synthetic (dotted line) spectra are 
displayed. Middle: Residual due to the CS contribution (observed - photospheric). Bot- 
tom: Normalized residual absorption, R = (1 - observed / photospheric). Vertical tick 
marks indicate the 4 absorption components of the Ca n K which can be identified in this 
spectrum, including the narrow IS absorption. In the bottom we show the 4 gaussians 
(dashed lines), whose sum (dotted line) excellently reproduces the R profile (solid line). 
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This analysis is applied to the lines with the strongest variability, which are good 
tracers of the CS gas properties: Ufi 4861 A, Hy 4340 A, U6 4102 A, He 3970 A, 
3889 A, Ca n K 3934 A, Ca u H 3968 A, Na i D2 5890 A and Na i D 1 5896 A. All Balmer 
lines with good SNR present transient absorption features (even the noisy H/c 3750 A, 
10'*' line in the Balmer series). They are also observed in many other lines (e.g. He i, Fe i 
and Fe ii), as already pointed out by Grinin et al. (2001), but a corresponding discussion 
is deferred to a future work. 

Underlying line emission is often significant in the Balmer lines, particularly in Hfi, 
Hy and H5; in fact. Ha is seen in emission (equivalent width ~ 10 A) in the simultaneous 
INT spectra. Ca n K and H also show little emission. In these cases R is underestimated 
and negative values can artificially arise when the synthetic photospheric spectrum is 
subtracted. This is especially bad for faint absorption components near the emission 
peaks since the latter appear as minima in the R plot. In order to avoid contamination by 
these emission components, we define a new "zero level", estimated by a linear fit to the 
minima observed in the R plot. Fig. 5.4 illustrates this approach. The consistency of the 
results obtained with different lines confirms its validity. Further improvements should 
be based on theoretical modelling or higher order fitting of the emission. 

5.4 Results 

Table 5.2 gives the radial velocity shift v, the velocity dispersion Av and the absorption 
strength R of each identified transient absorption, computed according to the procedure 
described in the previous section. Uncertainties affecting the values of Table 5.2 are 
difficult to quantify, especially in the case of blended components. We are confident, 
however, that the values are a good representation of the gas properties in the case of 
sharp, isolated features. The results allow us to identify trends, based on the similar and 
consistent behaviour of many different lines. 

Absorption components at different radial velocities - either redshifted or blueshifted 
- are found in the Balmer and metallic lines. Can and Nai components with the stellar 
radial velocity are also detected. Features in different lines with similar radial velocities 
are detected within each spectrum; thus, it is reasonable to assume that they form in 
approximately the same region. The absorptions with similar velocities appearing in 
different lines are called a TAC, which can be characterized by the average of the radial 
velocity of the individual lines. A total of 24 TACs can be identified in our spectra - 
17 RACs and 7 BACs. Fig. 5.5 plots the average radial velocity of each detected TAC 
as a function of the observing time. Error bars show the rms error and the number of 
individual lines used to estimate the average velocity is indicated. Fractional numbers 
come from the fact that a weight of 1/2 is assigned to those absorption line features that 
are less certain (the blended lines He and Ca n H, the weak line Hf and the Na i doublet 
affected by telluric lines), and weight 1 to the rest. Fig. 5.5 clearly shows that more than 
one TAC is present in each observed spectrum and that there is a radial velocity shift 
when comparing TACs in different spectra. The shift does not seem to be accidental; 
there appears to be a systematic temporal evolution of the TACs from one spectrum 
to the following one when their radial velocities are analysed. This is clearly evident 
in the behaviour of the TACs observed in the four spectra taken during the night J.D. 
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Figure 5.4: Underlying emission in the Ufi absorption line. Top: Normalized observed 

spectrum (solid line) and photospheric model (dotted line). The line profile is clearly 
modulated by an emission contribution at ^ +300 km/s in the line wings. Middle: Sub- 
tracted spectrum (observed - synthetic). Bottom: R curve. The vertical mark points to an 
absorption component which would be undetectable if the continuum level were not re- 
defined. This has been done by means of a linear fit to the peaks of the emission, shown 
as throughs in the R plot. This fit (solid line at the bottom) is the new zero level for the 
gaussians. 3 gaussians (dashed lines) whose sum (dotted line) provides the optimal fit to 
the R curve are also plotted at the bottom. 
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Table 5.2: Identified ti^ansient absorption components in the Balmer, Can and Nai lines of UX Ori. 
JD is the Julian Date (-2450000). Numbers in column 3 represent the event assigned to the particular 
absorptions (see text section 4). A "0" in that column corresponds to the narrow IS absorptions, while a 
"-" means that the absorption is not associated with a particular event. Colunms 4 to 6 give: v, the radial 
velocity of the transient absorption, Av, the FWHM and R, the normalized residual absorption. 
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2451210.5 with a time interval of around 4 hours from the first to the last one. Thus, the 
data reveal groups of TACs representing a dynamical evolution of the gas with which 
they are associated. We call each of these groups an "event". In total, 7 events are 
identified, of which 4 are redshifted and 3 blueshifted. The detailed temporal evolution 
of V, Av and R for each individual line absorption component of 2 of these events (1 
blueshifted and 1 redshifted) is shown in Fig.5.6. 

5.4.1 Kinematics 

In this section we describe the kinematic behaviour of the circumstellar gas as revealed 
by our observations. Since they cover only a limited amount of time, it is difficult to 
estimate if the trends we have observed are always present, or if they are specific to our 
observations, and a different choice of times, or a longer coverage of the star, would give 
a different picture. Only additional observations, and longer time coverage, can tell us 
what is the case. With this caveat, this database is nevertheless a good starting point for 
the discussion of the TAC phenomenon in UX Ori. 

Among the detected events, the infalling gas (RACs) shows the largest velocities. For 
instance, event #1 shows a maximum average radial velocity of about 200 km/s which is 
roughly a half of the value of the stellar escape velocity, ~410 km/s, while the maximum 
average velocity of the outflowing gas (BACs) is about 100 km/s. This observational 
result could be biased since we begin detecting outflowing gas at maximum velocities, 
which is not always the case for the infalling gas (e.g. events #4 and #7). We also note 
that blueshifted gas is detected when there are redshifted absorptions, but the opposite is 
not true; we do not know if there is a physical meaning behind this result. 

The dynamics of events #1 and #2 (RACs) and #3 and #5 (BACs) denotes a decel- 
eration of the gas. In event #4 the gas first accelerates and then decelerates; the data 
of #7 suggest an acceleration, while event #6 is only present in one spectrum. An esti- 
mate of the acceleration of the events over the time interval covered by the data is given 
in Table 5.3. Positive values mean acceleration while negative ones mean deceleration. 
Table 5.3 also gives the duration of each event, which is defined as the time it takes to 
reach zero velocity (for decelerating gas) or to go from zero to the maximum observed 
velocity for accelerating gas. The estimates suggest that the events typically last for a 
few days and that accelerations are a fraction of m s"^. It is interesting to point out that 
the acceleration and deceleration phases of event #4 are of the same order of magnitude. 
We also note that our data do not allow us to conclude that this event is unique, since we 
may have missed the accelerating part of the other events. In this respect, it would have 
been of great interest to follow the evolution of event #7, which represents infalling gas 
with a significant acceleration, departing from a velocity very close to the stellar velocity 
(though it is only well detected in Hfi and Hy). 

All the transient line absorption features are very broad (Table 5.2), as is evident in 
the raw spectra. We have computed for each TAC of each event the average value of Av, 
weighted as in the case of v (the rms dispersion of Av among the different lines in a TAC 
is typically ~25%). There is a tendency for Nai lines to be narrower than hydrogen and 
Can lines (AvNai - 0.5 AvH,caii)- This will not be investigated further in this paper, but 
it could give interesting clues to the physical conditions of the gas. Fig. 5.7 shows Av 
versus v and Av/v versus v for each detected TAC, which are two complementary views. 
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Figure 5.5: Observed transient absorption components in each of the UES spectra (characterized 
by the Julian date of the observations). Each point corresponds to one "TAC" and represents the 
mean velocity of the features with similar radial velocities detected in different lines. The radial 
velocity behaviour allows us to follow the temporal evolution of the gas causing the "TAC". 7 
different dynamic events are identified; each event is enclosed in a box with an identification 
number (#1, #2,....). Error bars show the rms error of the average velocity and the number 
of lines used to estimate the average velocity is indicated. Fractional numbers come from the 
weight attributed to individual lines (see text). For the sake of clarity, the spectra taken during 
January 3P' 1999 (JD 2451210.4) are expanded along the x-axis (This figure is available in color 
in electronic form). 
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time evolution of v (km/s), Av (km/s) and R (This figure is available in color in electronic 
form). 
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Table 5.3: Estimates of the acceleration and duration of each transient event. The time 
scale for the event's lifetime is estimated as t = |vniax/«l- The acceleration and deceler- 
ation phases of event #4 are shown separately, r for event #7 means the time needed to 
achieve Vmax departing from zero velocity. 



Event 


Type 


a(ms ^) 


T (days) 


1 


infall 


-0.8 


3.0 


2 


infall 


-0.1 


8.1 


3 


outflow 


-0.5 


2.1 


4+ 


infall 


+0.8 


2.6 


4- 


infall 


-0.8 


2.7 


5 


outflow 


-0.2 


3.2 


7 


infall 


+0.6 


0.3 



Each event with the corresponding TACs are enclosed in boxes. The diagrams show that 
events are well separated and that the velocity dispersion does not change drastically 
along their temporal evolution (perhaps with the exception of event #2). This behaviour 
is remarkable in event #4 which has both acceleration and deceleration phases. Thus, 
these results suggest that events could be characterized by a kind of "intrinsic" velocity 
dispersion throughout their lifetimes. 

5.4.2 R parameter 

Table 5.2 gives the parameter R for each absorption component. The ratio of the R values 
among different lines hardly varies for all 24 identified TACs. This is shown in Table 5.4, 
where we give the ratio of I^, US, H^, Ca n K and Na i D2 to Hy. Hy has been chosen 
as the reference line because it is clearly present in all events. To estimate the errors we 
have applied a sigma-clipping algorithm to reject bad points; the corresponding values 
are given in Table 5.4. The remaining lines do not have a well defined ratio to Hy; we 
think this has no physical meaning, since He and Can H are blended and the relatively 
weak Nai Dl is more aff'ected than Nai D2 by the strong telluric absorptions in this 
wavelength interval. 

The fact that the ratio of R between difl'erent lines is rather constant allows us to 
characterize each TAC by an "intensity" (< R >) in the following manner. The ratios 
Rime/Ruy of cach TAC are used to compute an "equivalent average R^y < R >". < i? > is 
estimated from UJ3, US, Can K (weight 1.0) and (weight 0.5). Fig. 5.8 shows <R> 
as a function of the TAC velocities and events are enclosed in boxes. No correlation is 
seen. However, as in the case of Av, the estimated < R > values do not seem to change 
drastically from TAC to TAC within individual events. 

5.5 Discussion 

The intensity and dynamics of the transient components are revealed by the spectra in 
such detail that a careful analysis in the context of specific models (magnetospheric ac- 
cretion models, to mention one obvious example) could be used to validate the theoret- 
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Figure 5.7: Av versus v (top) and the complementary |Av/vI versus v diagrams of all de- 
tected TACs. Events enclosed in boxes with dotted lines are infalls and these in dashed 
lines are outflows. The arrows show the dynamical evolution (acceleration/deceleration) 
of the events. Redshifted gas appears to be broader than blueshifted gas. The location of 
each event in the diagrams suggests that its velocity dispersion does not change substan- 
tially throughout the dynamical lifetime (This figure is available in color in electronic 
form). 
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Table 5.4: Ratios of the average R parameter of several lines to Hy. The ratios (column 
2) are estimated using all TACs in which the corresponding line absorption is present. A 
sigma-clipping algorithm has been applied in order to reject bad R values. The statistical 
error, the sigma-clipping value adopted and the fraction of the rejected values are given 
in columns 3 to 5. The theoretical value for the Balmer lines is shown for comparison 

(gfline/gfHy )■ 



Line 


Ratio 


Error 


clip, cr 


% Rej. 


Theor. 




1.2 


0.2 


2.5 


13 


2.67 


m 


0.8 


0.2 


2.5 





0.49 




0.6 


0.15 


2.5 





0.18 


CaiiK 


0.9 


0.2 


2.5 


5 




NaiD2 


0.14 


0.06 


2.0 


20 





0.8 



0.6 



0.4 



0.2 



1 1 1 1 1 1 1 1 1 


1 1 1 1 


1 1 




1 


r-^c^^^ #5 

- / / / 
/ / / 

- /// #4;-,;;--... 

- . ♦ . 


♦ « ♦ 

- ' ♦ - ' 


\^ 








-> 

«- 

"Hi"' 


. -♦ ' 


■ -♦ ~ 


1 1 1 1 1 1 1 1 1 


1 1 1 1 


1 1 




1 



50 100 150 200 

|v| (km/s) 
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ical scenarios. This is in itself a huge effort, well beyond the scope of this paper. There 
are, however, several points that we think are worth mentioning. Firstly, the data indicate 
that the events are the signatures of the dynamical evolution of gaseous clumps in the 
UX Ori CS disk. There do not seem to be large changes in the velocity dispersion and the 
relative line absorption strengths through the lifetime of each clump. This suggest that 
they are 'blobs' which basically preserve their geometrical and physical identity. In this 
respect, event #4 is remarkable since during its estimated life of ~6 days it shows both 
an acceleration and a deceleration. In addition, the data reveal other aspects concerning 
the nature and dynamics of the events which are indicated in the following subsections. 

5.5.1 Origin of the circumstellar gas 

Based on spectra qualitatively similar to those discussed here, Natta et al. (2000b) anal- 
ysed the chemical composition of one infalling event in UX Ori. Non-LTE models were 
used to estimate the ionization and excitation of Balmer (up to H(5), Can and Nai lines. 
Roughly solar abundances were found for the gas causing the redshifted features of that 
event. The present data suggest a similar nature for the gas causing the identified blobs. 
We detect evidence of infalling gas in the Balmer lines and in a number of metallic lines 
with similar velocities, which also display roughly constant strength ratios to the Balmer 
lines; these ratios are < 1 (see Table 5.4). Similar result are found for the outflowing gas 
(blueshifted absorptions). Thus, we can conclude that the CS gas in UX Ori is not very 
metal rich. In addition, the data suggest that the physical conditions of the outflowing 
and infalling gas are rather similar and they likely co-exist in space. 

A second point is that in both RACs and BACs the CS gas has a significant underlying 
emission, at least in the Balmer lines we have examined. This can be seen easily if we 
compare the values of the observed average ratio for the Balmer lines to the ratio of the 
opacity in the lines (see Table 5.4, columns 2 and 6). If the CS gas was just absorbing 
the photospheric flux, then the ratio of any two lines originating from the same lower 
level should be equal to the ratio of their opacities. This is clearly not the case (see 
also similar results in Natta et al. (2000b)). There are a number of reasons why this can 
happen. The first is that the observed absorption is due to a very optically thick cloud, 
whose projected size is smaller than the stellar surface. In this case, however, R should 
be the same for all the lines, roughly equal to the occulted fraction of the stellar surface, 
and this does not seem to be the case. In particular, the Nai Dl absorption is always 
weaker than the Nai D2 (RNaiD2 / RNaioi = 1.4 ± 0.4). A second possibility is that there 
is significant emission in the lines. 

R, for any given component, can be estimated from the ratio Fobs/^syn- This can 
be done following eq. 5 of Rodgers et al. (2002), which provides this ratio under the 
assumption of a stellar and a circumstellar contribution to the observed flux; the cir- 
cumstellar contribution is caused by a spherical occulting cloud and a more extended 
envelope. In our zero order approach we neglect the extended envelope, and we assume 
a small optical depth and black-body emission for the spherical occulting cloud. In this 
case, R can be written as: 



R = T 



[1 




(5.1) 
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where r is the optical depth in the line, and Tex is the black-body excitation tem- 
perature of the gas in the cloud. Rdoud is the projected radius of the CS cloud. Some 
"fiUing-in" of the absorption features (R < r) may occur if T^x ^ T^, and i?cioud ^ R* 
(e.g. Tex - 7000^ and Rdoud/R* = 1-6, which is the corotation radius for UX Ori, 
reproduce the observed ratios). 

5.5.2 Gas Dynamics 

As already pointed out, the discussion in the previous subsection together with the results 
of Natta et al. (2000b) indicates that the TACs in UX Ori do not arise from the evapora- 
tion of solid bodies, but in a substantially different scenario. The simultaneous presence 
of infalling and outflowing gas in the immediate vicinity of UX Ori is reminiscent of the 
predictions of magneto spheric accretion models, as developed by a number of authors in 
the last several years (see Hartmann (1998) for a basic account of this theory). All these 
models assume stationary conditions, and it is difficult to derive quantitative predictions 
on the time dependent behaviour of the CS lines. However, some simple points can be 
made to constrain the different models. 

Let us start with the infall motions we detect. The simplest case we can consider 
for comparison is that of a blob of gas moving along the field lines of a stellar dipole 
magnetic field, seen in absorption against the stellar photosphere. The infalling gas is 
practically contained in a region of the size of the corotation radius, which for UX Ori is 
Rco ~ 1.6 i?^. This is in agreement with the discussion in subsection 5.1. The small value 
of the corotation radius implies that the maximum poloidal velocity the blob will reach is 
of order- Vesc(l -^star/^co) ~ 220 kms"'. Since the UX Ori disk is seen almost edge-on 
(Voshchinnikov et al., 1988), the observed velocity shift will be of the same order, which 
is in agreement with what we observe. A simple calculation shows that the blob will be 
seen to accelerate from very low velocity to this maximum value in a time interval of 
about 1 day, until it hits the star. The acceleration seen in event #4 is reminiscent of such 
a behaviour, although the time scale is longer. On the contrary, the decelerating infall 
motions we see (events #1,2,4) require a modification of this simple model. A significant 
distortion of the velocity pattern can be caused by the rapid rotation of the star, which 
can create a sort of centrifugal barrier for the infalling gas. Using the expression of the 
poloidal velocity given by eq. (8.76) of Hartmann (1998) for the UX Ori parameters, we 
find a decrease of the maximum velocity (to about 180 km s^), but a similar acceleration 
pattern and time scale. Still, the interaction of the fast rotating magnetosphere with the 
accreting disk material should be considered in detail. MuzeroUe et al. (2001) have 
included rotation in a study of T Tauri Stars (TTS), where, however, it is much smaller 
than the escape velocity at the stellar surface and the change of the poloidal velocity 
is negligible. Similar calculations for cases where the stellar rotation is a significant 
fraction of the escape velocity would be very valuable. 

Outflowing gas is seen in events #3 and #5. In both cases we start detecting the out- 
flow at negative velocities of 100-80 kms \ and we see deceleration on time scale of 
2-3 days. Stellar-field-driven wind models (X-wind, see Shu et al., 2000, and refer- 
ences therein) predict that the wind leaves the disk surface close to the corotation radius 
on almost radial trajectories, accelerating rapidly to a speed of the order of 0.i,Rco ~ 350 
kms~^ This is significantly more than the maximum outflow velocity we observe, al- 
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though the time coverage of our spectra can bias this result since we always find the 
maximum velocity of the outflow in the first spectrum taken in both the October 98 and 
January 99 observing runs. Nevertheless, the deceleration we observe is not predicted 
by the models, unless coUimation of the wind into the direction of the magnetic axis 
(i.e., away from the line of sight to the observer) occurs before the maximum speed is 
achieved and takes place on time scales of few days. Winds driven by disk magnetic field 
(see Konigl & Pudritz, 2000, and references therein) tend to produce slower outflows 
with large velocity dispersions from a broad region of the disk. Goodson et al. (1997) 
discuss models where the interaction between the rotating stellar field and a disk field 
results in a time-dependent launching of a two component outflow, of which one (the 
disk wind) has some similarity to what we observe. However, also in these models it 
is not clear what could cause the observed deceleration. It is noteworthy that, once the 
acceleration phase is over, a wind can be decelerated by the effect of the stellar grav- 
ity (Mitskevich et al., 1993). In UX Ori, the deceleration caused by the central star is 
~ 0.2 - 0.5 ms~^ at distances of ~ 14 - 20 R^,, where the local escape velocity is ~ 
130 km s"' . An episode of outflow consistent with being gravitationally decelerated was 
observed in the TTS SU Aur by Petrov et al. (1996). 

The predictions of simple magneto spheric accretion models have been compared to 
the observed line variability in a number of TTS, with moderate success. The best case 
is that of SU Aur, a rather massive TTS with a rotation period of about 3 days. This 
star shows evidence of simultaneous infall and outflow motions in various lines, similar 
to those observed in UX Ori (Johns & Basri, 1995; Petrov et al., 1996; OUveira et al., 
2000). The infalling velocity is clearly modulated by the rotation period of the star, and a 
model with magnetically channelled accretion in a dipole field inclined by a few degrees 
to the rotation axis of the star may account for most of the observed infall properties. 
As in UX Ori, the outflow is very variable, with maximum velocity similar to that of 
the infalling gas and a tendency to decelerate on a time scale of days, with no obvious 
rotational modulation (Johns & Basri, 1995). Note that the timing strategy of our UX Ori 
observations is not suited to the detection of modulations on the time scale of the stellar 
rotation period, which is about 17 hr. 

UX Ori may provide a useful test-case. On one hand, it is likely that the structure of 
the stellar magnetic field is more complex than in TTS, since UX Ori (a A-type star of 2- 
3 Mq) lacks the surface convective layer that may create a simple magnetic field in TTS. 
On the other hand, one may take advantage of the fact that the average accretion rate of 
UX Ori is rather low (Tambovtseva et al., 2001). The CS gas is mostly seen in absorption 
against the continuum, when a sporadic increase of the accretion rate in the disk creates 
gas "blobs" which move along the magnetic field lines, tracing their pattern more clearly 
than in objects with higher accretion rates, where broad emission dominates the CS 
lines. We would like to mention that Pontefract et al. (2000) also suggest magnetospheric 
accretion as a promising model to explain the Ha spectropolarimetry data of the Herbig 
Ae star AB Aur. 

One final point to keep in mind in this context is that variations along the flow (both 
for infall and outflow) of the source function in the various lines may have an important 
effect on the comparison between observations and models, so that detailed calculations 
are required for validation of the models. 
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5.6 Conclusions 

The data presented in this paper allow us to analyse the spectroscopic behaviour of the 
CS gas disk around UX Ori on time scales of months, days and hours. Significant activ- 
ity in the CS disk is always present, which manifests itself in the continuous appearance 
and disappearance of absorption components detected in hydrogen and in many metal- 
lic lines. This activity is not related to substantial variations of the stellar photosphere. 
Blobs of gas experiencing infalling and outflowing motions are the likely origin of the 
transient features. Blobs undergo accelerations/decelerations of the order of tenths of 
m s"^ and last for a few days. Detectable changes in the gas dynamics occur on a time 
scale of hours, but the intrinsic velocity dispersion of the blobs appears to remain rather 
constant. No noticeable differences are seen in the properties of the infalling and out- 
flowing gas, although infalls might have larger velocity dispersion. The relative absorp- 
tion strength of the transient absorptions suggests gas abundances similar to the solar 
metallicity, ruling out the evaporation of solid bodies as the physical origin of the spec- 
troscopic features. We suggest that the data should be analysed in the context of detailed 
magnetospheric accretion models, similar to those used for T Tauri stars. 
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Abstract 

We present high resolution (/l/A/l = 49000) echelle spectra of the intermediate mass, 
pre-main sequence stars BF Ori, SV Cep, WW Wul and XY Per. The spectra cover 
the range 3800 - 5900 A and monitor the stars on time scales of months and days. All 
spectra show a large number of Balmer and metallic lines with variable blueshifted and 



69 



70 Chapter 6. Dynamics of the CS gas in BF Ori, SV Cep, WW Vul and XY Per 



redshifted absorption features superimposed to the photospheric stellar spectra. Syn- 
thetic Kurucz models are used to estimate rotational velocities, effective temperatures 
and gravities of the stars. The best photospheric models are subtracted from each ob- 
served spectrum to determine the variable absorption features due to the circumstellar 
gas; those features are characterized in terms of their velocity, v, dispersion velocity, Av, 
and residual absorption, i?max- The absorption components detected in each spectrum can 
be grouped by their similar radial velocities and are interpreted as the signature of the 
dynamical evolution of gaseous clumps with, in most cases, solar-like chemical com- 
position. This infalling and outflowing gas has similar properties to the circumstellar 
gas observed in UX Ori, emphasizing the need for detailed theoretical models, probably 
in the framework of the magnetospheric accretion scenario, to understand the complex 
environment in Herbig Ae (HAe) stars. WW Vul is unusual because, in addition to in- 
falling and outflowing gas with properties similar to those observed in the other stars, it 
shows also transient absorption features in metallic lines with no obvious counterparts 
in the hydrogen lines. This could, in principle, suggest the presence of CS gas clouds 
with enhanced metallicity around WW Vul. The existence of such a metal-rich gas com- 
ponent, however, needs to be confirmed by further observations and a more quantitative 
analysis. 

Key words. Stars: pre-main sequence - Stars: Circumstellar matter - Stars:individual: 
BF Ori, SV Cep, WW Vul, XY Per 

6.1 Introduction 

Observations reveal that the dynamics of the circumstellar (CS) gaseous disks around 
intermediate mass, young main sequence (MS) and pre-main sequence (PMS) stars is 
extremely complex. Variable absorption components detected in many lines of differ- 
ent elements and ions constitute good examples of such complexity. The kinematics 
and intensity strength of the absorption components contain relevant information on the 
physical properties of the gas. Further, their detailed characterization and analysis pro- 
vide clues and constraints on plausible formation mechanisms as well as on theoretical 
scenarios describing the structure and nature of the CS gaseous disks. 

The presence of metal-rich planetesimals in the young MS p Pic system has been 
inferred both observationally and theoretically in a long series of papers (e.g. Lagrange 
et al., 2000, and references therein). Summarizing, there are two main arguments on 
which this inference is based. Firstly, dust causing the far-IR excess and also seen in 
the Pic disk images may be second generation material continuously replenished by 
collisions of large solid bodies or slow evaporation. Secondly, transient spectral line ab- 
sorptions, usually redshifted, of different chemical species in a wide range of ionization 
states can be modelled in terms of the evaporation of km-sized bodies on star-grazing 
orbits. Star-grazing planetesimals have also been suggested to exist in the 51 Oph sys- 
tem, a star with an uncertain evolutionary status (Roberge et al., 2002; van den Ancker 
et al., 2001); also, blueshifted absorption in excited fine structure lines of C ii* at 1037 A 
and N ii* at 1085 and 1086 A have been used to infer the presence of ~1 m bodies in the 
Vega-type binary MS system cr Her (Chen & Jura, 2003). 
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Absorption features similar to those observed in p Pic have been observed towards 
many HAe stars (see Natta et al., 2000a, for the description of these stars as a subgroup 
of HAeBe stars), particularly in the UXOR-subclass (e.g. Grady et al., 1996) and, by 
analogy, they have been interpreted as indicative of the presence of large solid bodies in 
the CS disks around these PMS stars (e.g. Grady et al., 2000, and references therein). 
In principle, this interpretation is not in conflict with the accepted time scale for the 
formation of planetesimals of ^10"* yr (Beckwith et al., 2000), which suggests that plan- 
etesimals should already be present during the PMS phase of stars (~1-10 Myr). This 
explanation for the variable absorption features observed in HAe stars is, however, con- 
troversial and, in fact, Grinin et al. (1994) pointed out other alternatives, more concretely 
dissipation of dust clouds and the simultaneous infall of cool gas onto the star. Natta et al. 
(2000b) have analyzed the chemical composition of a strong redshifted event in UX Ori 
and shown it to have a solar-like composition. Instead of the planetesimal origin for the 
transient components, those authors suggested gas accretion from a CS disk. This re- 
sult is supported by Mora et al. (2002, from now on Paper I) who presented the analysis 
of a large series of high resolution optical spectra of UX Ori. Many variable absorp- 
tion events in hydrogen and metallic lines were attributed to the dynamical evolution of 
gaseous clumps with non metal-rich, roughly solar chemical compositions. In addition, 
Beust et al. (2001) have shown that the p Pic infalling planetesimal model would not pro- 
duce detectable absorptions in typical PMS HAe CS conditions. We also note that dust 
disks around HAe stars are primordial and can be explained in the context of irradiated 
PMS CS disk models (Natta et al, 2001). 

In this paper we present high resolution spectra of the HAe stars BF Ori, SV Cep, 
WW Vul and XY Per and perform an analysis similar to that carried out for UX Ori 
(Paper I). The spectra show very active and complex CS gas in these objects; many tran- 
sient absorption features in hydrogen and metallic lines are detected, indicating similar 
properties of the gas around these stars to those of UX Ori CS gas. In addition, the 
spectra of WW Vul show metallic features without obvious hydrogen counterparts; in 
this sense, this star presents a peculiar behaviour. The layout of the paper is as follows: 
Section 6.2 presents a brief description of the observations. Section 6.3 presents the re- 
sults and an analysis of the photospheric spectra and the CS contribution. Section 6.4 
presents a short discussion on the kinematics and strength of the variable features, and 
on the metallic events detected in WW Vul. Finally, Sect. 6.5 gives some concluding 
remarks. 



6.2 Observations 

High resolution echelle spectra were collected with the Utrecht Echelle Spectrograph 
(UES) at the 4.2 m WHT (La Palma observatory) during four observing runs in May, 
July and October 1998 and January 1999. 28 spectra were obtained: 4 of BF Ori, 7 
of SV Cep, 8 of WW Vul and 9 of XY Per (XY Per is a visual binary system with a 
Herbig Ae primary and a B6Ve secondary separated by 1 .4", the good seeing during the 
observations allowed us to fully separate both components, in this paper only the HAe 
primary star has been studied). The wavelength range was 3800-5900 A and the spectral 
resolution, /l/A/l, = 49 000 (6 km s~^). Wavelength calibration was performed using Th- 
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Ar arc lamp spectra. Typical errors of the wavelength calibration are ~5 times smaller 
than the spectral resolution. The observing log, exposure times and signal to noise ratio 
(SNR) values, measured at A - 4680 A, are given in Table 6.1. Further details of the 
observations and reduction procedure are given by Mora et al. (2001). For some spectra 
there are simultaneous optical photo-polarimetric and near-IR photometric observations 
(Oudmaijer et al., 2001; Eiroa et al., 2001). Table 6.1 also presents these simultaneous 
data. At the time of the observations the stars were close to their brightest state, BF Ori 
and XY Per, or at average brightness, SV Cep and WW Vul (Herbst & Shevchenko, 
1999; Eiroa et al., 2002). 

6.3 Analysis of the spectra and results 
6.3.1 The photospheric spectra 

Circumstellar absorptions with complex profiles and blended components are detected 

in hydrogen and metallic lines in all UES spectra of the four stars. The analysis of those 
spectral features requires the subtraction of the underlying photospheric spectra. Such 
subtraction is carried out following the method outlined in detail in Paper I. Briefly, 
Kurucz (1993a) model atmospheres are used to synthesize photospheric spectra. Four 
parameters are estimated: the heliocentric radial velocity (vrad), the rotation veloc- 
ity (vsinz), the eff'ective temperature (Tgfr) and the surface gravity (log^). It has been 
shown by Grinin et al. (2001) that the photospheric lines of BF Ori and WW Vul can be 
well reproduced using solar metallicity synthetic spectra. On the other hand, the pres- 
ence of CS components in most of the lines makes it very difficult to carry out a detailed 
abundance analysis. We have thus decided to assume that the four stars have solar metal- 
licities. The atomic line data have been obtained from the VALD database (Kupka et al., 
1999). The 'best' synthetic model, defined by the parameters listed before, is selected 
by comparing some appropriate faint photospheric absorption lines among the observed 
spectra and the synthetic ones. This is not straightforward because each star behaves dif- 
ferently, and the choice of pure photospheric lines in the spectra of such highly variable 
objects is not trivial. Therefore, slightly diff'erent, ad-hoc approaches for each star are 
needed. These approaches are discussed below and the stellar parameters giving the best 
synthetic spectra are given in Table 6.2. They are compatible within the uncertainties 
with the spectral types and rotational velocities quoted by Mora et al. (2001) and with 
the results by Grinin et al. (2001), who studied BF Ori and WW Vul. 

BF Ori: The photospheric lines are narrow and do not show a noticeable variability. 
V rad is estimated using faint photospheric lines and its rms error is low. Estimated 
errors are 250 K for Tgfr (the step in the synthetic spectrum grid created using Kuruzc's 
models), and 6 kms~^ for vsin/ (the spectral resolution); the logg values considered 
have been restricted to 3.5 and 4.0. Fig. 6.1 shows the excellent agreement between the 
synthetic model and the observed median spectrum. The extra absorption seen in the 
stronger lines is due to the circumstellar contribution. 

SV Cep: According to Finkenzeller & Jankovics (1984) the radial velocity of HAe 
stars coincides within a few km s"^ with the radial velocity of interstellar (IS) lines. The 
strong Na I D IS components in the line of sight of S V Cep allow us to determine a pre- 
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Table 6.1: EXPORT UES/WHT observing logs of BF Ori, SV Cep, WW Vul and XY Per. 
The Julian date (-2450000) of each spectrum is given in Col. 1. Column 2 shows the 

o 

exposure time in seconds. Column 3 gives the SNR at A - 4680 A. Columns 4 to 6 give 
simultaneous V, %Pv and K photopolarimetric data, where available. Typical errors are 
0. 10 in V, 0.05 in ^ and 0.1% in %Pv. 



BF Ori 


Julian date 


texp (s) 


SNR 


V 


%Py 


K 


1112.6324 


1800 


280 


9.65 


0.56 


7.91 


1113.6515 


2700 


190 


9.79 


0.75 


7.85 


1209.5542 


2700 


210 






7.76 


1210.4571 


2700 


280 




0.14 





SVCep 



Julian date 


texp (S) 


SNR 


V 


%Pv K 


950.6668 


1800 


100 






950.6893 


1800 


120 






1025.6260 


2700 


140 




0.96 - 


1025.6595 


1800 


130 






1026.6684 


2700 


140 






1113.4730 


2700 


70 


11.01 


1.05 - 



1209.3372 2700 170 - 



WW Vul 

Julian date texp (s) SNR V %Pv K 

950.6176 1800 130 - = = 

950.6413 1800 130 10.89 0.69 - 

951.6232 1800 120 - - - 

951.6465 1800 150 - - - 

1023.5186 1800 140 - - 7.37 

1023.5423 2700 190 - 0.40 - 

1112.3689 1800 110 10.77 0.37 7.44 

1113.3958 2700 120 11.03 0.65 7.50 



XYPer 



Julian date 


texp (s) 


SNR 


V 


%Py 


K 


1024.6728 


1800 


200 








1024.6967 


1800 


170 


9.04 


1.49 




1025.6948 


1800 


230 




1.55 




1025.7171 


1800 


270 








1026.7065 


1800 


220 








1112.4978 


1800 


270 


9.12 


1.65 


5.97 


1113.5483 


2700 


230 


9.05 


1.53 


5.99 


1207.3786 


2700 


190 








1209.3844 


1800 


260 


9.51 


1.58 


6.18 
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Table 6.2: Stellar parameters defining the 'best' synthetic Kurucz models for each star. 
See Sect. 6.3. 1 for a discussion on the uncertainties of Teff, logg and v sin /. 



Star 


i'eff 

(K) 


logg 


V sm I 
(kms" 


vrad 
(kms-^) 


BF Ori 


8750 


3.5 


37 


23.1 + 1.9 


SV Cep 


10000 


4.0 


225 


-11.9 + 0.8 


WW Vul 


9000 


4.0 


210 


-10.4+ 1.2 


XYPer 


8500 


3.5 


200 


8.3 + 0.6 



1.2 



.-^ 1 

t/o 

CD 

S 0.8 



0) 



0.6 - 



0.4 



1 I I I I r 




Observed spectrum 
Synthetic spectrum 

I I I I I I I L 



J_ 
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Figure 6.1: Synthetic (dashed line) vs observed spectrum (solid line) for BF Ori. The 
agreement between the photospheric and the observed median spectrum is remarkable 
over a large wavelength range. The extra absorption in the stronger lines (identified in 
the figure) is due to the circumstellar contribution. This example illustrates the need to 
use very faint lines to estimate the stellar parameters. (This figure is available in color in 
electronic form) 
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Wavelength (A) Observed spectra Synthetic spectra 

Figure 6.2: Synthetic vs observed spectra for SV Cep. Left: comparison of the in- 
termediate resolution Hor spectrum and two synthetic Hor profiles. The best result is 
obtained for Teff = 10000 K and log^ = 4.0 (blue dashed line). It was assumed that 
vsinz = 206 kms"^ (Mora et al., 2001). Right: Median UES spectra (solid lines) 
of several blends compared with the best synthetic spectra, with Tgff = 10000 K, 
vsin/ = 225 km s~^ (blue dashed lines). (This figure is available in color in electronic 
form) 



cise value of the radial velocity of this star. Its high temperature, rotation velocity and 
variability make it extremely difficult to identify photospheric lines in order to derive 
the Kurucz stellar parameters. We have circumvented this problem by using the EX- 
PORT intermediate resolution spectra (Mora et al., 2001), which have very high quality 
Ha profiles and are of great help in making a first estimate of Teff and log g. From 
the Ha wings and the almost complete absence of weak photospheric lines in the high 
resolution spectra, T^s - 10000 K and logg = 4.0 are estimated. It was assumed that 
vsini = 206 kms"^(Mora et al., 2001). The adequacy of this choice can be seen in 
Fig. 6.2, where the observed median SV Cep Ha line and 2 synthetic Balmer profiles 
(Teff = 10000 K, log^ = 4.0, vsin/ = 206 kms"' and Teff = 10000 K, log^ = 3.5, 
V sin / = 206 km s"') are shown. Using those values we have identified 7 blends of weak 
photospheric lines (absorption intensities < 4% of the continuum) with very low vari- 
ability. The blends are at ~ 3913, 4129, 4176, 4314, 4534, 5041 and 5056 A. A grid of 
photospheric spectra has been generated with different values of Teff and v sin i (log g was 
assumed to be 4.0 from the Ha analysis) and the rms differences between the synthetic 
and observed spectra have been estimated. The lowest rms difference is obtained for 
Teff = 10 000 K and v sin / = 225 km s"^ . The differences between the Ha wing profiles 
broadened to v sin i = 225 km s"^ and v sin / = 206 km s"' are negligible, so it was not 
needed to compute log^ again. Fig. 6.2 also shows the best synthetic and the observed 
median spectra for the selected blends. Uncertainties of ~500 K (the step of the Kurucz's 
models) and ~10% for Teff and v sin i , respectively, are estimated. 

WW Vul: V rad is estimated from the sharp Nai D IS lines. The star was very active 
during the observing runs and its spectra show a large number of broad and variable 
absorption features, e.g. redshifted features are superimposed on practically each photo- 
spheric line in the JD 1 1 13.40 spectrum, which pose severe difficulties for the selection 
of appropriate photospheric lines. Nevertheless v sin i has been estimated with the Mg ii 
4481 A blended doublet. Many pairs of Teff - \ogg values can reproduce the observed 
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Figure 6.3: The observed median UES spectrum of WW Vul (black continuous line) 
compared to the broadened synthetic one (blue dashed lines) in two different spectral 
regions. At the top of the figure, the unbroadened Kurucz model is shown with some 
line identifications (red solid line). (This figure is available in color in electronic form) 

spectra, though the agreement between the observed and synthetic spectra is generally 
poor. We have selected Ten = 9000 K, log g = 4.0 as representative of the WW Vul pho- 
tosphere because of its compatibility with the results by Mora et al. (2001), though large 
errors are likely, but not larger than around 10% in both Tgff and v sin /. Fig. 6.3 shows a 
comparison between the observed median spectrum and the broadened synthetic one in 
two spectral regions where a good fit is achieved; this figure also shows the unbroadened 
Kurucz's model with some line identifications. 

XY Per: vrad is estimated from the Nai D IS components. The photospheric lines 
are very broad but the high SNR of the spectra and the relatively low Tgff allow us to 
identify a large number of faint line blends in order to perform a precise estimate of Tgff, 
log g and V sin z (16 faint blends with absorption intensities < 8% of the continuum with 
very little CS activity could be identified). Errors of Tgff and log g are of the order of the 
step in the Kurucz models, 250 K and 0.5 respectively, while the error in v sin i is very 
low, <10%. The comparison between the synthetic and the observed median spectra of 
the 16 blends is shown in Fig. 6.4. 

6.3.2 The circumstellar transient absorption contribution 

Once the best photospheric spectrum of each star is determined, the circumstellar con- 
tribution to each observed spectrum can be estimated by subtracting the synthetic one. 
The residual spectra show transient absorption features, which can be characterized by 
means of the normalized residual absorption, defined as R = 1 - Fobs/^syn (Natta et al., 
2000b). The R profile of each line reflects the blending of several components. A multi- 
gaussian fit providing the radial velocity, velocity dispersion and absorption strength is 
used to identify the individual components (see Paper I, for details). Broad redshifted 
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Figure 6.4: Synthetic (dashed lines) vs observed median (continuous lines) spectra of 
XY Per for 16 spectral features with very low CS contribution. The stellar parameters 
of the synthetic spectrum are given in Table 6.2. (This figure is available in color in 
electronic form) 



and blueshifted absorptions at different radial velocities are found in the Balmer and 
metallic lines for all 4 stars analyzed. We apply the multigaussian fit to Balmer lines 
(ti/3 4861 A, Hy 4340 A, Ud 4102 A, He 3970 A, 3889 A), Can K 3934 A, Can H 
3968 A, Nai D2 5890 A and Nai Dl 5896 A, as well as to fainter metallic lines Fe ii 42 
multiplet (a6S-z6Po triplet: 4924 A, 5018 A and 5169 A), Tin 4444 A, Tin 4572 A, 
Fei 4046 A, Sen 4247 A and Cai 4227 A. We have chosen these ionic lines because 
they show significant CS variability and are relatively strong and isolated. Narrow IS 
components (mainly Na i, Ca ii and Fe ii) with the stellar radial velocity are also detected 
for all the stars. 

Consecutive spectra with a time delay of ~1 hour of SV Cep, XY Per and WW Vul 
were taken on several nights (see Table 1). These spectra were quite similar and the gaus- 
sian deconvolution of the R profiles essentially provides the same values for the fit pa- 
rameters; thus, any significant variation of the phenomena causing the transient absorp- 
tions is excluded on this time scale, at least during these observing periods. This result 
gives us confidence in the identification of the components and allows us to add the spec- 
tra taken during the same night in order to increase the SNR. Tables 6.6, 6.7, 6.8 and 6.9 
give the radial velocity shift v, the velocity dispersion Av and the absorption strength 
^max, the peak of the R profile, of each identified broad transient absorption component 
of the lines listed above for BF Ori, SV Cep, WW Vul and XY Per, respectively. Col- 
umn 1 gives the the corresponding Balmer or metallic line. Col. 2 gives the Julian Date, 
Col. 3 represents the event assigned to the particular absorptions. Col. 4 gives the radial 
velocity shift v. Col. 5 lists the velocity dispersion Av and Col. 6 gives the absorption 
strength Rmax- JD values in Tables 6.7, 6.8 and 6.9 correspond to the starting time of 
the first spectrum of each night. In the following, whenever a JD is given, 2 450 000 is 
subtracted. 

Tables 6.6, 6.7, 6.8 and 6.9 show that absorption components with similar radial ve- 
locities appear/disappear simultaneously in different lines, as is observed in UX Ori 
(Paper I). We assume that absorptions with similar velocities come from the same 
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gaseous clump, which can be characterized by an average radial velocity <v> (a Tran- 
sient Absorption Component or TAC). The time evolution of the TACs' velocity is re- 
ferred to as an event and represents the dynamical evolution of the gaseous clumps. We 
point out that there is an uncertainty in identifying TACs detected on different nights 
with the same gaseous clump; our assumption relies on the UX Ori results (Paper I, 
see also below). Figs. 6.5, 6.7, 6.9 and 6.11 plot <v> of the identified TAC versus JD. 
<v> is a weighted average in which the lines H/3, Hy, US, Can K, Feu 4924 A and 
Fen 5018 A have an arbitrarily assigned weight of 1 because of their higher intensity 
and non-blended nature. A weight of 1/2 is also arbitrarily assigned to the remaining 
lines in order to reflect that they are either blended (He, Can H and Fen 5 169 A), weak 
(Sen 4247 A, Fei 4046A and Cai 4227 A) or affected by telluric lines (Nai D2 & Dl). 
The weighted number of lines used in each average, which can be fractional because 
of the 1/2 weights, is plotted next to each point. The TACs have been grouped ac- 
cording to the event they represent; thus, the figures show the dynamical evolution of 
the gaseous clumps. Note that some events were only detected once (only one TAC). 
Figs. 6.6, 6.8, 6.10 and 6.12 show the R profiles of some selected lines for the four stars. 
The line absorption components and the corresponding event identifications are indi- 
cated. The lines are Ufi, Can K (except for BF Ori), Nai D2, Nai Dl and Feu 5018 A. 
Hy is shown for WW Vul (Fig. 6.10) since Hfi has a large underlying emission. The 
results for each star are presented in the following. 

BF Ori: 5 TACs grouped in 3 different events were detected in the spectra of BF Ori 
(Fig. 6.5). Event #1 is an accelerating redshifted event, first detected in the JD 1 1 12.63 
spectrum at approximately the stellar radial velocity. Absorbing gas is seen in Balmer 
and metallic lines of Na i, Fe n, Ti u, Sc u, and Ca i. The strongest H i lines (^S, Hy and 
US) appear to be saturated, i.e. i?max is very close to unity, and the metallic lines are also 
very strong. The Balmer lines are broader than the metallic ones. The parameters of at 
least some of these lines might be influenced by IS gas absorption (a careful look at the 
Nai D lines shows the presence of two peaks in the JD 1 1 13.65 spectrum, Fig. 6.6). The 
velocity dispersion of the lines tends to be larger when the event increases its velocity, 
while i?max values tend to decrease, though changes are modest. Event #2 is a strong 
redshifted decelerating event; the behaviour of its TACs is in general similar to those 
of #1. #3 represents blueshifted gas only detected in the last night of January 99 and 
is fainter than the redshifted ones. The R^^^ values of the Balmer lines are low, but 
they might be saturated, because i?max does not appear to decrease as we follow the 
series in what will be called the "expected Balmer decreasing trend" for optically thin 
gas (/?max(I^) > ^max(Hy) > /?max(H5) > i?max(H^)); also the relative intensity of the 
metallic lines with respect to the H i ones seems to be lower. 

The Can K line has non-photo spheric profiles with the simultaneous presence of 
redshifted and blueshifted components, but their radial velocities do not match the ab- 
sorptions observed in other lines (except the blueshifted TAC in JD 1210.45). 

SV Cep: Broad absorptions of H i, Ca n, Na i and Fe n are detected in the spectra of 
SV Cep, but no variability is found in Ti n, Sc u, Fe i and Ca i (unlike the other stars in the 
paper). The broad absorptions represent 10 TACs grouped in 8 different events: 3 of them 
correspond to outflowing gas and the remaining 5 to infalling gas (Figs. 6.7 and 6.8). 
Blueshifted gas shows small radial velocities, on average < 20 kms \ while redshifted 
components display velocities as high as 160 kms~^ Only one spectrum was taken in 
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Figure 6.5: Events in BF Ori. Each point corresponds to the radial velocity of one TAC 
and represents the average velocity <v> of the absorptions with similar radial velocities 
detected in different lines. <v> is a weighted average in which the higher intensity non- 
blended lines are assigned a weight of 1 and the rest 1/2. Error bars show the rms error 
of the average velocity; the numbers above the data points indicate the weighted number 
of lines used to estimate the average. Fractional numbers arise from the 1/2 weight 
attributed to some lines (see text). Redshifted events (infalling gas) are printed in red 
colour, while blueshifted events (outflowing gas) are in blue. (This figure is available in 
color in electronic form) 



80 Chapter 6. Dynamics of the CS gas in BF Ori, SV Cep, WW Vul and XY Per 




Figure 6.6: BF Ori R profiles. The normalized residual absorption profiles (R = I - Fobs 
/ Fsyn) of Ufi (left), Nai D2 and Dl (middle) and Feii 5018 A (right) are shown in the 
figure (black colour). The corresponding spectra are indicated in the vertical right axis 
(Julian and civilian epochs). The identified TACs (gaussian components, pink colour) 
and the reconstructed R profile fit (blue colour) are displayed. Event numbers are shown 
under the gaussians. The zero velocity interstellar components are marked as "IS". (This 
figure is available in color in electronic form) 
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Figure 6.7: Events in SV Cep. Details as for Fig. 6.5 (This figure is available in color in 
electronic form) 



May 98, October 98 and January 99, i.e. the events of these periods are composed of 
1 TAG only: these data represent isolated snapshots of the CS gas around SV Cep and 
no temporal evolution can be inferred from them. The 4 TACs detected in July 98 can 
be grouped in two events: #3 corresponds to redshifted gas with practically constant 
radial velocity, and #4 is gas observed at a velocity close to the stellar radial velocity 
or slightly blueshifted. In general, the Rmax values of the Balmer lines show the ex- 
pected Balmer decreasing trend and are broader and much stronger than the metallic 
lines, 7?max(H5) /i?niax(Fe II) > 3. There are, however, some exceptions. #1 shows rela- 
tively strong Fe ii and it is not clear that the Balmer lines show the expected decreasing 
trend. In #5 the Hi lines are not broader than the Feu ones, and the Balmer lines are 
probably saturated (but note that this event is very faint, and it could be an overinterpre- 
tation of the fit procedure). The strongest Hi lines are saturated in #7 and also in the 
JD 1026.66 spectrum of #3. There are anticorrelated changes in the Av and Ry^ax values 
of #3, but they show the same trend in #4. 

WW Vul: 15 TACs grouped in 9 different events are identified in the 5 spectra of 
WW Vul (Figs. 6.9 and 6.10). 4 events are seen in May 98. #1 is redshifted gas with 
saturated Balmer lines and strong Feu lines. Both Hi and Feu lines have similar Av 
and from one TAC to another i?max and Av show opposite trends. #2 corresponds to low 
velocity blueshifted gas clearly detected in the metallic lines but no counterpart in the hy- 
drogen lines is apparent (see the Fe ii 5018 A line in Fig. 6.10). #3 and #4 are blueshifted 
accelerating events. The only spectrum of July 98 (JD 1023.51) reveals 3 TACs: #5 is 
a broad, Av > 100 kms \ redshifted event only detected in metallic lines (the broad 
wing of the Feu 5018 A line profile in Fig. 6.10). #6 is a very low velocity blueshifted 
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Figure 6.8: SV Cep R profiles. Details as for Fig 6.6. (This figure is available in color in 
electronic form) 

component (v - -5 km s"^); this event is significantly broader in the H i lines (which are 
saturared) than in the metallic ones (the IS contribution cannot be separated from this 
low velocity event). #7 is a relatively narrow blueshifted event only detected in metallic 
lines and is clearly distinguished as a peak in the line profiles. Again, metallic red- 
shifted absorptions at v - 90 km s"' without H i counterpart are detected on JD 11 12.36 
(October 98). Similar metallic absorptions are also detected on JD 11 13.39, but on this 
date saturated hydrogen components with basically the same kinematic parameters (in- 
cluding the velocity dispersion) are present. We tentatively identify both TACs with the 
same event, #8, although we cannot exclude the possibility that the metallic absorptions 
detected on each date could be due to diff"erent gas. Finally, #9 is a strong, redshifted, 
decelerating event, identified in H i and in many metallic lines. In this case, the hydrogen 
lines are considerably broader than the metallic ones. The Hi lines seem saturated on the 
first night, while on the second one the expected Balmer decreasing trend is observable 
and on both nights the Fe ii lines are relatively strong. 

XY Per: 16 TACs grouped in 9 difl'erent events are detected in the 7 XY Per spec- 
tra. #1 is a faint redshifted event detected on the last two (out of three) July 98 nights. 
#2 corresponds to blueshifted accelerating gas; all three TACs of this event are strong 
and broad. #3 corresponds to a relatively narrow and faint blueshifted event only de- 
tected on JD 1025.69. #4 is also a narrow and faint event, centered at approximately the 
stellar radial velocity, only detected on JD 1026.71; IS absorption could contribute to 
the metallic Ry^ax values. #5 represents decelerating redshifted gas. #6 is a strong red- 
shifted event only detected on JD 1112.49. #7 corresponds to decelerating blueshifted 
gas; the expected Balmer decreasing trend is clearly present in both TACs and metallic 
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Figure 6.9: Events in WW Vul. Events marked with a "#" (1, 3, 4, 6, 8 and 9; red 
for redshifted events and blue for blueshifted events) are detected both in hydrogen and 
metallic lines. Events denoted with an "@" (2, 5 and 7; green) are only seen in metallic 
lines. The square point of #8 in JD 1112.37 correspond to TACs only observed in Feii 
and Till. Further details as for Fig 6.5. (This figure is available in color in electronic 
form). 
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Figure 6.10: WW Vul R profiles. Details as for Fig 6.6. (This figure is available in color 
in electronic form). 

lines are very strong. Both Av and R^ax increase very significantly from JD 1 1 12.49 to 
JD 1 1 13.55. #8 is decelerating redshifted gas detected on January 98 (note that the time 
interval between the two TAGS grouped in this event is 48 hours); in both TACs the 
metallic lines are very broad and strong: Feu lines are even broader and stronger than 
the Balmer ones which are faint and seem saturated. Finally, #9 represents accelerating 
blueshifted gas. 

6.4 Discussion 

The objects studied in this work and UX Ori are very much alike. All are PMS HAe stars. 
UX Ori, BF Ori and WW Vul are bonafide UXOR-type objects (e.g. Grinin, 2000), i.e. 
their light curves show high- amplitude variability (Am > 2.0 mag), Algol-like minima, a 
blueing efl'ect and an increase of the polarization when the object brightness decreases. 
SV Cep also shows UXOR characteristics (Rostopchina et al., 2000; Oudmaijer et al., 
2001). Further data are required before XY Per can be confidently classed as a UXOR- 
type object (Oudmaijer et al., 2001), although it does share the same complex and vari- 
able spectroscopic behaviour of the other stars. 

The time covered by the UX Ori spectra allowed us to analyze and identify the TACs 
as due to the dynamical evolution of gaseous moving clumps; this identification and 
dynamical evolution was particularly convincing in the case of a clump detected in four 
spectra taken within a time interval of four hours (Paper I). We rely on those results 
to ascribe to the same gaseous clump groups of H i and metallic absorption components 
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Figure 6.1 1: Events in XY Per. Details as for Fig 6.5. (This figure is available in color 
in electronic form). 



described in the previous Sect. We are aware, however, that this identification is doubtful 
in some cases, since the number of spectra for the stars in this paper is smaller and the 
time coverage is poorer than for UX Ori. With this caveat in mind, we will discuss these 
spectra as we have done in Paper I for UX Ori. We are quite confident that the main 
conclusions of this paper are not affected by the uncertainties with which some specific 
event can be identified. 



6.4.1 Kinematics 

Accelerating/decelerating blueshifted and redshifted events were detected in UX Ori, the 
events seemed to last for a few days and their acceleration rates were a fraction of a m s"^. 
Within the present time coverage limitations, the same is observed in the outflowing and 
infalling gas of BF Ori, WW Vul, SV Cep and XY Per. All 5 stars share the trend 
that infalling gas shows the largest velocities (the exception is event #4 of WW Vul in 
the JD 951.62 spectrum) and that blueshifted absorptions are detected when redshifted 
absorptions are present (the exception might be the JD 1024.67 spectrum of XY Per 
which shows a low velocity outflow but no infalling gas). Similar results are also present 
in the spectra of Grinin et al. (2001). 

Infalling gas appears to have larger dispersion velocities than outflowing gas. In 
addition, the H i lines are broader than the metallic ones in approximately 40 % of the 
identified TACs, as indicated, for example, by the fact that Avpen < 0.66Avhi (in UX Ori 
this trend was noticed in the Nai lines). The fraction of TACs with broader Hi lines 
varies from star to star. For example, in 12 out of 16 detected TACs in XY Per Hi and 



86 Chapter 6. Dynamics of the CS gas in BF Ori, SV Cep, WW Vul and XY Per 



H/S 



Call K 



Nal DS & Dl 



Fell 5018A 




-SOO -260 ZiO SOO -500 -2S0 2S0 500 -500 -250 250 500 -500 -250 260 500 

Radial velocity (km/s) 



Figure 6.12: XY Per R profiles. Details as for Fig 6.6. (This figure is available in color 
in electronic form). 
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Fell have similar velocity dispersion, and only 3 out of 10 TACs in SV Sep follow this 
trend. On the other hand, the Fe ii lines are broader than the H i lines in only two events, 
#5 of SV Cep and #8 of XY Per. It is worth noting that the event of XY Per has metallic 
^max values at least as large as those of the B aimer lines and that the SV Cep event is 
very weak (see above for this event). 

There seems to be a correlation between the dispersion velocity and the velocity of 
the TACs in the sense that TACs appear to be broader when the velocity increases. If 
we perform a linear regression (Av = A + B x |v|) to the whole set of data, we find a 
correlation coefficient of 0.68. There is also a suggestion of an anticorrelation between 
<AvFeii> and the < Rmax > values of the TACs during their evolution, i.e. events become 
fainter when they increase their dispersion velocity. This was also suggested in the 
case of UX Ori. For every event with more than 1 observation (14 in total), <AvFen> 
and < i?niax > have been normalized to the first observed values. If 2 'anomalous' 
events, which are the most uncertain identifications in the whole sample, are removed 
(WW Vul #8 - one TAC is observed only in the metallic lines while the second TAC 
is detected in both metallic and Hi lines -, and XY Per #7 - the Av and < R^aax > 
variations are much more extreme than those of all other detected events) we find a 
linear correlation coefficient of -0.66. Both correlations are highly significant, i.e., the 
probability of randomly obtaining such coefficients from 2 unrelated variables is <1% for 
the data sets considered. We have to point out, however, that the <AvFeii> vs < i?max > 
anticorrelation dissapears if the two 'anomalous' events are included in the statistics. 

6.4.2 Line intensity ratios 

Many absorption components in Tables 6.6 to 6.9 have i?max values close to unity, which 
suggests that they are saturated. This is not the case even for the strongest events of 
UX Ori (Paper I) (in fact, the ratio among the R^^x line values of the 24 UX Ori TACs 
does not vary much, which allowed us to estimate a line residual absorption average). 
Nevertheless, we have followed the procedure of Paper I to investigate whether a "fixed" 
ratio among the line absorption strengths might be present in those TACs which are most 
likely unsaturated. Thus, we have excluded from this exercise lines wich show signs 
of being saturated, i.e. those with Rj^a > 0.8 and TACs with Balmer lines of similar 
strength. The line Fe ii 4924 A has been taken as a reference as it shows the lowest sta- 
tistical errors (other lines which have been considered are H^, Ca ii K and Fe ii 5018 A). 

The ratio i?max, line /^max,ref. line (e.g. ^max,Hd/^niax,Fen4924A) haS been COmpUtcd for each 

line of every TAC, and later the mean <i?max,une/^max,ref .iine> has been estimated using a 
sigma-clipping algorithm to reject bad points. Table 6.3 gives <i?max,une/^max,ref.iine> for 
each star, together with statistical errors and line rejections (%). 

Values quoted in Table 6.3 can be used to assess whether the gas causing the variable 
absorptions is optically thin or not, simply by comparing the i^max ratios of lines belong- 
ing to the same element multiple! with their gf ratios (we recall that the TACs suspected 
of being saturated have already been excluded). We have analysed the Balmer, the Na i D 
doublet and the Fe ii 42 lines, for which US, Na i D2 and Fe ii 501 8 A have been taken as 
reference, respectively. i?max ratios have been computed following the above procedure 
and gf values have been taken from the VALD database (Kupka et al., 1999) for Hi and 
Nai. Fen lines do not have reliable experimental gf values, partly because Fen 5169 A 
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Table 6.3: Ratios of the average i?niax parameter of several lines to Fen 4924 A 
(< ^max,ime/^max,Feii4924A >) f^r each Star. The values correspond to lines in TACs which 
are most likely not saturated (see text). Values with no error mean that only one TAC is 
available. The sigma-clipping threshold adopted is 2.0 cr, except for the values followed 
by the symbol f, in which 1.5 cr has been used. The percentage of rejected lines is given 



in brackets. 


Line 


BFOri 
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5.23 ± 1.82 
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3.40 ± 0.83 (14%) 
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CaiiK 


5.79 
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2.37 ± 0.65 


2.20 ± 1.16 (8%) 


NaiD2 


1.06 + 0.40 


1.27 ± 0.40 


1.45 ±0.64 (17%)f 


0.76 ±0.12 (20%) 


NaiDl 


0.89 ± 0.39 


0.85 ± 0.33 


1.50 ±0.92 


0.48 ±0.09(11%) 


Fen5018A 


1.08 ±0.15 


1.41 ±0.14 


1.28 ±0.11 (15%) 


1.16 ±0.17 (14%) 


Feu 5 169 A 


1.18 ±0.15 


1.52 ± 0.30 


1.37 ±0.13 (15%) 


1.30 ±0.25 (7%) 


Tiu4444A 


0.57 ± 0.29 




0.25 ± 0.08 


0.30 ±0.14 (7%) 


Tiu4572A 


0.45 ± 0.07 




0.25 ± 0.08 


0.32 ±0.10 (7%) 


Scii4247A 


0.47 ± 0.09 




0.31 ± 0.06 


0.33 ± 0.23 


Cai4227A 


0.32 ± 0.06 




0.18 ±0.11 


0.26 ± 0.08 


Fei4046A 






0.16 ±0.06 


0.18 ±0.05 



o o o 

is blended with Mgi 5167 A, Fei 5167 A and Mgi 5173 A. Following the suggestion 
made by T.A. Ryabchikova and F.G. Kupka (private communication), we have used the 
semiempirical values computed according to Raassen & Uylings (1998) and available at 
http://www.science.uva.nl/research/atom/levels/levtext.html. The comparison between the 
estimated ratios and the theoretical gf shows that the Fe n 42 triplet is most likely satu- 
rated in BF Ori but not in the rest of the stars, and that the Na i D doublet is also probably 
saturated in BF Ori and WW Vul (Table 6.4). Concerning the Balmer lines (note we are 
referring to those TACs apparently unsaturated) their ratios are very different from the 
gf ratios. This is, in principle, similar to the case of UX Ori where the lines do not 
seem to be saturated in any of its events. Paper I suggested that the UX Ori results could 
be explained by underlying line emission caused by a spherical occulting cloud with a 
temperature Tgx ~ 7000 K and a radius of the order of the UX Ori corotation radius, 
RcIovlaIR* ~ 1-6 (see Rodgers et al., 2002, for details of the assumptions). However, only 
XY Per presents Balmer line ratios which could be adjusted using this scenario, namely 
gas at a Tg^ ~ 6600 K at approximately the corotation radius Rdoud/R* ~ 1-6. This is 
not applicable for the rest of the stars, where it is most likely that the Balmer lines are 
always saturated. 

6.4.3 Origin of the variable circumstellar gas clumps detected in hy- 
drogen and metallic lines 

Excluding the fact that redshifted events seem to be detected at larger velocities than 
blueshifted events (a fact deserving further observations and also a theoretical expla- 
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Table 6.4: Estimated < i?max,ime/^max,ref.iine > ratios among the lines of the Nai D dou- 
blet and the Feii 42 triplet for each star. The theoretical ratios (gfune / gfref.iine) and the 
line taken as reference are given in the last Col. The Feii 5169 A line is included for 
comparison purposes though its ratios are likely affected by the blend with Mg 1 5 167 A, 
Fei5167 Aand Mgi5173 A. 



Line 
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Reference 
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0.83 ±0.11 


0.73 ± 0.15 


0.81 ±0.12 


0.63 ± 0.07 


0.50 


NaiD2 


Fell 4924A 


0.94 + 0.16 


0.72 + 0.08 


0.79 + 0.07 


0.86 ±0.17 


0.69 


Fell 501 8A 


Feu5169A 


1.09 ±0.03 


1.07 ±0.15 


1.06 ±0.13 


1.07 ±0.19 


1.25 


Fell 501 8A 



nation), the gas is observed in the same absorption lines with similar R and velocity 
dispersions, i.e., there are no fundamental differences in the behaviour of the infalling 
and outflowing gas in all the observed stars, including UX Ori. This result suggests that 
their physical conditions are rather similar and that they probably originate at similar 
distances from the star. Roughly solar abundances were found by Natta et al. (2000b) 
for an event of redshifted gas in UX Ori, which was further supported by Paper I. This 
led to the conclusion that the clumps of CS gas in UX Ori are non metal-rich. Since 
high velocity gas is observed simultaneously in the Balmer and metallic lines in BF Ori, 
SV Sep, XY Per and in most of the detected events in WW Vul, a similar conclusion very 
likely holds also for the CS gas in these stars (but see next Sect, for the metallic events in 
WW Vul). Unlike the detected events in UX Ori, the gas in those stars is often optically 
thick, as suggested by the saturation of the absorption features. This might indicate that 
high density gas is more frequently observed in BF Ori, SV Sep, WW Vul and XY Per 
than in UX Ori. 

Paper I compares the dynamics of the gaseous clumps in UX Ori with the predictions 
of magneto spheric accretion models (see Hartmann, 1998, for a very good basic descrip- 
tion of this theory) and different wind models (e.g. Goodson et al., 1997; Shu et al., 2000; 
Konigl & Pudritz, 2000). The present data do not add any new substantial aspect to that 
discussion, only that, with very few exceptions, outflowing gas displays smaller veloci- 
ties than infalling gas. Thus, to avoid repetition we refer to that paper, stressing the need 
for further theoretical efforts to explain the complex circumstellar environment of HAe 
stars, at least to the level of understanding achieved for the less massive T Tauri stars. 

6.4.4 The intriguing case of WW Vul 

WW Vul is a very interesting case, somewhat different from the other stars. In all the 
spectra we have obtained, in addition to events detected in both metallic and hydrogen 
lines, as it is always the case in the other objects, we see also metallic absorption features 
(both blueshifted and redshifted) that do not seem to have a counterpart in the hydrogen 
lines. They appear as broad high velocity wings in the R profiles, e.g. events #5 and 
#8 in the JD 1023.52 and 1112.37 spectra, as well as relatively narrow distinct peaks, 
e.g. #7 in JD 1023.52 (Fig. 6.10). The observed outflowing events are narrower and 
have larger 7?n,ax values than the infalling ones. To ensure that the metallic events are not 
artificially introduced by the multigaussian fit we have proceeded in two different ways: 
1. We have carefully inspected the original spectra, i.e., before subtraction of the photo- 
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spheric components. Fig. 6.13 shows several lines as observed in the original (without 
photospheric subtraction) JD 1023.52 WW Vul spectrum. In particular, the relatively 
narrow, blueshifted, metallic event #7 is clearly seen on a simple visual inspection as an 
absorption peak with no distinctly separate hydrogen counterpart. We do note that there 
is hydrogen absorption covering the velocity range of event #7 and that the hydrogen 
lines are quite broad and likely to be saturated at this time. As a result, Gaussian fitting 
may not properly represent the intrinsic shape of the hydrogen components; however, 
without advance knowlegde of their intrinsic shape, Gaussian fitting is used here as a 
default procedure. 2. In general, our approach has been to use a number as small as 
possible of gaussians to reproduce the R profiles in a self -consistent way, avoiding to 
overfit the data. 

Nevertheless, since we were intrigued by the WW Vul behaviour, we have performed 
a number of additional numerical tests, to check as well as possible that the appar- 
ent lack of a hydrogen counterpart to some event seen in the metallic lines was not 
an artifact. First of all, we have added to the Hi identified components more gaus- 
sians which take into account the parameters of the absorptions only identifed in the 
metallic lines. This led to unphysical results, i.e., adding more gaussians to the Hi 
lines does not produce kinematic components similar to those observed in the metal- 
lic lines, but create spurious Hi components. For example, in JD 1023.52 we tried to 
fit 3 components to the H i R^xi profiles, i.e. the same number of gaussians obtained 
from the Fe n fit. The data obtained from the B aimer lines deconvolution are not self- 
consistent: i) strange, unrealistic Balmer intensity sequences are obtained, e.g. in event 
#7 ^max,H4- = -^max.H/S > -^max.Hy > -^max.Hd- ii) the intensity ratios between different events 
differ widely from line to line (e.g. H^ is about two times more intense in #5 than in 
#6, Hy is only 0.8 times as intense in #5 as in #6. We have then also tried to fix the 
radial velocity of the gaussian H i components to the values obtained from the Fe n lines 
in order to improve the fit, but the results were even worst: the method diverged for H/3 
(i.e. no fit could be found) and the sequence of intensities for the Balmer series in event 
#6 was unphysical (Rm^Mi = Rm^,us = 2.57?max,Hr, no Ufi component). 

One can also argue that the absence of Hi counterparts to the metallic events is 
a consequence of blending, which is unresolved by the multigaussian deconvolution 
procedure. In order to test such possibility, we have generated a composite synthetic R 
profile consisting of three gaussians with the same R^^ax and v as the Feu 5018 A JD 
1023.52 spectrum and Av = 150 km/s, which is the average velocity dispersion of the 
only H I component detected in that date. We have also added random noise to achieve a 
S/N ratio worst than that of the I^ and Hy lines of that spectrum; finally, we have applied 
our multigaussian deconvolution procedure to the synthetic profile. The result is that 
the three individual gaussians have been successfully retrieved. Note that the gaussian 
parameters we have introduced are test values. This result suggests that we would have 
been able of identifying H i counterparts of the metallic kinematics components if they 
would exist. Therefore, we consider that the existence of metallic absorptions without 
obvious H I counterparts is a rather firm result. 

There are no obvious differences, in terms of kinematics and absorption strength, 
between the WW Vul metallic events and the rest of the events detected in this star and 
in BF Ori, SV Cep, XY Per and UX Ori. Some remarkable differences appear, however, 
when the ion column densities causing the events are compared. Lower limits on the ion 
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Figure 6.13: Original spectrum (i.e. prior to the synthetic stellar spectrum subtraction) 
as observed in the JD 1023.52 spectrum of WW Vul. Wavelengths have been converted 
into radial velocities. Plotted lines are Hy, and the three lines of the Feu 42 triplet. 
The photospheric synthetic spectrum used to form R is shown in dotted lines for Hy and 
Fen 4924A (being very similar for the other lines). The synthetic spectra reveal the 
presence of emission in the wings of Hy and a little in the blue wing of Fell 4924A. 
The average radial velocity <v> of the events identified after the photospheric spectrum 
removal and R multigaussian fit is indicated using the same colour and label conventions 
of Fig. 6.9. There is 1 blueshifted TAC (#6, blue colour) observed both in metallic and 
hydrogen lines, note that Avhi » Avpen- There are also 2 metallic TACs (@5, @7, 
green) observed only in metallic lines. Event @7 can be clearly identified by visual 
inspection as a relatively narrow peak in the Feu lines. In events #6 and @7 <v> does 
not exactly correspond with the local minima in the Fe u spectrum because the whole R 
profile is used in the multigaussian fits. (This figure is available in color in electronic 
form). 
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Table 6.5: Column density estimates (according to Eq. 6.1) of ions excited to the energy 
ground level of B aimer and Fe n 42 triplet lines for events WW Vul #7 (metallic) and 
XY Per #8. 





WW Vul #7 


XY Per #8 


Line 


Na (cm-2) 


Na (cm-2) 


D 

Fe n 4924A 


> 1.1 • 10^"^ 


> 4.0 • 10^^ 


FenSOlSA 


> 8.1 • 10^^ 


> 2.1 • 10^^ 


Feii5169A 


> 7.8 ■ 10^^ 


>4.4- 10^^ 




< 2.9- W 


> 1.3-1013 


Hy 


< 1.2-1012 


> 4.2 - 1013 


H(5 


< 3.3 • 10^2 


> 9.6 • 1013 






> 2.9 • lOi^ 



column densities causing the absorptions can be estimated according to the following 
formula (Spitzer, 1978), which becomes exact when the gas is optically thin: 



4n€o MeC^WA 



(6.1) 



Where Na is the column density of ions in the ground state of the line, e is the 
elementary charge, is the electron mass, eo is the permeability of free space, c is 
the speed of light, / is the oscillator strength of the line, A is the wavelength of the 
line and is the equivalent width of the absorption component. The column densities 
obtained are similar, within an order of magnitude, for all the metallic events. An upper 
limit on the column density of the H i atoms in the Balmer energy level (n=2) can be 
estimated if we assume that TACs with an intensity lower than 3 times the noise level 
of the spectra cannot be detected. Table 6.5 gives the estimated values for the Fen 
lines and some Balmer lines in the case of event WW Vul #7. Event #8 in XY Per has 
extremely weak saturated H i features and strong metallic components (H i line strengths 
are lower than or just comparable to the Fen lines); in this sense, this event is the more 
similar to the metallic ones among all detected events with H i and metallic components. 
Table 6.5 gives lower limits to the column densities, estimated from the Balmer and 
Fe n lines of event #8 in XY Per. The values of Table 6.5 nicely shows the difference 
between the estimated column densities of ions excited to the energy ground level of 
the Feu 42 triplet and Balmer lines in both type of events, and it could point out to a 
fundamental difference on the nature and origin of the gas from which the absorptions 
rise. However, a definitive statement on this issue requires a deep analysis of the metallic 
events detected in WW Vul in order to estimate chemical abundances, or abundance 
ratios among elements. Such analysis, which is beyond the scope of this paper, needs a 
NLTE treatment of the spectra; such treatment would produce model dependent results, 
since a previous knowledge of physical quantities, as for example volume densities of 
the gaseous clumps and electron temperature, are required. Such quantities cannot be 
estimated from our data in a confident way. 
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6.5 Concluding remarks 

We have analyzed optical high resolution spectra of the HAe stars BF Ori, SV Cep, WW 
Vul and XY Per. These spectra monitor the stars on time scales of months and days and, 
as in the case of the previously studied UX Ori (Paper I), they provide observational 
constraints, which should be considered for any realistic scenario of the gaseous circum- 
stellar disks around intermediate-mass PMS stars. Our results and conclusions can be 
summarized as follows: 

1. The gaseous circumstellar environment of these stars is very complex and active. 
The spectra always show circumstellar line absorptions with remarkable variations 
in their strength and dynamical properties. 

2. Variable absorption features are, in most cases, detected simultaneously in hydro- 
gen and in many metallic lines with similar velocities. In each case, there are 
several kinematic components in each line, both blue-shifted and red-shifted with 
respect to the systemic velocity, denoting the simultaneous presence of infalling 
and outflowing gas. We attribute the variable features detected in both Balmer and 
metallic lines to gaseous clumps of solar-like composition, evolving dynamically 
in the circumstellar disks of these objects. In this respect, the disks around the stars 
studied in this paper are similar to the UX Ori disk. Following the conclusions of 
Paper I we suggest that these clumps and their dynamical evolution should be in- 
vestigated in the context of detailed magnetospheric accretion models, similar to 
those of T Tauri stars. 

3. The star WW Vul is peculiar and behaves differently from the other stars studied 
in this paper and also from UX Ori. It is the only star that shows, in addition 
to events seen both in metallic and hydrogen lines, similar to those observed in 
the other stars, also transient absorption components in metallic lines that do not 
apparently have any obvious counterpart in the hydrogen lines. This result, taken 
at its face-value, would indicate the presence of a metal-rich gas component in 
the environment of WW Vul, possibly related to the evaporation of soUd bodies. 
However, any such conclusion is premature. We think that a series of optical spec- 
tra with better time resolution (hours) and longer monitoring (up to around seven 
days), spectra in the far UV range - to analyze Lyman and metallic resonance lines 
- and detailed NLTE models of different CS gas environments are essential for 
further progress and to provide clues on the origin of these apparently metal-rich 
events, in terms of their appearance/disappearance statistics, dynamics, metallicity 
and nature. 
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Table 6.6: Identified Transient Absorption Components (TACs) in BF Ori. Column 
1 gives the corresponding Balmer or metallic line. Col. 2 gives the Julian Date 

(-2450000), Col. 3 represents the event assigned to the particular absorptions (see text 
Sect. 6.3.2), Cols. 4 to 6 give the parameters of each transient absorption as estimated 
from the multigaussian fit of the normalized residual absorption R: v, radial velocity, 
Av, FWHM, and i?max> the strength of the absorption (peak of the gaussian). "0" in Col. 
3 corresponds to the narrow IS absorptions, while "-" means that the absorption is not 
associated with a particular event. 



Line 


JD 


Event 


v(kms ') 


AvCkms"') 


«max 


Line 


JD 


Event 


V (km s ^ ) 


Av(kms-') 


^max 


H/S 


1112.6324 


[ 


-8 


184 


0.95 


Feu5169A 


1209.5542 


2 


50 


199 


0.48 


Hr 


1112.6324 




7 


183 


0.95 


Ti u 4444A 


1209.5542 


2 


38 


185 


0.14 


H6 


1112.6324 




4 


157 


0.93 


Tiu4572A 


1209.5542 


2 


52 


189 


0.13 


»( 


1 1 12.6324 




11 


120 


0.79 


Scn4247A 


1209.5542 


2 


21 


75 


0.22 


NaiD2 


1 1 12.6324 




-1 


66 


0.83 


NaiD2 


1209.5542 





1 


31 


0.58 


NaiDl 


1112.6324 




-2 


60 


0.75 


NaiDl 


1209.5542 





1 


22 


0.51 


Fen4924A 


1112.6324 




3 


83 


0.64 


Feu4924A 


1209.5542 





4 


20 


0.09 


FenSOlsA 


1112.6324 




5 


92 


0.71 


Fen5018A 


1209.5542 





4 


18 


0.11 


Fen5169A 


1112.6324 




2 


93 


0.76 


Fen5169A 


1209.5542 





3 


24 


0.12 


Tiii4444A 


1112.6324 




-2 


61 


0.33 


CanK 


1209.5542 




-95 


136 


0.33 


Tiii4572A 


1112.6324 




-3 


58 


0.35 


CanK 


1209.5542 




96 


172 


0.71 


Sen 4247A 


1112.6324 




-1 


66 


0.30 




1210.4571 


2 


38 


155 


0.95 


Cai4227A 


1112.6324 




3 


62 


0.23 


Hr 


1210.4571 


2 


25 


185 


0.97 


Hf 


1112.6324 




30 


132 


0.75 


H5 


1210.4571 


2 


29 


155 


0.95 


CanK 


1112.6324 




-52 


70 


0.54 


He 


1210.4571 


2 


32 


130 


0.79 


CanK 


1112.6324 




92 


172 


0.61 


Hi 


1210.4571 


2 


23 


145 


0.79 


CanH 


1112.6324 




-40 


53 


0.71 


NaiD2 


1210.4571 


2 


13 


83 


0.47 




1113.6515 




27 


194 


0.86 


NaiDl 


1210.4571 


2 


16 


78 


0.32 


Hr 


1113.6515 




33 


189 


0.86 


Feu4924A 


1210.4571 


2 


23 


110 


0.52 


HiS 


1113.6515 




33 


195 


0.79 


Feu5018A 


1210.4571 


2 


23 


115 


0.59 




1113.6515 




39 


139 


0.65 


Feu5169A 


1210.4571 


2 


22 


118 


0.65 


NaiD2 


1113.6515 




6 


57 


0.80 


Tiu4444A 


1210.4571 


2 


20 


87 


0.22 


NaiDl 


1113.6515 




5 


51 


0.69 


Tiu4572A 


1210.4571 


2 


12 


117 


0.21 


Fen 4924A 


1113.6515 




15 


98 


0.54 


Seu4247A 


1210.4571 




14 


98 


0.21 


FcnSOlsA 


1113.6515 




20 


111 


0.62 


Hfi 


1210.4571 


3 


-52 


62 


0.32 


Fen5169A 


1113.6515 




19 


116 


0.66 


Hr 


1210.4571 


3 


-67 


94 


0.19 


Tin4444A 


1113.6515 




17 


76 


0.24 


US 


1210.4571 


3 


-61 


90 


0.27 


Tiii4572A 


1113.6515 




13 


70 


0.24 


He 


1210.4571 


3 


-65 


43 


0.15 


Sen 4247A 


1113.6515 




16 


83 


0.23 


HC 


1210.4571 


3 


-89 


63 


0.15 


Ca]4227A 


1113.6515 




13 


51 


0.15 


Can K 


1210.4571 


3 


-77 


92 


0.47 


He 


1113.6515 




73 


172 


0.63 


CanH 


1210.4571 


3 


-67 


93 


0.47 


CanK 


1113.6515 




-74 


134 


0.28 


Fen4924A 


1210.4571 


3 


-57 


77 


0.08 


CanK 


1113.6515 




96 


167 


0.59 


Fen5018A 


1210.4571 


3 


-53 


55 


0.07 


H/J 


1209.5542 


2 


40 


218 


0.83 


Feu5169A 


1210.4571 


3 


-54 


48 


0.08 


Hr 


1209.5542 


2 


40 


229 


0.88 


Tiii4444A 


1210.4571 


3 


-42 


46 


0.09 


Hi 


1209.5542 


2 


40 


217 


0.82 


Tin4572A 


1210.4571 


3 


-34 


25 


0.04 


He 


1209.5542 


2 


58 


176 


0.60 


NaiD2 


1210.4571 





-1 


15 


0.41 




1209.5542 


2 


58 


184 


0.60 


NaiDl 


1210.4571 








14 


0.47 


NaiD2 


1209.5542 


2 


55 


152 


0.21 


Fen4924A 


1210.4571 





3 


10 


0.05 


NaiDl 


1209.5542 


2 


38 


135 


0.18 


Fen5018A 


1210.4571 





2 


13 


0.06 


Fen 4924 A 


1209.5542 


2 


46 


188 


0.36 


Fen5169A 


1210.4571 





2 


13 


0.06 


FenSOlsA 


1209.5542 


2 


54 


196 


0.44 


CanK 


1210.4571 




78 


146 


0.69 
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Table 6.7: Identified Transient Absorption 
for Table 6.6. 







JD 




V (km s ' ) 


Av(kms"') 




H/? 
np 


950.6668 


— 


35 


52 


0.58 




950.6668 




40 


45 


0.35 




950.6668 




40 


46 


0.33 


He 


950.6668 




31 


50 


0.38 




950.6668 




40 


80 


0.28 


CanK 


950.6668 


J 


25 


27 


0.42 


CanH 


950.6668 


I 


23 


16 


0.38 


NaiD2 


950.6668 


1 


26 


14 


0.25 


Fen 4924A 


950.6668 


I 


23 


25 


0.15 


Fen 5018A 


950.6668 




25 


27 


0.22 


Fen 5169A 


950.6668 




24 


28 


0.25 


np 


950.6668 


2 


-17 


73 


0.93 




950.6668 


2 


-11 


84 


0.83 


US 


950.6668 


2 


-13 


79 


0.67 


He 


950.6668 


2 


-17 


49 


0.42 


H<r 


950.6668 


2 


-22 


74 


0.35 


Ca nK 


950.6668 


2 


-8 


28 


0.73 


Ca nH 


950.6668 


2 


-6 


24 


0.66 


Fe n 4924A 


950.6668 


2 


-11 


17 


0.10 


Fen 5018A 


950.6668 


2 


-9 


20 


0.15 


Fe n 5 169 A 


950 6668 


2 


-10 


16 


20 


Na I D2 


950 6668 


Q 




16 


1 04 


Nai Dl 


950 6668 







13 


1 05 


H/? 


1025.6260 


3 


99 


276 


0.34 


Hv 


1025.6260 


3 


78 


217 


0.48 


H(5 


1025.6260 


3 


76 


219 


0.43 




1025.6260 


3 


77 


210 


0.31 


Can K 


1025.6260 


3 


40 


193 


0.51 


NaiD2 


1025.6260 


3 


71 


162 


0.13 


NaiDl 


1025.6260 


3 


68 


199 


0.10 


Fen 4924 A 


1025.6260 


3 


51 


197 


0.09 


Fen 5018A 


1025.6260 


3 


70 


245 


0.12 


Fe n 5 169 A 


1025.6260 


3 


60 


207 


0.12 


tip 


1095 6960 


4 


-6 


147 


0.67 


Hr 


1025.6260 


4 


-21 


139 


0.54 


w 


1025.6260 


4 


-16 


133 


0.48 


He 


1025.6260 


4 


15 


231 


0.67 




1025.6260 


4 


-10 


134 


0.40 


CanK 


1025.6260 


4 





29 


0.41 


CanH 


1025.6260 


4 


-2 


19 


0.66 


Fen 4924 A 


1025.6260 


4 


12 


33 


0.09 


Fen5018A 


1025.6260 


4 


12 


39 


0.14 


Fen5169A 


1025.6260 


4 


10 


39 


0.15 


NaiD2 


1025.6260 








17 


0.95 


NaiDl 


1025.6260 





-0 


14 


0.96 


US 


1026.6684 


3 


80 


202 


0.75 


Hr 


1026.6684 


3 


72 


190 


0.76 




1026.6684 


3 


69 


164 


0.70 


He 


1026.6684 


3 


87 


116 


0.56 


H( 


1026.6684 


3 


66 


152 


0.55 


Ca uK 


1026.6684 


3 


75 


115 


0.53 


NaiD2 


1026.6684 


3 


82 


130 


0.16 


NaiDl 


1026.6684 


3 


75 


127 


0.12 


Fen4924A 


1026.6684 


3 


84 


87 


0.11 


Fen5018A 


1026.6684 


3 


87 


95 


0.16 


Fen5169A 


1026.6684 


3 


86 


83 


0.17 


H8 


1026.6684 


4 


-8 


82 


0.59 


Hr 


1026.6684 


4 


-11 


68 


0.42 


IK 


1026.6684 


4 


-13 


60 


0.35 


He 


1026.6684 


4 


-18 


112 


0.59 




1026.6684 


4 


-10 


53 


0.21 



Components (TACs) in SV Cep. Details as 



Line 


JD 


Event 


V (km s ^ ) 


Av (km s ' ) 


^max 


CanK 


1026.6684 


4 


-6 


46 


0.49 


CanH 


1026.6684 


4 


-22 


20 


0.32 


Fen4924A 


1026.6684 


4 


1 


40 


0.13 


Fen5018A 


1026.6684 


4 


-0 


43 


0.19 


Fen5169A 


1026.6684 


4 





42 


0.23 


CanK 


1026.6684 





-3 


12 


0.29 


CanH 


1026.6684 





-2 


17 


0.85 


NaiD2 


1026.6684 








18 


0.97 


NaiDl 


1026.6684 





-1 


14 


0.99 


H/J 


1113.4730 


5 


138 


38 


0.04 


Hr 


1113.4730 


5 


113 


64 


0.12 


W 


1113.4730 


5 


114 


81 


0.12 


He 


1113.4730 


5 


97 


41 


0.13 


HC 


1113.4730 


5 


85 


20 


0.13 


Can K 


1113.4730 


5 


89 


85 


0.11 


Feii4924A 


1113.4730 


5 


87 


63 


0.04 


Fen5018A 


1113.4730 


5 


90 


292 


0.05 


Fen5169A 


1113.4730 


5 


67 


173 


0.04 


US 


1113.4730 


6 


-28 


125 


0.99 


Hr 


1113.4730 


6 


-21 


117 


0.94 


H6 


1113.4730 


6 


-22 


106 


0.86 


He 


1113.4730 


6 


-13 


85 


0.65 




1113.4730 


6 


-18 


95 


0.52 


CanK 


1113.4730 


6 


-19 


77 


0.56 


CanH 


1113.4730 


6 


-11 


65 


0.53 


NaiD2 


1113.4730 


6 


-26 


12 


0.14 


NaiDl 


1113.4730 


6 


-27 


30 


0.07 


Feii4924A 


1113.4730 


6 


-15 


78 


0.15 


Fen5018A 


1113.4730 


6 


-14 


70 


0.17 


Fen5169A 


1113.4730 


6 


-23 


77 


0.18 


CanK 


1113.4730 





-3 


15 


0.44 


CanH 


1113.4730 





-1 


12 


0.50 


NaiD2 


1113.4730 





-1 


16 


1.04 


NaiDl 


1113.4730 





-1 


14 


1.03 


Fen4924A 


1113.4730 





2 


13 


0.09 


Fen5018A 


1113.4730 





-1 


14 


0.11 


Feii5169A 


1113.4730 





-2 


23 


0.16 


H/J 


1209.3372 


7 


198 


349 


0.33 


Hr 


1209.3372 


7 


164 


315 


0.32 


m 


1209.3372 


7 


134 


334 


0.30 


He 


1209.3372 


7 


141 


298 


0.29 


H^ 


1209.3372 


7 


143 


313 


0.24 


CanK 


1209.3372 


7 


172 


269 


0.19 


Fen 4924 A 


1209.3372 


7 


136 


324 


0.05 


FenSOlsA 


1209.3372 


7 


156 


253 


0.07 


Fen5169A 


1209.3372 


7 


167 


278 


0.06 


H/3 


1209.3372 


8 


-4 


99 


0.97 


Hr 


1209.3372 


8 


3 


111 


0.77 


w 


1209.3372 


8 


6 


110 


0.60 


He 


1209.3372 


8 


6 


100 


0.51 


H^- 


1209.3372 


8 


9 


99 


0.46 


CanK 


1209.3372 


8 


12 


115 


0.40 


NaiD2 


1209.3372 


8 


40 


96 


0.07 


NaiDl 


1209.3372 


8 


33 


51 


0.06 


Fen4924A 


1209.3372 


8 


19 


76 


0.09 


Fen5018A 


1209.3372 


8 


16 


78 


0.13 


Fen5169A 


1209.3372 


8 


17 


74 


0.15 


CanK 


1209.3372 





-1 


18 


0.43 


CanH 


1209.3372 





1 


22 


0.77 


Na 1 D2 


1 211'),.; .572 





1 


211 


11.87 


Na 1 L) 1 


1209..i.i72 





1 


l.-i 


0.9J 
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Table 6.8: Identified Transient Absorption 
for Table 6.6. 



Line 


JD 


Event 


v(kms ^) 


Av (km s ' ) 






950.6176 




131 


95 


0.10 


Hy 


950.6176 




71 


162 


0.41 


H<5 


950.6176 


\ 


58 


174 


044 


He 


950.6176 


\ 


59 


164 


040 




950.6176 




69 


168 


0.38 


CanK 


950.6176 




74 


170 


0.34 


Fen 4924 A 


950.6176 




40 


186 


0.16 


Fe II 5018A 


950.6176 




45 


203 


0.20 


Fe II 5 169 A 


950.6176 




44 


207 


0.19 






2 


-13 


36 


17 


feu JUiiSA 


Qcn fil If, 
yju.oi lO 




-12 


41 


0.22 


Fen 5169A 


950.6176 


2 


-12 


41 


0.26 




yju.oi /o 


3 


-25 


1 10 


1.27 


Hy 


yju.o I/O 












^ JU.U 1 / u 


3 


-34 


106 


0.73 


He 


950.6176 


3 


-30 


87 


0.61 




950.6176 


3 


-31 


104 


0.62 


Ca II K 


Qsn fil 76 


3 


-28 


89 


61 


Ca II H 


950.6176 


3 


-11 


58 


0.98 


Nai D2 


950.6176 


3 


-34 


91 


0.16 


Nai Dl 


950.6176 


3 


-22 


59 


0.12 


np 


950.6176 


4 


-104 


58 


0.63 


Ht 
ny 


950.6176 


4 


-97 


55 


0.43 


Hi5 


950.6176 


4 


-90 


55 


0.39 




950.6176 


4 


-82 


49 


0.25 


Ca II K 


950.6176 


4 


-86 


42 


0.31 


Ca II H 


950.6176 


4 


-74 


55 


0.63 


Fen 4924 A 


950.6176 


4 


-67 


45 


0.09 


Fen 501 8A 


950.6176 


4 


-68 


46 


0.12 


Feu 5169A 


950.6176 


4 


-66 


41 


0.13 


Ca II K 


950.6176 





_4 


23 


0.21 


Ca II H 


950.6176 


Q 


_[ 


7 


0.09 


Nai D2 


950.6176 





-2 


22 


0.86 


Nai Dl 


950.6176 





_1 


16 


0.80 


li0 


951.6232 




93 


75 


0.33 


Hy 


951.6232 




72 


109 


0.53 


HS 


951.6232 




42 


136 


0.68 


He 


951.6232 




70 


91 


0.34 


CanK 


951.6232 




30 


142 


0.57 


NaiD2 


951.6232 




59 


62 


0.14 


NaiDl 


951.6232 


\ 


67 


43 


0.12 


Fe II 4924A 


951.6232 




45 


112 


0.18 


Feu5018A 


951.6232 




44 


111 


0.26 


Fen 5169A 


951.6232 


^ 


46 


113 


0.27 


CanK 


951.6232 


2 


-21 


53 


0.36 


Ca II H 


951.6232 


2 


-12 


46 


0.74 


NaiD2 


951.6232 


2 


-14 


59 


0.50 


NaiDl 


951.6232 


2 


-10 


58 


0.36 


Fen4924A 


951.6232 


2 


-20 


42 


0.12 


Fen 501 SA 


951.6232 


2 


-20 


52 


0.21 


Fe II 5 1 69A 


951.6232 


2 


-23 


66 


0.33 


np 


951.6232 


3 


-30 


134 


1.10 


Hv 
ny 


951.6232 


3 


-39 


136 


0.96 


US 


951.6232 


3 


-53 


107 


0.70 


He 


951.6232 


3 


-21 


129 


0.79 




951.6232 


3 


-49 


52 


0.21 


Ca II K 


951.6232 


3 


-64 


41 


0.31 


CauH 


951.6232 


3 


-57 


50 


0.57 


Fen4924A 


951.6232 


3 


-40 


135 


0.13 


Fen 501 sA 


951.6232 


3 


-63 


127 


0.12 


Fe II 5169A 


951.6232 


3 


-91 


84 


0.10 


np 


951.6232 


4 


-174 


63 


0.14 


Hy 


951.6232 


4 


-170 


81 


0.19 


HS 


951.6232 


4 


-156 


78 


0.18 


H^ 


951.6232 


4 


-137 


119 


0.10 


CanK 


951.6232 


4 


-122 


113 


0.17 


CanH 


951.6232 


4 


-139 


70 


0.10 


CanK 


951.6232 





-3 


19 


0.14 


Ca II H 


951.6232 





-4 


9 


0.07 


NaiD2 


951.6232 





-0 


15 


0.54 


NaiDl 


951.6232 








12 


0.62 


Fen4924A 


951.6232 





11 


11 


0.04 


Fen5018A 


951.6232 





11 


8 


0.06 


Fen5169A 


951.6232 





8 


8 


0.03 


H< 


951.6232 




10 


163 


0.75 


CanK 


1023.5186 


5 


41 


106 


0.46 


CanH 


1023.5186 


5 


28 


28 


0.28 


Fen 4924 A 


1023.5186 


5 


37 


115 


0.08 


Feii 5018A 


1023.5186 


5 


50 


132 


0.11 


Feii5169A 


1023.5186 


5 


40 


127 


0.11 


Ti u 4444A 


1023.5186 


5 


63 


113 


0.02 




1023.5186 


6 


-11 


157 


1.31 


Hy 


1023.5186 


6 


-9 


155 


1.09 


H5 


1023.5186 


6 


-9 


131 


0.99 



Components (TACs) in WW Vul. Details as 



Line 


JD 


Event 


V (km s ^ ) 


Av (km s ' ) 


^max 


He 


1023.5186 


6 


-3 


99 


0.98 


Hf 


1023.5186 


6 


-9 


101 


0.85 


CanK 


1023.5186 


6 


-11 


62 


0.66 


CanH 


1023.5186 


6 


-6 


46 


1.15 


Fen 4924 A 


1023.5186 


6 


-6 


34 


0.32 


FenSOlsA 


1023.5186 


6 


-6 


40 


0.41 


Fen5169A 


1023.5186 


6 


-5 


39 


0.45 


Tiu4444A 


1023.5186 


6 


-6 


36 


0.05 


Tiu4572A 


1023.5186 


6 


-7 


37 


0.05 


CanK 


1023.5186 


7 


-46 


28 


0.41 


CanH 


1023.5186 


7 


-41 


31 


0.85 


NaiD2 


1023.5186 


7 


-30 


25 


0.50 


NaiDl 


1023.5186 


7 


-23 


34 


0.35 


Fen 4924 A 


1023.5186 


7 


-35 


24 


0.30 


FenSOlsA 


1023.5186 


7 


-35 


23 


0.33 


Fen5169A 


1023.5186 


7 


-35 


24 


0.40 


Tin4444A 


1023.5186 


7 


-34 


21 


0.04 


Ti II 4572A 


1023.5186 


7 


-34 


22 


0.04 


NaiD2 


1023.5186 





-1 


22 


0.93 


NaiDl 


1023.5186 








14 


0.86 


Fen4924A 


1112.3689 


8 


86 


162 


0.09 


FenSOlsA 


1112.3689 


8 


96 


195 


0.13 


Fen5169A 


1112.3689 


8 


113 


248 


0.11 


Tin4444A 


1112.3689 


8 


91 


163 


0.03 


Tin4572A 


1112.3689 


8 


78 


224 


0.03 




1112.3689 


9 


7 


156 


0.98 


Hy 


1112.3689 


9 


17 


159 


0.97 


H5 


1112.3689 


9 


17 


152 


0.99 


He 


1112.3689 


9 


11 


174 


1.04 


Hf 


1112.3689 


9 


24 


145 


0.93 


CanK 


1112.3689 


9 


23 


148 


0.78 


NaiD2 


1112.3689 


9 


21 


83 


0.59 


NaiDl 


1112.3689 


9 


22 


72 


0.45 


Fcii4924A 


1112.3689 


9 


18 


77 


0.32 


Fen5018A 


1112.3689 


9 


17 


82 


0.39 


Fen5169A 


1112.3689 


9 


18 


85 


0.45 


Fer4046A 


1112.3689 


9 


32 


76 


0.03 


Tin4444A 


1112.3689 


9 


14 


57 


0.08 


Tiii4572A 


1112.3689 


9 


16 


62 


0.09 


Sen 4247A 


1112.3689 


9 


35 


79 


0.08 


Ca I 4227A 


1112.3689 


9 


17 


66 


0.03 


CanK 


1112.3689 





-7 


25 


0.16 


CanH 


1112.3689 





-5 


56 


0.89 


NaiD2 


1112.3689 





-2 


14 


0.51 


NaiDl 


1112.3689 





-1 


13 


0.62 


Fcu4924A 


1112.3689 








13 


0.11 


Feu5018A 


1112.3689 





1 


12 


0.12 


Fen5169A 


1112.3689 





-0 


13 


0.15 


Scn4247A 


1112.3689 





2 


31 


0.07 


Hy 


1113.3958 


8 


129 


248 


0.42 


H<5 


1113.3958 


8 


106 


233 


0.53 


He 


1113.3958 


8 


85 


257 


0.62 


H< 


1113.3958 


8 


106 


220 


0.47 


CanK 


1113.3958 


8 


108 


239 


0.46 


NaiD2 


1113.3958 


8 


82 


187 


0.14 


NaiDl 


1113.3958 


8 


129 


185 


0.11 


Fen 4924 A 


1113.3958 


8 


80 


236 


0.18 


Fen5018A 


1113.3958 


8 


90 


263 


0.24 


Fen5169A 


1113.3958 


8 


116 


291 


0.24 


Fei4046A 


1113.3958 


8 


127 


116 


0.04 


Tiu4444A 


1113.3958 


8 


74 


231 


0.06 


Tiu4572A 


1113.3958 


8 


107 


253 


0.06 


Sen 4247 A 


1113.3958 


8 


123 


140 


0.06 


Ca 1 4227 A 


1113.3958 


8 


134 


109 


0.03 




1113.3958 


9 


-4 


124 


0.96 


Hy 


1113.3958 


9 


-4 


121 


0.73 




1113.3958 


9 


-11 


109 


0.66 


He 


1113.3958 


9 


6 


68 


0.44 




1113.3958 


9 


-9 


102 


0.65 


CanK 


1113.3958 


9 


-4 


63 


0.49 


CanH 


1113.3958 


9 


-2 


58 


0.79 


NaiD2 


1113.3958 


9 


4 


48 


0.74 


NaiDl 


1113.3958 


9 


5 


43 


0.76 


Fen4924A 


1113.3958 


9 


4 


43 


0.34 


FeuSOlsA 


1113.3958 


9 


4 


47 


0.39 


Fcu5169A 


1113.3958 


9 


4 


49 


0.46 


Fei4046A 


1113.3958 


9 


3 


40 


0.06 


Tin4444A 


1113.3958 


9 


6 


36 


0.09 


Tin4572A 


1113.3958 


9 


1 


48 


0.09 


Sen 4247 A 


1113.3958 


9 


6 


55 


0.12 


Ca 1 4227 A 


1113.3958 


9 


6 


46 


0.10 


w 


1113.3958 




251 


149 


0.20 
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Table 6.9: Identified Transient Absorption Components (TACs) in XY Per. Details as 
for Table 6.6. 



Line 


JD 


Event 


V (km s ^ ) 


Av(kms ') 


/JniiLK 


Line 


JD 


Event 


V (km s ^ ) 


Av(kms"') 


/^max 


H/J 


1024.6728 


2 


-9 


171 


0.92 


Ca nK 


1026.7065 


2 


-66 


91 


0.55 


Hy 


1024.6728 


2 


-4 


171 


0.92 


CanH 


1026.7065 


2 


-51 


101 


0.55 


m 


1024.6728 


2 


-3 


163 


0.88 


NaiD2 


1026.7065 


2 


-36 


102 


0.26 


Hf 


1024.6728 


2 


2 


147 


0.72 


NaiDl 


1026.7065 


2 


-43 


64 


0.18 


Hf 


1024.6728 


2 


-5 


146 


0.70 


Fen4924A 


1026.7065 


2 


-52 


83 


0.37 


Cau K 


1024.5728 


2 


-7 


142 


0.77 


Fen5018A 


1026.7065 


2 


-56 


66 


0.43 


CanH 


1024.6728 


2 


-27 


104 


0.63 


Fen5169A 


1026.7065 


2 


-54 


82 


0.45 


Nai D2 


1024.6728 


2 


-26 


100 


0.16 


Tin4444A 


1026.7065 


2 


-52 


66 


0.09 


Nai Dl 


1024.6728 


2 


-19 


92 


0.10 


Tin4572A 


1026.7065 


2 


-53 


68 


0.10 


Fe n 4924A 


1024.6728 


2 


-15 


118 


0.25 


Sen 4247 A 


1026.7065 


2 


-54 


68 


0.06 


Fen5018A 


1024.6728 


2 


-14 


115 


0.32 




1026.7065 


4 


-3 


54 


0.42 


Fen5169A 


1024.6728 


2 


-16 


114 


0.35 


Hy 


1026.7065 


4 





41 


0.24 


Tin4444A 


1024.6728 


2 


2 


156 


0.06 


m 


1026.7065 


4 


2 


38 


0.19 


Tiii4572A 


1024.6728 


2 


-11 


144 


0.06 


He 


1026.7065 


4 


-4 


38 


0.21 


Sen 4247A 


1024.6728 


2 


-3 


133 


0.04 


Hi 


1026.7065 


4 


4 


42 


0.20 


CanK 


1024.6728 





-0 


18 


0.23 


Can K 


1026.7065 


4 


-8 


48 


0.63 


CanH 


1024.6728 





-1 


25 


0.44 


CanH 


1026.7065 


4 


-4 


33 


0.49 


NaiD2 


1024.6728 





-0 


19 


0.93 


Feu4924A 


1025.7065 


4 


-3 


26 


0.11 


NaiDl 


1024.6728 








17 


0.98 


Feu5018A 


1026.7065 


4 


-1 


46 


0.28 


Fen 4924 A 


1024.6728 





-5 


27 


0.06 


Fen5169A 


1026.7065 


4 


-5 


29 


0.15 


Fen5018A 


1024.6728 





-6 


23 


0.06 


Tin4444A 


1026.7065 


4 


-1 


33 


0.04 


Fen5169A 


1024.6728 





-4 


27 


0.07 


Tin4572A 


1026.7065 


4 


5 


50 


0.05 


Sen 4247A 


l()2l.(i72» 





-7 


9 


0.01 


Seu4247A 


1026.7065 


4 


12 


52 


0.03 


H/5 


1025.6948 


1 


70 


72 


0.37 


NaiD2 


1026.7065 





-0 


17 


0.89 


Hy 


1025.6948 


1 


74 


58 


0.16 


NaiDl 


1026.7065 





-0 


17 


1.02 


H,S 


1025.6948 


1 


66 


93 


0.21 


H/} 


1112.4978 


5 


131 


149 


0.27 


He 


1025.6948 


1 


60 


59 


0.23 


Hy 


1112.4978 


5 


83 


206 


0.41 


Fen 4924 A 


1025.6948 


1 


59 


66 


0.05 


HS 


1112.4978 


5 


118 


173 


0.30 


FenSOlsA 


1025.6948 


1 


73 


66 


0.05 


H( 


1112.4978 


5 


127 


182 


0.20 


Fen5169A 


1025.6948 


1 


65 


57 


0.04 


Can K 


1112.4978 


5 


119 


177 


0.29 


Tin4444A 


1025.6948 


1 


34 


113 


0.06 


Fc II 4924A 


1112.4978 


5 


147 


139 


0.07 


Tiii4572A 


1025.6948 


1 


44 


79 


0.04 


FcuSOlsA 


1112.4978 


5 


123 


190 


0.13 


Sen 4247A 


1025.6948 


1 


37 


88 


0.04 


Fcu5169A 


1112.4978 


5 


147 


218 


0.11 


H/J 


1025.6948 


2 


-10 


90 


0.90 


Ti u 4444A 


1112.4978 


5 


122 


209 


0.02 


Hy 


1025.6948 


2 


-22 


147 


0.96 


H/J 


1112.4978 


6 


30 


129 


0.77 


H<S 


1025.6948 


2 


-27 


133 


0.90 


Hy 


1112.4978 


6 


18 


136 


0.58 


He 


1025.6948 


2 


-15 


96 


0.63 


HS 


1112.4978 


6 


14 


142 


0.73 




1025.6948 


2 


-11 


141 


0.70 


He 


1112.4978 


6 


16 


189 


0.84 


CanK 


1025.6948 


2 


-20 


130 


0.78 


H^ 


1112.4978 


6 


17 


139 


0.59 


CanH 


1025.6948 


2 


-20 


94 


0.73 


Can K 


1112.4978 


6 


35 


104 


0.47 


Na]D2 


1025.5948 


2 


-28 


104 


0.24 


NaiD2 


1112.4978 


6 


40 


119 


0.15 


NaiDl 


1025.6948 


2 


-29 


89 


0.16 


NaiDl 


1112.4978 


6 


29 


132 


0.10 


Fen4924A 


1025.6948 


2 


-35 


115 


0.27 


Fen4924A 


1112.4978 


6 


26 


124 


0.25 


FenSOlsA 


1025.6948 


2 


-29 


117 


0.34 


Fen5018A 


1112.4978 


6 


22 


114 


0.26 


Fcn5169A 


1025.6948 


2 


-34 


113 


0.38 


Feu5169A 


1112.4978 


6 


18 


121 


0.28 


Ti II 4444A 


1025.6948 


2 


-44 


84 


0.05 


Fci4046A 


1112.4978 


5 


39 


100 


0.03 


Tiii4572A 


1025.6948 


2 


-42 


88 


0.06 


Ti n 4444A 


1112.4978 


5 


32 


117 


0.06 


Sen 4247A 


1025.6948 


2 


-42 


69 


0.03 


Tin4572A 


1112.4978 


5 


38 


139 


0.08 


H/J 


1025.6948 


3 


-88 


82 


0.65 


Sen 4247 A 


1112.4978 


6 


36 


137 


0.07 


Hy 


1025.5948 


3 


-89 


51 


0.17 


Ca 1 4227 A 


1112.4978 


6 


23 


66 


0.05 


as 


1025.5948 


3 


-85 


46 


0.14 


US 


1112.4978 


7 


-26 


57 


0.58 


Hi 


1025.6948 


3 


-68 


58 


0.14 


Hy 


1112.4978 


7 


-24 


44 


0.34 


CanK 


1025.6948 


3 


-94 


40 


0.17 


m 


1112.4978 


7 


-21 


37 


0.23 


CanH 


1025.6948 


3 


-85 


51 


0.28 


H( 


1112.4978 


7 


-16 


40 


0.15 


FenSOlsA 


1025.6948 


3 


-77 


28 


0.03 


CanK 


1112.4978 


7 


-34 


33 


0.53 


Fen5169A 


1025.6948 


3 


-79 


29 


0.04 


CanH 


1112.4978 


7 


-25 


43 


0.90 


CanK 


1025.6948 





1 


16 


0.27 


NaiD2 


1112.4978 


7 


-23 


19 


0.31 


CanH 


1025.6948 





1 


14 


0.32 


NaiDl 


1112.4978 


7 


-22 


18 


0.18 


NaiD2 


1025.6948 





-1 


19 


0.87 


Fen4924A 


1112.4978 


7 


-22 


24 


0.18 


NaiDl 


1025.6948 





-1 


18 


0.96 


Fcii5018A 


1112.4978 


7 


-22 


24 


0.24 


Fen 4924A 


1025.6948 





1 


23 


0.15 


Fcu5169A 


1112.4978 


7 


-22 


22 


0.27 


Fen5018A 


1025.6948 





-0 


20 


0.15 


Tin4444A 


1112.4978 


7 


-21 


20 


0.04 


Fen5169A 


1025.6948 





-1 


23 


0.18 


Tin4572A 


1112.4978 


7 


-24 


23 


0.04 


Tin4444A 


1025.6948 





1 


20 


0.04 


CanK 


1112.4978 





-2 


30 


0.48 


Ti II 4572A 


1025.6948 





1 


28 


0.03 


Ca nH 


1112.4978 





2 


20 


0.61 


Sen 4247A 


1025.6948 





2 


22 


0.02 


NaiD2 


1112.4978 





1 


19 


0.92 


H/J 


1026.7065 


1 


69 


75 


0.24 


NaiDl 


1112.4978 





1 


16 


0.97 


Hy 


1026.7065 




62 


81 


0.23 


US 


1113.5483 


5 


92 


53 


0.08 


US 


1026.7065 




55 


96 


0.24 


Hy 


1113.5483 


5 


76 


107 


0.26 


He 


1026.7065 




65 


59 


0.10 


HS 


1113.5483 


5 


75 


123 


0.28 


H( 


1026.7065 




51 


72 


0.20 


H<r 


1113.5483 


5 


82 


111 


0.21 


CanK 


1026.7065 




49 


27 


0.10 


CanK 


1113.5483 


5 


76 


51 


0.11 


Fen4924A 


1026.7065 




30 


76 


0.10 


Fen4924A 


1113.5483 


5 


85 


85 


0.12 


Fen5169A 


1026.7065 




24 


82 


0.12 


Fen5018A 


1113.5483 


5 


79 


136 


0.14 


Ti II 4444A 


1026.7065 




32 


113 


0.04 


Fen5169A 


1113.5483 


5 


87 


94 


0.14 


Ti 11 4572A 


1026.7065 




62 


68 


0.02 


Tiu4444A 


1113.5483 


5 


82 


115 


0.03 


Sen 4247 A 


1026.7065 




58 


51 


0.02 


Tin4572A 


1113.5483 


5 


94 


83 


0.03 


US 


1026.7065 


2 


-58 


124 


0.79 


US 


1113.5483 


7 


-16 


149 


0.99 


Hy 


1026.7065 


2 


-48 


117 


0.78 


Hy 


1113.5483 


7 


-23 


133 


0.93 


HS 


1026.7065 


2 


-52 


105 


0.76 


HS 


1113.5483 


7 


-25 


126 


0.89 


He 


1026.7065 


2 


-61 


51 


0.40 


He 


1113.5483 


7 


1 


133 


0.51 


Hf 


1026.7065 


2 


-51 


84 


0.57 


H^ 


1113.5483 


7 


-19 


117 


0.73 
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Table 6.9: (Cont.) 



Line 


JD 


Event 


v(kms ^) 


Av(kms ^) 


Rmax 


CanK 


1113.5483 


7 


-18 


107 


071 


CaiiH 


1113.5483 


7 


-22 


83 


0.81 


NaiD2 


1113.5483 


7 


-4 


64 


0.37 


NaiDl 


1113.5483 


7 


-3 
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Chapter 7 
Conclusions 



In this Chapter, the most relevant conclusions derived from this thesis are summarized. 
Most of them have already been given in Sections 5.6 and 6.5, included in Paper I (Chap- 
ter 5) and Paper II (Chapter 6), respectively. Some considerations about the spectral 
database used (Chapters 2 and 3) and the measurements of v sin i (Chapter 4) are also 
included. 

1. The excellent quality echelle spectra used in this thesis belong to a large database 
which includes a large number of objects extensively monitored. This unique set of 
data has allowed a very detailed study, without precedents, of the TACs in UXOR 
stars. All the spectra are reduced and available, under request, to all astronomers 
who ask for them. 

2. A systematic study of the projected rotational velocities, v sin i, for all the 49 stars 
in the sample observed with the UES echelle spectrograph has been performed. 
These measurements, apart from having a high scientific value on their own, have 
been of great help in the determination of the physical parameters for the stars in 
this thesis. 

3. Optical high resolution spectra of the HAe UXOR stars BF Ori, SV Cep, UX Ori, 
WW Vul and XY Per have been analysed. These spectra monitor the stars on time 
scales of months, days and hours. Significant activity in the CS disks is always 
present, which manifests itself in the continuous appearance and disappearance of 
absorption components detected in hydrogen and in many metallic lines. These 
components display a wide range of variability in radial velocity, velocity disper- 
sion and intensity. This activity is not related to substantial variations of the stellar 
photosphere. 

4. Variable absorption features are, in most cases, detected simultaneously in hydro- 
gen and in many metallic lines with similar velocities. In each case, there are 
several kinematic components in each line, both blue-shifted and red- shifted with 
respect to the systemic velocity, denoting the simultaneous presence of infalling 
and outflowing gas. We attribute the variable features detected in both B aimer and 
metallic lines to gaseous clumps of solar-like composition, evolving dynamically 
in the circumstellar disks of these objects. 



99 



100 



Chapter 7. Conclusions 



5. Blobs undergo accelerations/decelerations of the order of tenths of ms~^ and last 
for a few days. Detectable changes in the gas dynamics occur on a time scale of 
hours, but the intrinsic velocity dispersion of the blobs appears to remain rather 
constant. It has been found a positive correlation between the velocity of the 
clumps and their velocity dispersion. No noticeable differences are seen in the 
properties of the infalling and outflowing gas, although infalls generally display 
larger velocity and velocity dispersion. Strong indications of a negative correlation 
between the dispersion velocity and the absorption intensity, during the dynamical 
evolution of the gaseous blobs, have been found. 

6. The relative intensity of the absorptions for several line multiplets has been stud- 
ied. The diff"erent saturation levels presented by the multiplets reveal the diversity 
of physical conditions in the CS disks of the objects. Most of the gas clumps 
are detected both in hydrogen and metallic lines. The relative absorption strength 
of the transient absorptions have been analysed in the context of NLTE models. 
The results suggest that the gas abundances are similar to the solar metallicity, 
ruling out the evaporation of solid bodies as the general physical origin of the 
transient absorptions in UXOR stars. It is suggested that these clumps and their 
dynamical evolution should be investigated in the context of detailed magneto- 
spheric accretion models, similar to those of T Tauri stars. The models should be 
non-axi symmetric, in order to explain the observed diversity of temporal evolution 
behaviours. 

7. The star WW Vul is peculiar and behaves diff"erently from the other stars stud- 
ied in this thesis. It is the only star that shows, in addition to events seen both 
in metallic and hydrogen lines, similar to those observed in the other stars, also 
transient absorption components in metallic lines that do not apparently have any 
obvious counterpart in the hydrogen lines. This result, taken at its face- value, 
would indicate the presence of a metal-rich gas component in the environment of 
WW Vul, possibly related to the evaporation of solid bodies. However, any such 
conclusion is premature. We think that a series of optical spectra with better time 
resolution (hours) and longer monitoring (up to around seven days), spectra in the 
far UV range - to analyse Lyman and metallic resonance lines - and detailed NLTE 
models of diff"erent CS gas environments are essential for further progress and for 
providing clues on the origin of these apparently metal-rich events, in terms of 
their appearance/disappearance statistics, dynamics, metallicity and nature. 

8. Finally, the whole set of data and relations obtained in this thesis provide obser- 
vational constraints, which should be considered for any realistic scenario of the 
gaseous circumstellar disks around intermediate-mass PMS stars. 



Chapter 8 
Conclusiones 



En este capitulo se enumeran las conclusiones mas relevantes que se derivan de esta tesis. 

La mayor parte de ellas yahan sido expuestas en las secciones 5.6 y 6.5, correspondientes 
a los artfculos Paper 1 (capitulo 5) y Paper 11 (capitulo 6), respectivamente. Se incluyen 
ademas unas consideraciones acerca del conjunto de datos empleado (capitulos 2 y 3) y 
de las medidas de v sin i (capitulo 4). 

1. Los espectros echelle empleados en esta tesis forman parte de una base de datos 
de excelente calidad, amplia cobertura temporal y gran numero de objetos. Esto 
ha permitido efectuar un estudio de los TACs en estrellas UXOR con un grado de 
detalle sin precedentes. Todos los espectros estan reducidos y son accesibles, bajo 
peticion previa, a los investigadores que lo soliciten. 

2. Se ha Uevado a cabo un estudio sistematico de las velocidades de rotacion proyec- 
tadas segun la Imea de vision, v sin i, para el conjunto de las 49 estrellas de la 
muestra observadas con el espectrografo echelle UES. Estas medidas, ademas de 
tener un gran valor cientffico por si mismas, han facilitado la determinacion de los 
parametros ffsicos para las estrellas de esta tesis. 

3. Se han analizado los espectros visibles de alta resolucion de las siguientes estrellas 
HAe con comportamiento UXOR: BE Ori, SV Cep, UX Ori, WW Vul y XY Per. 
Se ha efectuado un seguimiento de estos objetos en escalas temporales de meses, 
dias y horas. Los espectros muestran una intensa actividad en los discos CS de 
estas estrellas, caracterizada por la aparicion y desaparicion de componentes de 
absorcion en Imeas de hidrogeno y metalicas. Estas componentes presentan una 
gran variabilidad en velocidad radial, anchura e intensidad. Esta actividad no esta 
relacionada con cambios apreciables en la fotosfera estelar. 

4. Generalmente, se observan en cada espectro varias componentes cinematicas en 
cada Imea, desplazadas tanto al rojo como al azul. Esto supone la existencia si- 
multanea de gas eyectado y en cafda hacia la estrella. En general, las componentes 
de absorcion detectadas simultaneamente en distintas Imeas, se pueden agrupar en 
funcion de la velocidad radial de las mismas. Se ha atribuido el origen y variabi- 
lidad de estas absorciones a la evolucion dinamica (aceleracion o desaceleracion) 
de condensaciones de gas, en el seno de los discos CS de las estrellas estudiadas. 
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5. Estas nubes de gas experimentan aceleraciones o desaceleraciones de valor en 
tomo a decimas de ms"^. La escala temporal en la cual aparecen y desaparecen 
estas condensaciones es de unos pocos dfas, aunque es posible detectar cambios de 
velocidad a las pocas horas. La anchura de las absorciones, relacionada con la dis- 
persion de velocidades de la nube de gas, permanece aproximadamente constante 
con el tiempo. Se ha encontrado, sin embargo, una correlacion positiva entre la 
velocidad de las condensaciones y su dispersion de velocidades. No se ban encon- 
trado grandes diferencias entre las propiedades del gas desplazado al rojo o al azul, 
aunque las mayores velocidades y dispersiones de velocidades se suelen presentar 
en el gas desplazado al rojo. Se ban encontrado fuertes indicios de que, durante la 
evolucion de una nube de gas, el aumento de la dispersion de velocidades conlleva 
una disminucion de la intensidad de las absorciones y viceversa. 

6. Se ha estudiado la intensidad relativa de las absorciones para Imeas pertenecientes 
a varios multipletes. Las distintas condiciones de saturacion de los mismos, mues- 
tran la diversidad de condiciones ffsicas presentes en los discos CS de los obje- 
tos estudiados. La gran mayoria de las condensaciones de gas se detectan tanto 
en Imeas metalicas como de hidrogeno. Un analisis detallado mediante modelos 
NLTE de este fenomeno, sugiere que las abundancias qufmicas del gas son simi- 
lares a la solar. Esto descarta la evaporacion de cuerpos solidos como el mecanis- 
mo general de produccion de las absorciones CS en estrellas UXOR. Todo apunta 
a que el analisis de los resultados de esta tesis se debe realizar en el contexto de 
modelos detallados de acrecion magnetosferica. Se requiere para ello el desarroUo 
de modelos teoricos con un grado de detalle similar al alcanzado para CTTSs, pero 
que tengan en cuenta la naturaleza no estacionaria de estos fenomenos. 

7. La estrella WW Vul es distinta del resto de objetos estudiados. En WW Vul, al 
igual que en las demas estrellas, se pueden encontrar eventos detectados tanto en 
Imeas metalicas como de hidrogeno. Pero ademas, en WW Vul se ban encontrado 
componentes metalicas sin contrapartida clara en las Imeas de hidrogeno. Este re- 
sultado podria estar relacionado con la presencia de nubes de gas ricas en metales, 
quizas originadas por la evaporacion de cuerpos solidos. A pesar de su interes, el 
hecho de encontrar estas componentes anomalas no es una prueba suficiente de la 
presencia de FEBs en WW Vul. Se necesita la combinacion de un fuerte esfuerzo 
observacional y teorico para poder arrojar luz sobre este aspecto. Por una parte, 
se deberia efectuar una ambiciosa campana de observaciones espectroscopicas en 
el visible, con una cobertura temporal amplia (unos siete dfas) y densa (espectros 
cada varias horas), complementada con espectros en el rango UV lejano, que per- 
mitiesen analizar las Imeas de resonancia metalicas y de Lyman. Por otra parte, 
se tendrian que desarroUar modelos NLTE representatives de diferentes entornos 
gaseosos CS y diversas condiciones ffsicas. De esta manera se podrfa efectuar un 
progreso real que proporcionara claves acerca del origen de estos eventos aparente- 
mente ricos en metales, tanto de su estadfstica de aparicion y desaparicion, como 
de su dinamica, metalicidad y procedencia. 

8. Por ultimo, el conjunto de todos los datos y relaciones obtenidos en esta tesis 
supone un conjunto de restricciones observacionales, que deben ser tenidas en 
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cuenta en los futures modelos de discos de gas circunestelares en estrellas PMS 
de masa intermedia. 
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